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Résumé. — On a étudié les spectres de R. M. N. du 3 chloro-propénil-benzène, de l’alcool cinna- 
mique et de l’anéthol en les analysant comme systèmes ABX, et ABX,. On utilise, pour les fré- 
quences et les intensités, des formules basées sur des fonctions formées par de simples produits. 
Le système ABX, est étudié en fonction du champ magnétique appliqué. On détermine les cons- 
tantes de couplage et le déplacement chimique pour chacun des composés indiqués. Les valeurs 
Han se confirment en les utilisant pour le calcul des spectres pour deux fréquences 
différentes. 


Abstract. — The proton magnetic resonance spectra of 3-chloro-propenyl-benzene, cynnamic 
alcohol and anethol are studied and analyzed as ABX, and ABX, systems. Formulas for fre- 
quencies and intensities are obtained using simple product functions as basic functions. The 
ABX,system is studied as a function of the applied field. The chemical shift and spin-spin cou- 
pling parameters are determined for the above substances and the values obtained are checked, 


calculating the spectra for two different spectrometer frequencies. 


I. Introduction. — Dans le présent travail on a 
déterminé l’ensemble des paramètres caractéris- 
tiques des spectres R. M. N. de haute résolution de 
plusieurs substances de structure : 


HeCLOH 
®, ® substitué 


] 


RC 4x 
HÉESSS ; 


Y CH, X 


I 


Ceci a pour but d’obtenir les valeurs des cons- 
tantes de couplage entre les différents protons, et 
les comparer à celles obtenues par d’autres auteurs 
[1, 2, 3, 4] pour des substances de structure sem- 
blable. On a analysé également le comportement 
des spectres ABX, pour différentes valeurs du 
champ magnétique ou, ce qui revient au même, à 
différentes fréquences ; nous nous proposons de 
vérifier ainsi que l’approximation X est maintenue 
dans les cas étudiés et que l’étude de la variation 
du groupe X, en fonction du champ extérieur per- 
met, comme dans le cas ABX, [4], de déterminer 
les signes relatifs des constantes de couplage Jax 
et Jex. On a vérifié que l’approximation ABX,; 
_ suit les mêmes lignes générales que dans le sys- 


tème ABX;, bien que SBX, < SBX,. Dans les cas 
présents l’interprétation du groupe AB est rendue 
difficile parce que Jar — das. 


II. Théorie des spectres ABX,. — Bien que le 
cas ABX, ait déjà été étudié, en employant des 
fonctions de symétrie (1.2), on a trouvé plus conve- 
nable le traitement suivant, basé sur des fonctions 
formées par de simples produits. Suivant la nomen- 
clature de B. P, S. [5], on considère le sys- 
tème ABX, comme un système de quatre spins 
dans lequel il y a deux protons équivalents indis- 
cernables (X,). Dans le tableau I sont indiquées 
les fonctions de base et les éléments diagonaux de 
matrice. Les seuls éléments non diagonaux sont : 
HS = Hs =1/2JisCommeilnya2pss 
mélange entre fonctions de base avec différents Zx, 
dans le groupe X,, les niveaux sont dégénérés ; 
cette dégénérescence n’influence pas les fréquences 
mais on doit en tenir compte en calculant les inten- 
sités relatives des raies. Il y a seulement mélange 
entre les états qui ont le même #; (composante 
dans la direction z du spin total) et le même Zx 
(spin total du groupe X:) [5]. Ceci simplifie la 
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TABLEAU I 
FONCTIONS DE BASE ET ÉLÉMENTS DIAGONAUX DE MATRICE POUR LE SYSTÈME ABX>2 
Re PORN I OUT Fr d. 
1 œoœ(occ) Étape Af2{va + vB) + 1/4(Ja8 + 2J4x + 2JBx) 
2 œo(ocfà) 1 4/2(va + v8) + 1/4(J48) 
3 œo(Bf) SE 0 A [2(va + vB) + 1/4(Ja8 — 2Jax — 2JBx) 
k aB(æo) 1 1 A2(va — v8) + 114(-— Jar + 2Jax — 2JBx) 
6) œ (x) 0 0 1 J2(va = VB) + 1/4 (— J'1B) 
6 aB(B6) —14  —1  1f2{va — v8) + 1/4(- Jar — 2J4x + 2J8x) 
7 Bal) äk 1 1/2(— va + VB) + 1[4(— Jap — 2J4x + 2J8x) 
8 Bor(æf) Gt A[2(-— va + v8) + 1/4(-— Jap) 
9 Ba(B6) 14  —1A4 1/2 va + v8) + 1/4(— Ja + 2Jax — 2JBx) 
10 BB(xx) 1 0 A[2(— va — vB) + 1/4(JaB — 2J4x — 2J8x) 
11 BB (of) .—17 1/2 va— vs) + 1/4(Ja8) 
12 BB(BB) 4 2  Af2— va — ve) + 114(Jas + 2ax + 2JBx) 
TABLEAU II | matrice, qui peut être divisée en trois sous-ma- 
ae ter trices du second ordre qui se résolvent de la façon 
RES A PR DAS usuelle. Le tableau IT donne les fonctions d’onde 
stationnaires et les énergies pour les états de mé- 
n Fe En lange et, en résultant, trois paramètres auxi- » | 


es LR liaires D;;. Les équations ont été simplifiées par le 2 
afloæ)] + sin 0, [Baluæ)] — 14/48 + Dar choix d’un système de coordonnées dans lequel : « 
aBlaB)] + sin 0, [Bx(aB)] —1/4Ja8 + D5e  Vx — 0. Le tableau [TT donne pour le système ABX, # 


6” cos 0, [xB(BB)] + sin 0, [B«(BB)] —1/4J41r8 + Ds les fréquences de transition et les intensités rela- 

7, — sin 6, [aBlux)] + cos 0, [Ba(ao)] —1/4Jap — Dar Vives. 

8” — sin 0, [«B(xB)] + cos 0, [BalxB)] — 1/4Jap — Dé AE x; f 
Mere PS cs HI. Variation de la partie X, du spectre avec le « 

Ÿ sin°0s [us(BB}E-F cos 6: 1Ba 186] LABS champ magnétique. —— L’examen du tableau III 

où : montre : que les transitions 13 et 14 du groupe X, 

Dan = 112 Jar? + [(va — vn) + (Jax — Jax)]? } 1? sont d'intensité maximum ; que la fréquence est . 

LOS PATTES Er indépendante des à, c’est-à-dire du champ magné- . 
5 — 1/2{ Jam? + (va — vs) } tique appliqué et que la différence entre les fré- 

Dis = 1/2 Jar? + [(va — v8) — (Jax — Jax)]? } /? quences est égale à |Jax + Jsx]. Comme dans le 


TABLEAU III 


FRÉQUENCES DE TRANSITION ET INTENSITÉS RELATIVES POUR LE SYSTÈME ABX; 


TRANSITION A ÉNerGies INTENSITÉ RELATIVE 
1) 1-7’ 1 [2(va + v) — 1/2J41p — 1/2(J18 + Jpx) — Dr 1 — sin 26, 
2) 2-8’ A [2(va + vB) — 1/2J48 — Dé 2(1 — sin 0;) 
3) 3-9” 1 [2(va + VB) —1/2J48 + 1[2(Jax + JBx) 1 — sin 20, 
4) k7-10 1 [2(va 7 YB) + 1/2J18 — 1 [2(J41x + JBx) ee Dy 4 + sin 206, 
5) 4-41 l'a ve) et I2Ta De 2(1 + sin 20,) 
6) 6’-12 1f2(va + vB) + 1/2Jap + 1/2(Jax + JBx) — D 4 + sin 20, 


7) 1-47 L2(va + VB) — 11/2748 —1/2(Jax + JBx) + Dy 1 + sin 26, 
8) 2-5’ 1 [2(va + v8) —1/2Ja8 + Du 2(1 + sin 20,) 
9) 3-6” 1f2{va + vB) — 11/2418 + 1/2(J4x + JBx) + Do 1 + sin 20, 
10) 7-10 1/2(va + v8) + 1/2Ja8 — 1/2(Jax + JBx) + Dar 1 — sin 20, 
11) 8/-11 A 2va + vs) 44/2038 © Da 2(1 — sin 20,) 
12) 97-12 12(va + V8) + 1/2J48 + 1/2(Jax + Jex) + Do 1 — sin 26, 
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TABLEAU III (suite) 
Transirions X:, 


13) 1-2 1/2(J4x + JBx) 2 

14) 2-3 1[2(Jax + Jax) pi 

45) 4-5’ DEUr 2 cos? ( 

16) 5-6 D;8 — Ds 2 cos? (0, — 
: 17) 7-8 — Dy + D; 2 cos? | 

18)8"-9° — D5s + Dso 2 cos? (0, — 

19) 10-11 — 1/2(Jax + Jex) 2 

20) 11-12 — 1/2(J4x + JBx) 2 


COMBINAISONS 


} 4-8" Das + Dis | 2 sin? (0, — 0;) 
22) 5-9’ DD 2 sin? ( 
15755: — Dis — Dy 2 sin? | 
24) 8-6’ — D59 — D; 2 sin? ( 
où : 


: J'1B J Jap 
PE A EN LE à D AUS SLT 
PMR OaB + (Jax + JBx) PER AB (es SaB — (Jax + JBx) 


cas ABX, [4] la position des autres raies du d t) 
groupe X, dépend de Sas, de la valeur du champ t (Me/s) 
appliqué et des signes relatifs de Jax et Jpx. Pour 

- des valeurs suffisamment grandes du champ (ou ” 
pour das > Jar et [Jax — Jex] > 0) le groupe X2 400 
se simplifie en se rapprochant du quadruplet donné 
par la théorie simple du premier ordre. Les raies 
mobiles se superposent et forment un doublet d’in- 
tensité presque semblable à celui donné par les raies 
fixes (fréquences 13 et 14) avec une séparation : 
égale à la valeur absolue de |[J4x — Jpx|. Quand pü 


Jax et Jsx sont de même signe 
(Jax + Jex > Jax — Jx), 
les raies fixes sont les raies extérieures de l’en- | | 


semble X,, quelle que soit la valeur de A5. Si Jax 


et Jpx ne sont pas de même signe, pour des valeurs 
assez élevées de 71,, les raies mobiles se placent à « 
l’extérieur des raies fixes. 

On peut préciser ce qui précède par les figures Îa il 


et 1b où l’on a calculé l’aspect du groupe X, en 
fonction de la fréquence de résonance des protons 
pour des champs différents et pour deux ensembles 
de paramètres : Pour la figure 1a, 
dan 0,25 pi pion. Jan = 14,0°c/s. ; 
Di — 7,0 c/s. ; Tor — + 0,0 cjs. 20 


Pour la figure 10, JBx — — 0,6 cjs. Il est \ 
intéressant de noter que, comme dans le cas 
ABX, [4], on observe ici une paire de raies qui 
coalescent en une seule pour une certaine valeur /\ 
du champ. Cette valeur de A, peut être obtenu en 0 


égalant à zéro les fréquences correspondant aux DESSERTE ST TER TE VS 
transitions 15 ou 17 et 16 ou 18. On obtient ainsi : 


1) Dir — Dis a 0 4 *% Ds = Ds 


RTC 


Fic. laet 1 b. — Aspect du groupe X, du spectre ABX, 

VJ as? + [va — vs) + (Jax — JBx)]|° pour différentes fréquences de mesure. Dans les cas con- 

; sSidérés #01n — D;25) D. Demers Jabt— 1e 01050 

— V Jap? + (va — vB)? Jax = 7,0cjs.; a) Jpx = + 0,6 c/s.; b) Jex = —0,6c/s. 
A A | 
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# x (va — vB) + (Jax — JBx) = + (va — 8) 


*'x €) Jax = JBx 
b) 2{(va — vs) = — (Jax — Jex) 
; { 
nYx VA — VB = 6a8.Ho = Te (Jax — JBx) 


——— (Jax— JBx) où : Hÿ — H, de croisement. 
AB 


+*+ Ho = 
2) Dig — Dés = 0 xx Dss Es D59 
#*x VJaB? + (va — vB)? 
= VJ ar? + [(va —8) 


#*x + (va — VB) = (va — vB) — (Jax — JBx) 


(JaAx 
x L a) Jax = JBx 
* * © b) 2{va — v8) = Jax — Jax 
1 
#*x (va — vB) = 648 Ho — 2 (Jax — JBx) 


1 


% 
Ho s4 261B 


(Jax — JBx) 


On observe ainsi un croisement dans les raies 
de la figure 1 pour une seulement des paires de 
raies mobiles, celui-ci se produisant pour les tran- 
sitions 16 et 18 quand Jax > Jex et pour 15 et 17 
Si Jax << Jx. Pour le 3 chloro-propénil-benzène on 
doit noter le croisement à 14,0 c/s, mais comme la 
largeur de$ raies est relativement grande, on 
l’observe déjà à 25 Mc/s (fig. 3a et 3b). 


IV. Interprétation des spectres. — Dans le cas 
où le spectre se rapproche du premier ordre (fig. 2) 

les paramètres nécessaires pour le calcul des fré- 
” quences d’un système ABX, s’obtiennent de façon 
similaire à celle développée pour le système ABX.. 
On suppose que la partie X, du spectre est inexis- 
tante, et les raies des quadruplets se rejoignent et 
laissent, au lieu de 12, quatre raies correspondant 
au spectre AB pur à partir desquelles on calcule J'18 
et Ôap. L'attribution des raies aux transitions 
correspondantes données par le tableau III n’offre 
pas de difficultés si J << à quels que soient les signes 
relatifs de Jax et JBx (fig. 2a et 2b), bien qu'il y 
ait une inversion de l’ordre des raies du groupe B. 
Quand J — à, l’aspect du spectre est bien distinct ; 
en effet quelques raies, par exemple celles corres- 
pondant aux transitions 3 et 12, sont d'intensité 
pratiquement nulle (fig. 2c et 2d). Dans les 
cas ABX, étudiés, les raies correspondant au 
spectre AB pur peuvent se déterminer à partir des 
transitions 2, 5, 8 et 11. À l’encontre du cas 
où J € 5, il n’y a pas ici d’ambiguité dans l’ordre 
des attributions du groupe AB, puisqu’un change- 
ment de signe de J4x ou Jsx ne modifie pas l’ordre 
des fréquences où apparaissent les raies corres- 
pondant aux transitions du tableau III (fig. 2c 
et 2d). Par conséquent, quand J = à, il n’y a pas 
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d’altération importante dans les intensités rela- … 
tives des raies du groupe B. Les valeurs des para- … 
mètres auxiliaires D;; s’obtiennent expérimenta- … 
lement à partir des différences des raies | 
2Dy7 = 17 —1—=10 —4,2D33 = 8 —2— 11 —5 54 
2D39 — 6 — 9 + Jar. À partir des expressions » 
algébriques correspondant aux D. on calcule la . 


1 2 0e4 RAC 
89 
POSER 

FiG. 2. — Partie AB du spectre du type ABX, : 
a) Sam — 1,00 p. p. m.; Jax et Jpx de même signe ; 
b) Sas — 1,00 p. p. m., Jax et JBx de signes opposés. 
c) ap — 0,25 p. p. m., Jax et Jpx de même signe ; 
d) Sam — 0,25 p. p. m., Jax et Jpx de signes opposés. 
Dans tous les cas : [Jax| — 7,0 c/s. ; [JBx| — 0,6 c/s. ; 


Halte 0e ps 


valeur de [Jax —— Jsx|. [ax —- Jex| étant déjà 
connu à partir du groupe X,, on peut déter- 
miner Jax et Jrx avec leurs signes relatifs (La 
théorie des ABX, ne permet pas de déterminer quel 
est le paramètre de signe négatif dans le cas de 
signes différents.) Avec cet ensemble de valeurs 
pour les paramètres, on caleule tout le spectre et en 
ajustant un peu les valeurs on obtient un bon 
accord avec les résultats expérimentaux. On n’a 
jamais observé de raie de combinaison. 


V. Partie expérimentale. — Les spectres ont été 
pris à deux fréquences distinctes : à 25 Mc/s avec 
un spectromètre modèle Trüb-Tauber et à 40 Me/s 
avec un Varian 4 300 en utilisant la technique et 
l’étalonnage habituels en haute résolution. Les sub- 


Noz 


stances étudiées sont des liquides commerciaux 
d’une grande pureté. On a même utilisé, dans cer- 
tains cas, des solutions dans S,C afin de diminuer 
la viscosité et obtenir ainsi des raies plus étroites, 
établissant par là que S,C agit comme un solvant 
inerte. 


VI. Application. — (a) 3 chloro-propénil-benzène : 
Les deux hydrogènes oléfiniques non équivalents 
donnent la partie AB du spectre tandis que les 
hydrogènes méthyléniques forment la partie X,. 
Les spectres obtenus à 40 et 25 Me/s sont donnés 
dans les figures 3a et 3b avec les spectres calculés. 


PE # 


3 ‘ 
VB YA 


Fic. 3. — Spectre du 3 chloro-propénil-benzène * 
a) 40 Mcjs ; b) 25 Mc/s. La raie marquée r est une bande 
latérale de rotation. Les spectres théoriques ont été cal- 
culés à partir des paramètres déterminés à une seule 
fréquence. On observe la constance de la séparation des 
composantes les plus importantes du groupe X, et la 
coalescence des raies mobiles. 


Comme les valeurs employées pour déduire les para- 
mètres sont des moyennes de plusieurs spectres, 1l 
n’est pas possible d’observer graphiquement une 
concordance parfaite entre les valeurs calculées et 
un seul des spectres employés. On observe qu’en 
accord avec le déplacement des raies du groupe X; 
en fonction du champ de la figure 1, le contour 
de X, se profile à 40 Mc/s comme un quadruplet 
et à 25 Mc/s comme un triplet. 

L'ensemble des paramètres proposés pour cette 
substance, basé sur la meileure concordance avec 
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le spectre moyen, est de Jin 112860) 
Jax = 7,0 cjs. se Jex = 1,0%c/825 
da = 0,28 p. p. m. ; 1/2(va + ve) — 2,22 p. p. m. ; 
VX — 0. ; 
b) Alcool cinnamique : La figure 4 représente les 
spectres calculés et observés à 40 Mc/s. Les para- 


ä 


Ÿ 


X 


F1G. 4. — Spectre de l’alcool cinnamique. 
Fréquence : 40 Meys. 


mètres proposés sont Tir = MAP ORCIE 
Jix= = 4, Tec ex 0 5S2c/ 
CAB = 0200 0p. Doi: 
A f2(va + vs) = 0,427 p. p. m. 


c) Anéthol : La figure 5 représente Les spectres 
calculés et observés de cette substance à 40 Meys. 


| 


16 


100S 
Ve Va mr ar es Ÿ 


X 


Fic. 5. — Spectre de l’anéthol. Fréquence : 40 Mc/s. On 
observe dans le groupe X;, l’existence d’un doublet 
d’écartement plus grande que le doublet principal, ce 
qui prouve l’existence d’une constante de couplage de 
signe négatif. 


Ce spectre est du type ABX, avec J + 3 et a été 
calculé à partir des formules données antérieu- 
rement [4]. La présence dans le groupe X, de raies 
extérieures au doublet d’intensité maximum, 
prouve l’existence de signes opposés pour Jax 
et /sx. Les paramètres proposés sont : 


TRE ER A OM MON A ÈS 
Jsx = — 17 cJs. ; din — 0,314 P. Pp. M. ; 
L2(va + vs) = 0,085 p. p. m. 
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TABLEAU IV VII. Discussion. — Le tableau IV montre les. 
; valeurs des constantes de couplage Jan pour les. 
RÉSULTATS OBTENUS PAR DIVERS AUTEURS b t ù die d ] È t t il Se t 
ei ee SO ne substances étudiées dans le présent travail, jointes 
SEMBLABLES A CELLES ÉTUDIÉES aux résultats d’autres auteurs [1, 2, 3, 4, 7] pour 
DANS LE PRÉSENT TRAVAIL des substances de structure similaire. Comme on a 
prouvé qu’un changement de signe de J;, ou Jr: 
M2 “V2  J23 H3  Jf3 ne modifie pas le spectre, on les a supposé positives, 
Le ce qui correspond à un travail théorique de » 
# 2 Karplus [6]. L’analyse spectrale des substances 
a) [2] Sc=c (1 131 SENTE étudiées indique que J,, et J,, sont de signe opposé. : 
m* con, EI 1841 Her Les résultats obtenus d’une part par Bothner- 
HR" APE os By [7], d'autre part par Sheppard et col. [2] pour 
b) [2] AS R]70 74 2 les substances de même type (analysées dans deux 
ts #20) © champs comme système ABC,) ont montré que J;3 
RAC A ouJ,, sont les constantes négatives ; nous avons 
cf Rs M Sr he donc admis également que, dans les cas étudiés ici, 
à Nu J:3 est négatif. Il est intéressant d'observer que, 
) Ne sauf pour l’alcool cinnamique, l’unique alcool dis- : 
CO0H CH, 4, ponible e cette structure, les valeurs des para- 
Dr aE 4? mètres sont tous d’un même ordre de grandeur qui 
e) {4 Je 149 +, est caractéristique de la structure commune à 
Gi Cha toutes ces substances. 
CH,0OC He Les spectres à 40 Mc/s ont été obtenus à l’Institut 
f) fe er 15,4 700 de Physique de l’Université d’'Uppsala (Suède) et 
«y ETES ceux à 25 Mc/s au Laboratoire de Spectroscopie 
P P 
ENTER Hertzienne à la Sorbonne (Paris) ce dont les auteurs 
9) Dre 6) 14,8 AIR tiennent à remercier Monsieur le Directeur 
1 . Siegbahn et M. le Professeur R. Freymann 
LE P' K. Siegbahn et M. le Prof R. Frey 
Se CEA ; 
h) CRE CT 15,0 4227°0$ Manuscrit reçu le 20 septembre 1960. 
; th CH, OH 
OCH, (2) 
H 
i) et 15,4 668 -172 
«y CH, (3) 
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THÉORIE DE LA PHOTOEXCITATION DANS LES SOLIDES 


Par A. MEESSEN, 
C. P. N. L. Université de Louvain 


Résumé. — La théorie de l’absorption de la lumière par photo-excitation est résumée, puis 
approfondie. On montre qu’il faut distinguer entre les transitions « obliques » et « verticales ». 
Elles suivent des règles de sélection différentes et mettent en jeu des processus différents. Elles 
sont aussi caractérisées par des probabilités de transitions distinctes, et il apparaît que la pro- 
babilité de transition oblique, n’est pas, comme on l’avait cru, indépendante de l’énergie. 

Le calcul de l’ordre de grandeur et de la distribution spectrale de l’absorption par transitions 
obliques et verticales conduit aux conclusions suivantes : 1. L’absorption due aux transitions 
verticales est environ cent fois plus forte que celle due aux transitions obliques. — 2. L’absorption 
fondamentale dans les isolants et semi-conducteurs est due aux transitions verticales. — 3. Les 
transitions verticales sont interdites pour les métaux (alcalins). — 4. L’accord entre l’absorp- 
tion de photoexcitation par transitions obliques, observée et calculée, est obtenu pour un ordre 
de grandeur correct du paramètre de variation du potentiel Ua. 


Abstract. — The theory of optical absorption by photoexcitationis reviewed and further inves- 
tigated. It is shown that a distinction is to be made between « oblique transitions » and « vertical 
transitions». They correspond to different selection rules and imply different physical processes. 
They are also characterized by distinct transition probabilities, and it appears that oblique tran- 
sition probabilities, is not energy independent as it has been assumed in the past. 

The calculation of magnitude and spectral distribution of the absorption by oblique and ver- 
tical transitions permits the following conclusions : 4) Vertical transitions give an absorption 
that is about a hundred times greater than the oblique transition absorption. 2) The fundamental 
absorption in insulators and semiconductors is a vertical transition absorption. 3) Vertical 
transitions are not allowed for (alkali) metals. 4) Agreement between observed and calculated 
photoexcitation absorption by oblique transitions is achieved with the right order of magnitude 


for the potential variation parameter Ua. 


Introduction. — La théorie de l'effet photo- 
électrique nécessite : 1) la description du processus 
d’excitation d’un électron du solide consécutive à 
labsorption d’un photon, c’est ce que nous appe- 
lons le processus de photoexcitation et 2) l’ana- 
lyse du processus d’émission, c’est-à-dire de la 
sortie du solide de l’électron ainsi excité. 

La théorie de la photoexcitation fait partie de 
la théorie de l’absorption optique des solides. On 
sait que les ondes électromagnétiques sont amorties 
dans les métaux, par suite de la conductibilité « qui 
exprime dans quelle mesure les électrons du 
métal sont mis en mouvement par le champ élec- 
trique de l’onde. Avec l’indice de réfraction com- 
plexe n' — n(l — ix), qui inclut le coefficient 
d'extinction x, on peut écrire pour l’amplitude de 
l'onde 


mn \. ”- 
AR Va DRE 


Le coefficient d'absorption de l'intensité lu- 
mineuse sera donc 


« = 2n xwfce ou «x = krofc 
puisqu'on montre facilement que nxw — 2ro. On 


peut done aussi exprimer le pouvoir d'absorption, 
à partir du coefficient de conductibilité o exprimé 


en s-! ou en @Q°.s/g.cm?. L’absorption de con- 
ductibilité a pû être décrite dans le cadre des 
théories classiques [1], mais la comparaison de la 
théorie avec l’expérience n’est satisfaisante que 
pour l’infrarouge et non pas pour la lumière vi- 
sible et surtout ultraviolette. Dans cette région 
spectrale on doit tenir compte de transitions quan- 
tiques bande à bande pour les électrons du solide. 
Des études sur ces transitions de photoexcitation 
ont été effectuées par de Kronig [2], Bethe [3], 
Wilson [4], Butcher [5] et Dexter [6]. L’applica- 
tion à la théorie de l'effet photoélectrique de 
volume a été faite par Mayer et Thomas [7]. 


La probabilité de transition. — Pour calculer : 
le coefficient d'absorption «(v) ou le coefficient de 
«conductibilité » o(v), il suffit de calculer le nombre 
s(v) d'électrons excités par la lumière de fréquence v 
par s et par cm du métal, puisque le théorème 
de Poynting prévoit pour l’énergie électromagné- 
tique absorbée par s et par cmÿ 


G(V) EH 2E=ES (NE AV: 
Soit s(k; v) la probabilité de transition par s et 


par em sous l’action de la lumière de fréquence y, 
pour un électron dont le nombre d’onde initial 


est k. Comme le nombre d'électrons par em et par 
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élément de volume de l’espace des X est égal à 
2/8 x? d°k, on aura 


.S(v) = fs, v) dt. 


Les électrons du solide sont décrits par une 
fonction d’onde de Bloch, normalisée au volume 
de périodicité V : 

Var) = VU un) .eifer 


avec 
Limétdr=n et ff, lumsltdr= 1.1 


Plus loin nous sous-entendons l’indice n qui 
caractérise la bande d’énergie, en considérant une 
transition e(k) — e’(x’). Les fonctions u(r) pos- 
sèdent la périodicité du réseau, que nous suppo- 
sons à symétrie cubique. 


ux(T) = ux(r + ma), M = (Ma, Mo, Ma). 


Cette périodicité permet d'écrire le dévelop- 
pement en série de Fourier 


2 5 = = ITS 
uE(T) — D ag eiGer 


où 2r/a est le vecteur de base du « réseau réci- 
proque », tandis que a est le vecteur de base du 
« réseau ». On aura aussi un développement sem- 
blable pour le potentiel 


U(r) = Ÿ Ugeiêr. 
G 


Remarquons que l'énergie e(k) et que la fonc- 
tion u; sont périodiques dans le réseau réciproque : 
e(R+ G) = e(X) et uf+a/7) — uw). 
La probabilité de transition est donnée par 
Q2 


SLR, s) = 2€ KA > 1e PE] 


avec Q — v’(k) — e(k) — hv. Nous obtenons donc 
la conservation de l’énergie €’ — v + hv puisque 
pour un temps { assez grand, il reste 

AR, ) = |< FIRE ZE 8(Q). 


L’hamiltonien d'interaction entre l’électron et 
le photon est 


si le photon est décrit par le potentiel vecteur 
À = à Ajettot+t.7) 


a est le vecteur de polarisation et n — noce est 
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le vecteur d’onde. L'élément de matrice contiendra … 
deux termes, en effet 
< k|H|k> = me fe 40 einer pE(a.V) VE dr 


mc 


eh See Ce 
EE EN ie ei + nr vX [av 
mceV ‘Jr x 


+ 18. k] ur dr 
On peut négliger n par rapport à k et 4’. De plus 
on peut se ramener à une intégration sur une seule 
cellule élémentaire de volume a. Il suffit de poser 
en effet r —r, + 2m. a, la somme étant étendue 
aux À cellules élémentaires du volume V,r, étant 
pris dans une cellule élémentaire à partir d’un des . 
sommets de la cellule cubique. On obtient à cause 
de la périodicité de u(r) : 


<KIHR> = AE euF—Fira 
mceV 
ie ei%—F).r0 us [a.V - ia. k] U% d? To- 
a 


L'élément de matrice ne sera appréciable que 
si les exponentielles de la somme sont égales à : 
Punité, c. à d. pour #4 — k — G —2rn/a. Ceci 
exprime la conservation de l’impulsion. En re- 
marquant encore que le potentiel scalaire est nul 
pour la radiation, et que À = cE/]w, on a 
ER : e—48.r US [a.v + ia.k] ux dr. 


mo a Ja 


L'ETAT 


Les transitions verticales et les transitions 
obliques. — La transition de photoexcitation est 
donc régie par les relations de conservation de 
l’énergie et de l’impulsion 


= e+hv et K=k+G—=R+ 7. 


Nous distinguons (fig. a) les transitions ver- 


ticales, avec k°— k et les transitions obliquss, 
avec k' = k + G(G Z 0). On a toujours identifié 
ces deux types de transitions, en se basant sur la 
périodicité en G de l’énergie e, et en ne considérant 
que la première zone de Brillouin. Cette distinction 
est cependant essentielle. 

1) Ces deux types de transition représentent 
des phénomènes physiquement distincts. Les tran- 
sitions obliques correspondent à une réflection de 
Bragg de l’électron sur le réseau cristallin en même 
temps qu’il absorbe l’énergie Av du photon inci- 
dent. Dans les transitions verticales par contre il 
n'y a pas de réflection de Bragg. L’interaction 
entre l’électron et le photon, demande évidemment 
toujours la présence d’un troisième corps pour 
réaliser la conservation de l’énergie et de l’impul- 
sion, Tamm et Schubin [8] en déduisaient pour 
l'effet photoélectrique de volume, que la photoex- 


Se 


N°3 


citation doit nécessairement être accompagnée 
d’une réflection de Bragg sur le réseau cristallin. 
Cela ne serait vrai que pour un électron parfaite- 
ment libre, décrit par une parabole d’énergie. 
En fait les transitions verticales sont permises, 


Etk)ev 
Na(1,1,0) 


Tr verticales 


F1G. 1. — Schéma pour la photoexcitation 
par transitions verticales et par transitions obliques. 


à cause de la périodicité des diagrammes e(k), ce 
qui est d’ailleurs aussi une conséquence du fait 
que l’électron est soumis au potentiel périodique 
du réseau. Notons encore qu’un électron parfai- 
tement libre ne pourrait pas non plus être excité 
dans une transition oblique, avec réflection de 
Bragg, puisque cet électron ne ressent pas la 
périodicité du potentiel. Nous trouverons d’ail- 
leurs plus loin que la probabilité de transition est 
nulle dans ce cas. 

2) En utilisant l’approximation des électrons 
quasi-libres, on peut mettre en évidence la diffé- 
rence entre ces deux processus de photoexcitation 
d’une façon très simple. 
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Nous supposons donc que 


ñ? k2 
2m 


et e(k') — 


e(k) — 


pour les valeurs de Æ et £’ qui interviennent dans 
les transitions étudiées. Cette approximation est 
assez justifiée pour les métaux alcalins. La fi- 
gure la donne en effet le diagramme e(k) pour la 
direction (1, 1, 0) du Na, d’après les calculs de 
Slater [9], Wigner et Seitz [10]. Nous y avons 
également indiqué le niveau de Fermi p et le 
niveau du vide W. L’échelle des énergies y est 
donnée en eV. 

Dans le cas des transitions obliques, on obtient 
pour la condition de conservation de l’énergie et 
de l’impulsion, lors de l’absorption de la lumière 
de fréquence v 


2 LE = 
Q(X) = LR + CNE EN 
avec 


A2 G2 
TITI 


—1) avec Va — 


La surface Q(k4) — 0 est donc le plan perpen- 
diculaire au vecteur G, à la distance k, de l’ori- 
gine, dans l’espace des 4. Ce plan est représenté 
à la figure 1. Puisque les électrons sont dans leur 
état initial à l’intérieur de la sphère de Fermi 
de rayon kr, au moins à00K,on a —%r < ko< + kr 
donc 


hv, = hvall — 2kr/G) < hv < hvs = hall + 2kr/G). 


A chaque valeur de G, correspond donc une région 
d'absorption comprise entre hv, et Av. Comme 
on a pour les métaux alcalins Ær — (6rx2)1/3/a 
on obtient pour le seuil d'absorption 


Rv, — 150 sl — 1,24 /n) 


si Av, est en eV et a en ampères. Mott et Jones [11] 
ont montré que les seules valeurs acceptables de 
n = (n1, ñ2, na) sont telles que n, + nr: + mm 
soit un nombre pair. On peut done avoir n — /2 
pour (1,1,0); n = 2 pour (2, 0, 0); nr —1/6 pour 
(2, 1, 1) ; etc. Le tableau I donne les valeurs de 
hv et de Av, pour les métaux alcalins et différentes 
valeurs de n. 

Il est important de remarquer que les électrons 
excités par la lumière de fréquence v dans des 
transitions obliques, possèdent un spectre d'énergie 


€ = € + hv < €’ <u + hv 
avec 
ñn2 kè 
2m 


ha 


(4 — AvJhva)? = 


€o — 7 
on peut exciter en effet des électrons, dont le 
nombre d’onde k, < k < kr. À la figure {a on 
trouve les deux transitions obliques extrêmes pour 
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une fréquence v. €, est l’énergie minimum des 
électrons qui peuvent être excités. 

Dans le cas des transitions verticales, la condition 
de conservation d’énergie et d’impulsion est vé- 
rifiée pour 

A2 k2 ñ?2 k2 


HE AE Re 


—"hv —10 


Q(k) = & — € 


si nous avons deux paraboles opposées comme dans 
le cas de la figure 1-a, caractérisées par les masses 
effectives m, et m,. La surface Q(k) — 0 est donc 
maintenant la sphère de rayon k, comprise dans la 
sphère de Fermi. À cause de la périodicité de 
e(k), on voit que le seuil d’absorption est encore 
kv, comme pour les transitions obliques. Mais on 
remarque que les électrons excités pour une fré- 
quence v ne correspondent maintenant qu’à une 
seule énergie €’ — € + hv. 

3) La probabilité de transition s(k, v) est diffé- 
rente suivant qu’on. a des transitions verticales 
ou des transitions obliques. Pour les transitions 
obliques on a l’élément de matrice 


M = Mo = ferai [av + ia fludr 


et pour les transitions verticales on a l’élément de 
matrice 


M = M, = fuÿa.v ux dr + ia. [ui ux dr, 


Pour estimer ces éléments de matrice il faut 
connaître les fonctions d’ondes uz et uy. Dans le 
cas des transitions verticales, on voit que l’élé- 
ment de matrice en a.V est nul si les deux états 
sont de même parité, par exemple deux états s. 
Nous aurons donc 


M M4=ia.kfuurdrk pour Ar +1. 


Si par contre les deux états sont de parité op- 
posée, un état s et un état p par exemple, le premier 
terme sera prédominant. 


RE M = fu a.v Ur dr Cite pour Te": 


Cette remarque a été faite par Dexter [6] quine 
considère que l’élément de matrice des transitions 
verticales pour l’absorption optique dans les semi- 
conducteurs. Dans le calcul de l’absorption op- 
tique des métaux on n’a pas fait [1] [5] de distinc- 
tion entre M, et M,. On a aussi laissé tomber 
l’élément de matrice M,, en admettant [3] l’or- 
thogonalité des fonctions ug(r). Ainsi le calcul de 
Butcher [5] ne fait intervenir que la première 
partie de l’élément de matrice M4. Nous remar- 
quons cependant que la seconde partie de W, 
ne peut pas être négligée, même si on admet 
l’orthogonalité des fonctions ux(r), à cause du 
facteur exponentiel. 
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Comme la fonction d’onde 4,x{r) vérifie l’équa- 
tion . 


ñ? A 
(55 + 0) dur — eo) ÿnz 
on aura pour les fonctions u,z(r) 
RA  in2— r? k? 
[is Le EE + OU] umé = enlh) un 
m m 2m 


Puisque cet opérateur est hermétique, on a 
l’orthogonalité des fonctions u,zx(r) et uyx(r) 


PRE 


qui sont en fait les fonctions qui interviennent 


dans M+ où 4 — k: Mais pour avoir cette ortho- 
gonalité, il faut encore que le potentiel U soit le 
même pour l’électron dans l’état initial et dans 


l’état final. Si nous supposons donc que l’électron ” 


soit faiblement lié dans l’état initial et que l’élec- 
tron soit libre dans l’état final, nous aurons 
M} 7 0. Si nous considérons par contre un élec- 
tron libre ou un électron faiblement lié dans un 
même potentiel aussi bien pour son état initial 
que final, nous aurons M}; — 0. 

Comme on peut toujours poser le développement 
en série 


uk (r) —Dageitr et uy (F) = Zageiêr 
on peut écrire pour les éléments de matrice 
ia. k[Z a_g.aa] a 


M- = ia. G(D a_g.aalat avec 


Ma 
CEE 00! et 


Ma = ia.(G +%) Zae.2e +0] avec GÆ0. 
G? : 


Ces éléments de matrice sont donc bien diffé- 
rents. Pour les estimer nous supposons que l’élec- 
tron soit lié dans son état initial au potentiel 
périodique 


U=Y Ugeid:r 


et que l’électron soit libre dans son état final. 
Pour. l’électron libre on aura ux{(r) = 1, il en 
résulte que M_ — 0. Si le potentiel U ne présente 
que de faibles variations, on peut appliquer un 
calcul de perturbation pour calculer la fonction 
d’onde dans l’état initial. On trouve 

+. UG pe, 


dt 2m 0 


Avec UE 0 


avec k'=k+G. 


Donc 


On obtient donc finalement pour les éléments de 


matrice 


Mx+ = ia.kas; M = 0: Ma = — ia. F TE, 
V 


M4 est donc proportionel à k' = k + G et non pas 
à G comme il résultait du calcul de Butcher. La 


8:23 


probabilité de transition n’est donc pas indépen- 
dante de l’énergie de l’électron excité [5] [71]. 
Aussi notre probabilité de transition sera égale 
au double environ de la probabilité de transition 
de Butcher, près du seuil Av. 


Le coefficient d’absorption pour les transitions 
obliques. — Nous avons vu que le coefficient 
d'absorption « est proportionnel au coefficient © 
qui sera donné par l’expression suivante 

ANS 


Et gs 0, v) d°x 


6(v) 


n° e 2 
_ 2h (5 JAUT 8(Q) d° k. 


Considérons maintenant un élément de surface 
dS sur la surface Q(k) — 0. Comme on aura dans 
la direction normale à cette surface dQ = v,Q.dx, 
on peut écrire d 4 — dS dQ7/Vz Q. Il en résulte 
que 

AC 
_ 2hvm?r2 4 


1M? 


ro ds. 


o(v) 


Pour les transitions obliques, on aura un effet 
de polarisation dans l’absorption optique si on a 
un monocristal avec une direction bien définie 
de G. En effet M, —a.k". Pour la lumière pola- 
risée Îl G, on aura donc M, = G-k, > G-kr > kr 
et pour la lumière polarisée L G, on aura 
MA; — ki < kx. Pour la lumière non palarisée ou 


2 10 sec”! 


Cr621 


(o) 
1 2 J 
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quand la direction de G est quelconque, on peut 
prendre la valeur moyenne 


up = Mai = À (T5) 8 


hy 


Pour les transitions obliques, la surface Q(k) — 0 
est le plan k,. 


: 2 
Q(H = (GE + 2ke Ch et VO = EC. 


Posons 4’? — k1? + (G-k,)? et dS — 2rk1 dky 
avec 0 < ky SV — #2 : et'soit set sx —e—e, 
les énergies correspondantes à k’et à Æ,. On aura, 
si on intègre sur €’ : 


U% e? m fu + Av)? — (se, + 2] 
= nl LUE : 


On peut mettre © sous une autre forme, si on 
intègre sur 1 : 


s() Uëe?m ee AIT = 


| &4n G 78 UNE 
[* — hv;) (hve — hv) + 2(3hva — hv) | 
8(Ava)? 


Au dernier facteur près, nous obtenons ainsi 
la formule calculée par Wilson [4] et Butcher [5], 
qu’on obtient si on utilise G au lieu de k’ dans M4. 

Les points expérimentaux, auxquels nous allons 
comparer nos résultats théoriques, sont repris de 


A 


hr(eV) 


Fi. 2. — Points expérimentaux et courbes théoriques pour l’absorption de photoexcitation, exprimés par le 


coefficient de conductibilité, pour quatre métaux alcalins. 
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Butcher [5], qui a déduit l’absorption de conduc- 
tibilité classique [1] de l’absorption observée ex- 
périmentalement. La figure 2 donne ces valeurs 
«expérimentales » de o(v) pour les métaux alcalins. 
Les courbes en trait plein sont les courbes théo- 
riques o(v) calculées d’après les formules précé- 
dentes, en adaptant la valeur de U, et en utili- 


sant la plus petite valeur possible » — \/2. 

L’allure générale des courbes est semblable à 
celle des courbes de Butcher, mais nous trouvons 
un accord avec les points expérimentaux pour des 
valeurs de U% de l’ordre de 2 à 3 eV, tandis que 
Butcher doit introduire des valeurs de l’ordre 
de 0,3 à 0,5 eV. 

Mais il est possible de calculer la valeur de U4 
pour les métaux alcalins par une méthode indé- 
pendante. C’est la méthode de Bethe [3] que nous 
avons appliquée au cas des métaux alcalins. 


Ua = 2,28 ZU3 F'(b)[aÿ avec b — Z—-V3n/[2a 


où U, est en eV et a en ampères. La fonction 
F'(b) est le facteur de forme connu de la théorie 
de la diffusion des rayons X. Le tableau I contient 
les valeurs calculées ainsi pour U, avec les trois 
premières valeurs de n. On voit que les valeurs 
calculées pour n — 1/2 sont de l’ordre de grandeur 
de 4 eV, ce qui est en bon accord avec nos valeurs 
« expérimentales ». 

Les valeurs du seuil d'absorption Av, indiquées 
dans le tableau I, montrent qu’on doit tenir 


TABLEAU I 

n Na K Rb Cs 

a(A) 4,24 5,25 5,62 6,05 
hva(eV) v2 16,7 10,9 9,5 8,2 
2 33,4 21,8 19,0 16,4 
V6 50 32,6 28,5 24,6 

hv,(eV) V2 2,07 1,35 1,18 1,02 
2 12,7 8,6 SE) 6,4 
v6 24,7 16,1 18,1 12,4 
UaleV) V2 4,5 k,1 5,3 6,0 
2 3,1 3,0 4,0 4,8 
V6 2,4 2,5 3,4 4,1 
UaleV) exp. 2,0 22 2,6 3,2 
hv,(eV) v2 31,0 20,4 17:8 15,3 
u(eV) 3,1 2,1 1,7 1,5 
W(eV) 5,3 4,3 3,7 3,4 


compte, plus loin dans l’ultraviolet, des autres 
valeurs de n. Ceci explique sans doute la diver- 
gence entre les points expérimentaux et la courbe 
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théorique pour le Cs puisqu’un second spectre 


d'absorption commence à 6,4 eV. Mais notons de 
plus qu’on n’est pas fixé dans ce cas, sur la vali- 
dité de l’approximation de l’électron libre, qui a 
servi dans tous ces calculs. ; 

Nous pouvons conclure que l’absorption optique 
des métaux alcalins, est très bien décrite en sup- 
posant une photoexcitation par des transitions 
obliques. 


Le coefficient d’absorption pour les transitions 
verticales. — Le coefficient o(v) est encore donné 
par la même formule de départ. Mais maintenant 


2 
ME = M4 = 0. 


Puisque la distribution de k est uniforme pour 
les électrons qui peuvent être excités, la polarisa- 
tion de la lumière sera sans influence. 

La surface Q(k) — 0 est maintenant la sphère 
de rayon # et dS — k? do. Si nous avons deux 
paraboles d’énergie opposées par leur sommet, 
comme dans la figure 1a, caractérisées par des 
masses effectives m, et m, avec y = m;fm2 on 
aura 


= 2 
' 2 Er 
OUT CPE (4 + y) k + AE — hy 
Q ñ2 . 
et VRQ = HER 
Finalement 

ñ e2 29/2 mu? e2 (AE — hy)#? 
c(v) NE BR PS ——— 

hv rm;(1 + y) (1 + y)? 77? hv 


La variation de o(v) est donc bien différente 
pour les transitions verticales et pour les tran- 
sitions obliques. Si nous calculons o(v) pour le cas 


du Na:avec-AF= 95 °m, =met ms 120 
on obtient 
UE 3/2 
M _ en s—1, 


On obtient ainsi des valeurs de o(v) qui sont de 
l’ordre de 100 fois plus grandes que les valeurs 
obtenues pour le Na dans le cas des transitions 
obliques. Nous pouvons conclure que les transi- 
tions verticales doivent être interdites pour l’ab- 
sorption optique ds métaux alcalins. Ce résultat 
est assez évident puisque nous pouvons admettre 
l’orthogonalité des fonctions wy, parce qu’en fait 
l’électron ressent le même potentiel dans son état 
initial et final. 

Mais cela n’est pas toujours vrai. Dans les semi- 
conducteurs et isolants [6], on observe en effet 
une absorption fondamentale, dûe aux transitions 
des électrons de la bande de valence à la bande 
de conduction, qui est décrite par une formule 
du type 


(ms [m)/2 C—5/2 1,54 105 (hv — AE)Y? 


œ(v) = Ts en cm—1, 


No 3 


où la valeur de C dépend de la configuration 
des diagrammes d’énergie. Quand on a deux para- 


—., 


Ê —..— 


Cs, Sb 


Morgulis 
[o) Burton 
A Wallis 

me ne Spicer 


—..— TJaft Philipp 


x(v)=1,1-1 


0° (hv=-16)7? 
by 
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boles qui sont opposées par leur sommet, avec leurs 
concavités proches l’une de l’autre, on a C = 1 + +. 
Cette formule correspond à l’élément de matrice 
M — M}, —k. On obtient une dépendance en 
(Av — AË)V?/Jhy quand on a des transitions avec 
Ar = --1et. M= .M2= Cie, 

La figure 3 donne le coefficient d’absorption «(v) 
pour labsorption optique du semiconducteur 
Cs3Sb. On y trouve les résultats expérimentaux 
de Morgulis [12], Burton [13], Wallis [14], Spi- 
cer [15], Taft et Philipp [16] ainsi que la courbe 
théorique 


a(v) = 1,1 406 (Av — 1,6)%2Jhv 


ce qui correspond très bien à l’ordre de grandeur 
prévu théoriquement. Les déviations des points 
expérimentaux par rapport à cette courbe, sont 
dûes à des impuretés du réseau. Notons d’ailleurs 
que l’absorption des isolants et semiconducteurs 
peut être compliquée par des transitions « indi- 
rectes », qui correspondent à une photoexcitation 
avec absorption ou émission simultanée de pho- 
nons de vibration du réseau [6]. 

De toute façon nous pouvons conclure que les 
transitions verticales ne sont pas interdites dans 
les isolants et semiconducteurs. L’orthogonalité des 


fonctions ux n’y est donc pas réalisée, ce qui de-. 


mande que l’électron ressent un potentiel différent 
suivant qu’il est dans la bande de valence ou dans 
la bande de conduction. Les transitions obliques, 
avec réflection de Bragg ne sont pas interdites 
non plus, mais leur intensité est négligeable par 
rapport aux transitions verticales. 

Nous remercions M. le Professeur de Hemptinne, . 


1 2 3 4 5 hytev) pour ses suggestions et intéressantes discussions. 
F1c. 3. — Déterminations expérimentales du coefficient ; | 
d'absorption pour le semiconducteur Cs,Sb. Manuscrit reçu le 23 septembre 1960. 
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QUELQUES PROPRIÉTÉS DE SYMÉTRIE INTRINSÈQUES VIS-A-VIS DU GROUPE D'INVERSION 


EN THÉORIES DES PARTICULES ET LE RAPPORT SPIN-STATISTIQUE 


Par K. H. TZOU, 


Institut Henri-Poincaré, Paris. 


Résumé. — Dans les théories générales des particules ayant pour équations de champ celles 
du type (1)-(5), on démontre d’abord des propriétés de symétrie intrinsèques de certains opérateurs 
matriciels en rapport avec chacune des opérations d’inversion, A l’aide de ces propriétés intrin- 
sèques, il est montré que les inversions antiunitaires-unitaires ou celles unitaires-antiunitaires 
peuvent toutes déterminer la statistique des particules élémentaires suivant leur spin, si l’on 
postule l’invariance des règles de quantification des champs par rapport à l’une quelconque de 
ces opérations. Les inversions unitaires-unitaires ou celles antiunitaires-antiunitaires, au contraire, 
sont incapables de la détermination des statistiques. On remarque en particulier que l’inva- 
riance CPTM, comme l’invariance CPT, n’est pas une invariance automatique en théorie quan- 
tique des champs ; elles n’y sont assurées que si les particules élémentaires obéissent à leur bonne 
statistique respective suivant leur spin. 


Abstract. — In the general theories of particles with field equations of the type (1)-(5), we 
study at first some intrinsic symmetry properties of certain matrix operators in connection with 
each of the inversion operations. With the aid of these intrinsic properties, it is shown that the 
antiunitary-unitary inversions or the unitary-antiunitary ones can all determine-the statistics of 
elementary particles according to their spin, if we postulate the invariance of the field quanti- 
zation rules under any one of these operations. The unitary-unitary inversions or the antiunitary- 
antiunitary ones, on the contrary, are incapable for the determination of particle statistics. ‘It is 
particularly remarked that the CPTM invarianee, like the CPT invariance, is not an automatical 
invariance in quantum theory of fields ; they are assured in this case only when elementary par- 
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ticles obey their respective good statistics according to their spin. 


1. Introduction. — Dans des articles précé- 
dents, nous avons étudié l’invariance du champ 
de spin 1/2 par rapport à un groupe de symétrie 
composé d’inversions de l’espace (P), du temps (T'), 
de la masse () et de la charge (C), groupe d’in- 
version à seize éléments (opérations d’inversion) [1]. 
Nous y avons précisé en particulier l’identification 
des états propres d’une particule de spin 1/2 
engendrés par chacune des seize opérations de ce 
groupe. 

Cette étude et ses conclusions se généralisent 
facilement aux autres théories des particules. La 
généralisation en est particulièrement évidente et 
directe quand les équations des champs sont 
exprimées toutes dans une expression semblable à 
l’équation de Dirac, à savoir, 


Pu(du — ieay) (x, m) + mY(x, m) = 0. (1) 


Les matrices f, obéissent à des règles non com- 
mutatives qui sont en général différentes suivant 
le spin des particules. Par exemple, 


Bu Bv + Pv Pu — 20uv (2) 
en cas de spin 1/2, et 
Ba Bu By + By Bu Ba — au Bv + uv Ba (3) 


en cas de spin 0 ou 1. Dans une théorie où les 
spins 0 et 1 sont superposés d’une certaine façon 


(champ vectoriel général sans condition supplé- 
mentaire) [2], nous avons aussi, au lieu des rela- 
tions (3), les règles moins restrictives que ces 
dernières, 


Ep(Bx Bu By — dau By) = 0 (4) 


où p désigne toutes permutations des indices Auv. 
Pour un certain champ général composé de spins 
8, 8-1, s-2, ..., 1/2 (ou 0), les règles les plus géné- 
rales que peuvent vérifier de toute façon les ma- 
trices P, sont 


HE ENTER 
Iin=0 s s S 


al NI NI NI 
£ Zp dits Otratta -.. Ô 


Uon—iton Puonsi +: Puosr1 — 


où s° — $ + 1/2 pour s — demi-entier-impair, et 
! 


$ —$ pour s — entier. p désigne maintenant 
toutes les permutations des indices 


M1 Ma Pa “+ Mos+i: 

Dans ce qui suit, pour raison de simplicité sans 
pourtant perdre la généralité, nous supposerons 
que les B, sont des matrices hermitiennes. 

Aux paragraphes suivants, nous allons traiter 
de la question de l’invariance des règles de quanti- 
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fication des champs par rapport à chacune des 
opérations d’inversion. À ce sujet, nous ne nous 
bornerons plus au seul. champ de spin 1/2 et nous 
étendrons notre étude à tous les spins selon les 
théoriés qui ont pour équations de champ celles du 
type (1)-(5). Nous ne nous intéresserons pas à la 
structure détaillée de ces théories des particules. 
Mais les relations de commutation des matrices 
Bu, (2)-(5), seront très utiles dans nos considéra- 
tions ultérieures. ÿ 

Dans la théorie quantique des champs, l’inver- 
sion du temps T est une opération toute parti- 
culière par le fait que T engendre une transforma- 
tion antiunitaire des vecteurs d'état Q(t) du sys- 
tème des champs [3], tandis que la transformation 
de Q(t) induite par P, M ou C est toujours une 
transformation unitaire. En effet, dans la repré- 
sentation de Schrüdinger par exemple, les vecteurs 
d'état satisfont à l’équation de Schrodinger, 

ù 


iS Qalt) — H Qu(é) (6) 


où À est l’opérateur hamiltonien du système des 
champs considérés, fonctionnelle hermitienne des 


opérateurs d'émission et d'absorption. Sous l’opé- 
ration 7°, nous avons évidemment 


Où + J Q, J' a Je (7) 
H'= 3H 3, 


où a’ désigne l’état propre des champs correspon- 
dant à l’état a dans la référence (r, — t), et H 
l’opérateur hamiltonien correspondant. La rela- 
tion (7) représente une transformation antiunitaire 
des vecteurs d’état sous l’inversion 7’, transforma- 
tion qui est naturellement la même aussi dans 
toutes autres représentations (celles d’Heisenberg, 
d'interaction, etc.) que dans la représentation de 
Schrodinger. 

Soient © et À deux grandeurs quelconques des 
champs. L'élément de matrice de Q par exemple 
entre états a et best, selon (7), 


< alQlb > = < b’|Q'la >, (8) 

où 
"= 0 J* (9) 
est l’opérateur correspondant à Q dans la réfé- 


rence (r, —t). D’après (8) et (9). il est évident 
que, sous l’opération 7, 


< a|QR|b > =+-< W|(QR)a > = < b'|RQ'a >. (10) 

Quand il s’agit d’une transformation unitaire des 
vecteurs d’état Q sous P par exemple (ou sous © 
ou M), l’invariance des éléments de matrice 
s’exprime, au lieu de (8), par 


< alQlb > = < «91 >, (11) 


QUELQUES PROPRIÉTÉS DE SYMÉTRIE INTRINSÈQUES 


148 
où 
ET 01, 
(07 = 44 OX sf = T4, à 


Dans ce cas, au lieu de (10), nous avons 
<a a|QR|b LE ee < a”|(QR)"1b" > = LE a”|Q"R"|b" > 
(12) 


Les relations (10) et (12) seront très utiles dans 
les considérations ultérieures des commutateurs 
où anticommutateurs des opérateurs des champs. 

Selon leur façon d’opération à la fois sur les 
opérateurs des champs et sur les vecteurs d’état 
comme nous avons précisé plus haut, les seize 
éléments du groupe d’inversion se divisent en 
quatre ensembles différents : 


AZI, P, M, PM; 

B=C; CP, CM/CPM;: 

V= AT FÉPT MT PMTS 
Z=Bliæ= CL CP CMISCPMT. 


En effet, les opérations À sont des inversions 
unitaires-unitaires dans le sens qu’elles engendrent 
chacune une transformation unitaire des-opéra- 
teurs des champs et une transformation unitaire 
aussi des vecteurs d’état. Au même sens, les opé- 
rations B sont des inversions antiunitaires-uni- 
taires, les YŸ antiunitaires-antiunitaires, et enfin 
les Z unitaires-antiunitaires. On remarque donc 
que P, T, C'appartiennent à des catégories d’opé- 
ration différentes. 

Dans ce qui suit, nous allons montrer que 
l’invariance des règles de quantification des 
champs (1) par rapport aux inversions À et Y 
n’apporte aucune information sur là statistique 
des particules, mais elle exige pourtant l’unitarité 
de la matrice représentative de chacune de ces 
opérations sur les opérateuis des champs. Quant 
aux inversions B et Z, au contraire, l’invariance 
des règles de quantification par rapport à l’une 
quelconque d’entre elles fixe automatiquement la 
forme propre de ces règles, ce qui établit le rapport 
entre le spin et la statistique des particules. De 
plus, tout comme dans le cas des À et Y, cette 
invariance exige-t-elle aussi l’unitarité de la ma- 
trice représentative de chacune des inversions B 
et Z et encore aussi celle de la-matrice 0 définis- 
sant les champs adjoints. 

Le rapport spin-statistique des particules élé- 
mentaires est une question beaucoup étudiée anté- 
rieurement, notamment par Pauli [4] [5]. Pauli a 
examiné ce problème en détail en partant des 
postulats divers [4], par exemple ceux que l’éner- 
gie totale d’un champ quantifié et la mesure de la 
probabilité physique doivent être non négatives. 
Mais, il s’agit par là des propositions qui, à priori, 
ne font pas partie de la structure interne des 


Â4k 


théories des particules elles-mêmes. Cependant, 
Pauli et Belinfaute [5] ont examiné aussi la possi- 
bilité de déterminer les statistiques des particules 
par postuler l’invariance C. Leur conclusion est 
affirmative. Par l’invariance C ou par celle sous 
toute autre opération d’inversion, il s’agit bien 
des conditions sur la structure même des théories 
des particules. 

De son côté, Schwinger a étudié la possibilité de 
déterminer les statistiques des particules en postu- 
lant linvariance des sue de quantification par 
l’inversion du temps 7 [3]. Mais Watanabe a 
ensuite montré que ce n’est pas l’invariance T 
mais l’invariance CT qui permet de déterminer les 
statistiques [6]. Sa conclusion, pourtant, n’est pas 
en contradiction avec celle de Schwinger. En fait, 
c’est bien effectivement l’invariance CT mais non 
lPinvariance T qu’a considérée Schwinger, si T est 
défini comme l’opération antiunitaire-antiunitaire 
de Wigner [7] au lieu de l’opération unitaire-anti- 
unitaire définie par Racah et Pauli [81, celle-ci 
étant en réalité l’inversion CT. 

Dans ce qui suit, nous allons considérer l’inva- 
rilance par rapport à toutes les seize inversions en 
vue d'examiner la possibilité de la détermination 
des statistiques des particules. D’une telle étude, 
nous croyons que, d’une façon exhaustive, est 
résolu le problème sur la détermination des statis- 
tiques des particules élémentaires par l’invariance 
vis-à-vis des inversions. En ce qui concerne les 
invariances C, T, CT, nos conclusions seront en 
accord avec celles de Pauli et Belinfante. de 
Schwinger et de Watanabe. 

On remarque enfin que, dans les études anté- 
rieures sur le rapport spin-statistique, un postulat 
fondamental a toujours été essentiel, postulat que 
nous pouvons appeler « condition de causalité ». 
Selon cette condition, deux observables quelcon- 
ques commutent toujours, s’ils sont localisés à deux 
points de l’espace-temps séparés d’une distance du 
genre espace. En d’autres termes, deux observables 
de ce genre sont mesurables simultanément, car la 
mesure de l’un ne doit causer aucune perturbation 
à celle de l’autre, étant donné qu’il n’existe pas de 
propagation physique ayant une vitesse supérieure 
à la vitesse de la lumière. C’est là une condition 
physique fondamentale de la mesurabilité, qui 
restera essentielle aussi dans la présente étude sur 
les règles de quantification des champs. 


2. Propriétés de symétrie intrinsèques de cer- 
taines matrices en rapport avec les inversions dans 
les théories des particules. — Dans le cadre général 
des théories des particules ayant pour équations 
de champ celles du type (1)-(5), nous allons 
d’abord démontrer certaines propriétés intrin- 
sèques relatives à leur structure en rapport avec 
les inversions. À ce sujet, nous devons bien distin- 
guer les inversions unitaires des inversions anti- 
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unitaires vis-à-vis des opérateurs des champs et 
les traiter séparément. 

D'abord, les champs adjoints Ÿ doivent vérifier 
les équations du type 


Bu Qu + ieay) 


où l’indice supérieur T désigne la transposition des 
matrices. 

En comparant (13) à l’ 
les champs adjoints par 


PTx, m) — mET (man) = 10 08) 


équation (1), on définit 


Fix, m) = Y'(x, m) 0, (14) 


0 étant une matrice non singulière vérifiant les 
conditions 


OP — — 60, 064 — B0, (1—1,2,38). (15) 

_ On peut alors poser en général 
OÙ = a0—1,  (a* = a > 0), (16) 
6700), (b, arbitraire). (17) 


Cela fait, nous pouvons maintenant examiner 
la question de l’invariance sous les diverses opé- 
rations d’inversion. 

19 INVERSIONS UNITAIRES. — Soit 

a= I, P, M, PM, CT, CPT, CMT, CPMT. 


L’invariance de l’équation (1) sous l’inversion a 
s'exprime par (voir [1]) 


Var, m) = Te F(a(x, m)), (18) 
To Bu La = dB, (19) 


où 0% — + 1 ou — 1 suivant a et u, et l’opéra- 


tion a(x, m) change (x, m) selon la définition de 
l’inversion a. De (19), on peut poser en général 
(Ts, matrices non singulières) 


Pot fa lo VUS ee A0 (20) 
Ta = cale (caca = 1) (1). (21) 

Puis, de (15) à (19), nous tirons 
01 = 00, (= 1). (22) 


On peut démontrer que la valeur de «, ne dépend 
que de l’opération a et du spin des particules consi- 
déré, mais non de la représentation des matrices 
Bu Prenons une autre représentation par Bu que 
nous supposons matrices hermitiennes aussi. On a 


évidemment 
Bu = upyurt, (ut = ut) 


(*) Répétée une seconde fois, l'opération a sur (18) nous 
donne 


Voa(x, m) = Ta Va(a(x, m)) = TS Ÿ(x, m), 


d’où, en général, T2? — eiv, 


| 
| 


N°3 
Dans cette nouvelle représentation, en général, 


= yu0ut, mi = Gaulau-t, (4, Ca, constantes). 


Alors 


’ RE 
RS OA Lo 2 


c’est-à-dire que 
(23) 


La conclusion reste valable, même si les 8, 
ou Bou toutes les 8, et B, ne sont pas hermi- 
tiennes. 

Si nous nous référons par exemple à la théorie 
générale des particules à spin par méthode de 
fusion de L. de Broglie [9], la valeur de «, est facile 
à déterminer pour chacune des a (2). Nous la 
donnerons pour chaque inversion a dans un tableau 
à la fin de ce paragraphe. 

Finalement, selon (18), (17), (20) et (22), nous 
avons 


Ua — Age 


Pr, mn Fax, mn) — PR (a(x, m)) To (24) 
Ka —= Xg a (25) 
29 INVERSIONS ANTIUNITAIRES. — Soit 


b= C, CP, CM, CPM,T, PT, MT, PMT. 


L’invariance des équations (1) et (13) par rap- 


port à l’inversion b s'exprime par (voir par 
exemple [1]) 

; Pô(æ, m) = Ty VT(b(æ, m)), (26) 

To Bu To = 09 Bu, (27) 


où également 3 — + 1 ou — 1 suivant bet u, et 
l’opération b(x, m) change les coordonnées (x, m) 
selon la définition de l’inversion b. D’après (27), on 
peut poser en général (l'}, matrices non singu- 
Jières) 

(28) 


D =h16,. (= fo > 0), 


FAN (29) 


D'autre part, de (15), (17) et (29), on tire facile- 
ment 


ep T6, (e5 = 1). 


Ty, 07 15 — 0, (30) 


(2) Dans la théorie de L. de Broglie, un champ à 2s 
indicess pinoriels, Féos...ær eSt en général une certaine 


superposition des spins s, s —1,..., 1/2 (ou 0). Le champ 
adjoint est défini avec 
| A ME Ye 
et la conjugaison de charge C, par exemple, définie par 
CICM:CEr "10, 
CT RO SN UE 1, 2, à, 2e). 


[ci 7e sont les matrices de Dirac agissant sur l’indice 
spinoriel «x, et C® est la matrice représentative corres- 
pondante de la conjugaison de charge. Dans ce cas, il est 


_ évident que, par exemple, 


dd ten 1e 
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Comme dans le cas des inversions unitaires, les 
valeurs de &, et «, sont ici aussi indépendantes de 
la représentation des matrices B,. En effet, dans 
la nouvelle représentation avec les 6,,, on a 


D yuOu—-1, Ty =yulyur, (x, @&, constantes). 


Il est alors évident que 


LA 
XD — Xp, 


4 —\ €}. (31) 

Cette propriété intrinsèque de «, et e, est d’ail- 
leurs indépendante de l’hermiticité ou non-hermi- 
ticité des matrices B, ou ,. Leurs valeurs ne 
dépendent que de l’inversion b et du spin des parti- 
cules considéré. Si l’on se réfère par exemple à la 
théorie des particules par méthode de fusion [9], il 
est facile de déterminer &; et «, pour chacune des 
inversions b(?). Nous en donnerons les valeurs dans 
un tableau plus loin. 

Finalement, selon (26), (16), (17), (28) et (30), 


nous avons 


Vox, m) = Vb(x, m) = x» VT(b(x, m)L5*, (32) 
Xp = bp Afb. (33) 

Puis (?) 
F{b(x, m)) = e5xo" To Po (x, m). (34) 


TABLEAU 


VALEURS DE &g,@b ET €p SUIVANT L’INVERSION ET LE SPIN S 


INVERSION a INVERSION &p - € 
I 1 C (— 1)?s (—1)?s 

12 1 CP (— 1)?s 1 
M (— 1)?s CM 1 (— 1)?s 
PM : (— 1)?s CPM 1 (— 1)?s 
CT (— 1)?s 4} 1 (— 1)?s 
CPT (— 1)?s PTE 1 (—1)?s 

CMT 1 MT (— 1)?s il 
CPMT 1 PMT (— 1)? (—1)% 


(8) Répétée une seconde fois, l'opération b sur (26) nous 
donne 


Pov(x, m) = Ty FO (b(x, m)) = xp cp Vx, m), 


à ue _ 2 
d’où, en général, ep €p à = 1,oulaf|—= 1" Donc a? ff = 1, 
condition qui ne détermine pas les valeurs de a et de fs. 


En général, on ne peut pas poser x — e de sorte que 
L' 


bb — Ÿ, même quand b est une opération indépendante 
de P et de 7’. En cas de C par exemple, xg = eç nous 
conduit à poser afo — 1, condition bien réalisable. Mais 
en cas de CM, on aurait afom = (— 1)°5, ce qui est irréa- 
lisable pour spins demi-entiers-impairs. 

La relation inverse de (34) s’écrit 


or, m) = (eo x5 1) (xs 5) Ty PT (b(x, m)), 


relation qui est compatible à (26), car ep, &p, a, fo sont 
toutes des constantes réelles. 
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3. Invariance des règles de quantifiquation par 
rapport aux inversions. — Pour les règles de 
quantification des champs, nous pouvons les 
écrire d’une façon générale suivante : 


Yalæ, m) Pe(a’, m) — oFe(z’, m) Walz, m) 


= 1Aup(d, m) Afx — x’), (35) 
PVelz, m) Fe(z’, m) — ©’ Ye(x’, m) Fe(x, m) = 0, (36) 
Pal, m) Pe(a’, m) — ©*—1Fe(x’, m) Polx, m) — 0. (37) 


Dans ces relations, « et w’ sont deux constantes 
non nulles dont les valeurs déterminent la statisti- 
que des particules. On peut démontrer que w'? — 
(done &'#*71 = &'* — w'), et de la compatibilité 
des trois relations, on tire encore &’ = « ({). 
D’autre part, par covariance relativiste et condi- 
tion de causalité (voir $ 1), À doit être bien une 
fonction de la distance æ—x', fonction delta 
covariante dont les propriétés sont bien connues. 
(Voir par exemple [5]). Nous rappelons pourtant 
ici ses quelques propriétés de symétrie dans l’es- 
pace-temps, qui nous seront utiles, 


A(— v,i) = Afr,t), Ar, — 1) — 
A(— x) = — A(x) 


— Ar, t) (38) 


A(d, m) est un opérateur différentiel qui, de 
moindre ordre x possible en 8, à,, doit vérifier la 
condition de compatibilité avec les équations des 
champs, à savoir 


(Bu du + Mjay /\yel(d, m) = Ap(d, m) (E 


À(d, m) étant un autre opérateur différentiel qui ne 
comporte plus de facteur C1 — m?. Ainsi par 
exemple, A(d, m) = — 6,2, + m pour spin 1/2, et 
— (8,9,)? — mB,2, pour spin 0 ou 1 ou pour une 
superposition des deux spins selon les règles (4). 
D’après les relations non commutatives (2)-(5) 
des matrices B,, on peut établir une règle générale 
sur l’ordre minimum n de l’opérateur /\(d, m), 


TK m°), 


N—VpaIr, pour. S\—.entier, (39) 
n — impair, pour s — demi-entier-impair. | 
(*) De (36) par exemple, nous tirons évidemment 
Tele”, m) Felt, m) — wo! (x, m) plz’, m) = 0, 
d’où w—!— «’ par compatibilité, puis w’* = w’. D’autre 


part, mettons Ÿ’ — Y{x’, m) par exemple. Par des calculs 
simples, on montre alors que 


(CAR AA A Le MER TR OT DA 
ET RD Li DE Al tee 
== DOUTE OP Fe) Fe + DUT TR — o'VeTe) 
= i(1 — ww’) Agy(D”, m) A(x” — x) T4 
+ iw— &) Aey(d’, mm) À (x —x) Yp 
ET CN MIE Tr 
La compatibilité de (36) à (35) exige 
Ho 0, LOÙ @2,— "0 — 1 aussi 


donc wow’ = 1, 
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Après ces précisions sur les règles de quanti- 
fication (35)-(36)-(37), nous pouvons maintenant 
examiner l’invariance de ces règles par rapport aux 
inversions. Pour ce faire, nous nous rappelons qu’au 
fond, en théorie quantique, on a affaire avec les 
éléments de matrice. Dans ce qui suit, quand nous 
considérons les règles de quantification, nous 
pensons donc aux éléments de matrice des opéra- 
teurs en question, bien qu’on perte toujours les 
signes < >. 

Commençons d’abord par les inversions À, qui 
sont des opérations unitaires-unitaires. D’après 
(12), (18), (24) et (35), nous avons, pour les nou- 
veaux opérateurs des champs Ÿ4 et Ÿ, les règles 
de quantification suivantes : 


Paz, m) VÉ(x, m) — oVÉ(z’, m) Val, m) 
= %4 l'aay la68 (Fy(A(æ, m)) Fs(A(x’, m)) 
— wVs(A(z’, m)) Fy(A(x, m))] 
= ixa(Ta À (A@,m)) 12" )as À (A(x—2) 
(40) 


= KA NA A\ag(d, m) A(z FA æ); 


car, selon (19), (39) et (38). on démontre facilement 
que 
14 A (, m), 


1 pour 1, P, 
(— 1)?s pour M, PM, 


T4 À (AO, m)) Ta — 
(41) 


NA = AA = 
A(A(x — z')) = Alx — a’). 


D’après (40), la seule condition de l’invariance 
des règles de quantification par les inversions A est 


FAUNAËE— CIE 
d’où 


lat de (42) 


Cela signifie l’unitarité des matrices [°,, condi- 
tion nécessaire pour l’invariance des règles de 
quantification par rapport aux inversions A. Mais 
cette invariance À n’apporte aucune information 
sur la constante w, qui reste indéterminée. 

Considérons ensuite les inversions antiunitaires- 
unitaires B. Selon (12), (26), (32) et (35), les règles 
de quantification des nouveaux opérateurs des 


\1rB , . : 
champs W° et WF, s’écrivent maintenant 
PEL m) PE a, m) — FE (x, m) VB(x, m) 


= — cul roy lrô6 (Vs(B(z’, m)) Fy(B(x, m)) 


— w—1Ÿ,(B(x, m)) Vs(B(x’, m))] 
= — 1oxg(l'B AT(B(— d, m)) T3 Je 
A(B(x" — x)), SO @: 


LOXZ Ur: A\ag(d, m) A(x ai + z'), (43) 


I 


PEN EC EN 


NAS 


car il est facile de démontrer, avec (27), (39) et 
(38), que 


LB AT(B(— D, m)) L&° = n8 À (0, m), 


es à se il HOUTIC AC: Lu 
PRO Er CM CP 
A(B(z&— x)) = A(æ — x) = — A(x — x’). 


Ainsi, selon (43), de l’invariance des règles de 
quantification par rapport aux inversions B, on 
tire les conditions 


OXB NB = 1. (45) 


Mais ©? = 1% — 1, donc x£ — 1 aussi. Comme a 
et j# ne peuvent être négatifs, il faut évidemment 
que 


afg = 1. (46) 


Dans ce cas, 


© — xp n8 = ap nB — (—1)?s. (47) 

On établit ainsi le rapport entre spin et statis- 
tique des particules élémentaires en postulant 
linvariance vis-à-vis de l’une quelconque des 
inversions B. L’invariance B est d’ailleurs sure- 
ment réalisable, si les matrices 0 et l', sont toutes 
unitaires (a = fr — 1). 

Sous les inversions antiunitaires-antiunitaires Ÿ, 
nous avons, selon (10), (26), (32) et (35), pour les 
nouveaux opérateurs des champs, WY et Ÿ”, les 
règles de quantification suivantes : 


Vale, m) FE (x’, m) — © FE (x, m) Va(x, m) 
— 47 Dyay l'rsl Val Y(z’, m)) FyY(x, m)) 
— wV,(Y(x, m)) Fa(Y(x’, m))] 
= iy(Tr ATY(— D, m)) Tr')ag A(Y(x’ — à) 


= y ny /\ag(d, m) A(x — x), (48) 
car, avec (27), (39) et (38), on démontre que 
Dr AT 2, m)) 7° = ny À (, m), 
1 : (49) 
HAT DOUTER ETS, 
COUT; A)fs pour MT, PMT, 
A(Y(z — x)) = — A(r'— x) = A(z — »). 


Selon (48), l’invariance des règles de quantifica- 
tion par les inversions Y est bien établie, si seule- 
ment 


Ainsi 


XY NY = 1. 


ie ile (50) 
condition qui est réalisable, si les matrices 0 et l'y 
sont unitaires (a = fy = 1). 

L’invariance Ÿ ne nous procure aucune informa- 
tion sur la valeur de la constante «. Elle ne peut 
donc pas déterminer la statistique des particules. 
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Étudions enfin l’invariance des règles de quanti- 
fication par rapport aux inversions unitaires- 
antiunitaires Z. D’après (10), (18), (24) et (35), les 
règles de quantification des nouveaux opérateurs 


2 ES pes 
des champs, Ÿ? et W”, s’écrivent évidemment 


Pl'ax, m) V£ x’, m) — oT (x, m) YZ(æ, Im) 
= — o%7 Lzoy Lzo8|V(Z(æ, m)) Fa(Z(a’, m)) 
— 1 Vs/Z(x’, m)) F,(Z(x, m))] 
= — ioxg(Tz À (Z(®, m)) 17” )a8 A(Z(x — x), 
sl Gt 6), 


= 10XZ NZ /\ag(d, m) A(x — x’), (51) 


car, avec (19), (39) et (38), il est facile de démontrer 
que 

Lz A(Z(OD, m)) F1 = nz À (D, m), 

1 Dour CT, CPT, (52) 
(— 1)%s pour CMT, CPMT, 

A(Z (x — æ’)) — — Afx — x). 


ns = (—1)#soz = 


Selon (51), l’invariance des règles de quantifica- 
tion par rapport aux inversions Z n’est réalisable 
que Si 


GT = ©XZ nz = À, (53) 
d’où xZ = 1, et puis 
az 1 (54). 
Alors finalement 
o = Xz nz = ag nz = (—1)?s. (55) 


L’invariance Z détermine donc la statistique des 
particules élémentaires. 

On note d’autre part que les règles (36) et (37) 
sont automatiquement invariantes par rapport à 
toutes les inversions À, B, Y et Z, car w’? — 1. Par 
condition de compatibilité, nous devons pourtant 
toujours avoir ©’ — w, (voir ({)), ce qui complète 
la détermination de la statistique des particules 
suivant leur spin par l’invariance B ou Z. 

Dans ce que nous avons conclu plus haut, nous 
remarquons en particulier le fait que l’invariance M 
est complètement incapable pour la détermination 
de la statistique des particules. Il en est de même 
de l’invariance T, comme l’a déjà remarqué 
Watanabe [6], et aussi de l’invariance PT étudiée 
par Pauli [10]. Mais la conjugaison de charge C 
peut bien déterminer le rapport entre spin et sta- 
tistique, comme l’ont conclu dans une étude Pauli 
et Belinfante [5]. Enfin, il est intéressant de remar- 
quer que l’invariance CPT ou même l’invariance 
CPMT est une condition suffisante dans la déter- 
mination de la statistique des particules élémen- 
taires suivant leur spin. Cela montre en même 
temps que CPT et CPMT ne sont pas des inva- 
riances automatiques en théorie des champs quan- 
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tifiés. Les invariances CPT et CPMT n’y sont 
assurées que dans le cas où les particules de spin 
entier obéissent à la statistique de Bose et les 
particules de spin demi-entier-impair à la statis- 
tique de Fermi, et où les matrices 0, lopr ?t 
lompr sont des matrices unitaires. (A propos de 
l’invariance CPT en rapport avec les statistiques, 
voir aussi Pauli [10]). 


4. Conclusion. — Dans le cadre général des 
théories des particules ayant pour équations de 
champs celles du type (1)-(5), nous avons d’abord 
précisé des propriétés de symétrie intrinsèques de 
certains opérateurs matriciels en rapport avec 
chacune des inversions. Ainsi, la matrice 0 définis- 
sant les champs adjoints possède, vis-à-vis de 
chacune des inversions unitaires, une certaine 
symétrie qui ne dépend que de l’opération d’inver- 
sion et du spin considérés, mais non de la représen- 
tation des matrices f,.. Il en est de même, quand il 
s’agit des inversions antiunitaires, de 0 et encore 
des matrices représentatives de ces opérations 
d’inversion aussi. 

A l’aide des propriétés de symétrie intrinsèques 
de ces opérateurs matriciels, nous avons étudié la 
possibilité de l’invariance des règles de quantifica- 
tion des champs par rapport à chacune des seize 
opérations du groupe d’inversion. Il à été montré 
que l’invariance par les inversions unitaires-uni- 
taires (7, P, M, PM) ou par celles antiunitaires- 
antiunitaires (T, PT, MT, PMT) n'apporte 
aucune information sur la statistique des particules 
élémentaires, mais elle impose cependant l’unitarité 
à la matrice représentative de chacune de ces 
opérations. Au contraire, l’invariance par rapport 
aux inversions antiunitaires-unitaires (C, CP, CM, 
CPM) ou à celles unitaires-antiunitaires (CT, CPT, 
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CMT, CPMT) introduit non seulement la 
condition de l’unitarité aux matrices représenta- 
tives de ces inversions et aussi à la matrice 0, mais 
elle détermine encore d’une façon décisive la sta- 
tistique à laquelle doivent obéir les particules 
élémentaires suivant leur spin. Cette invariance ne 
permet que deux statistiques, celles bien connues 
de Bose et de Fermi. On établit ainsi le rapport 
spin-statistique en postulant simplement l’inva- 
riance parrapport à l’une quelconque des inversions 
antiunitaires-unitaires ou unitaires-antiunitaires. 
Au point de vue mathématique, c’est donc l’anti- 
unitarité de ces inversions qui joue le rôle décisif 
dans la détermination de la statistique des parti- 
cules. Maisles inversions deux fois antiunitaires en 
deviennentincapablesencore. Au fond du problème, 
c’est l’invariance C qui est effectif et décisif dans la 
détermination des statistiques des particules. 

Par la présente étude, nous croyons avoir, d’une 
façon exhaustive, résolu le problème de la 
détermination des statistiques des particules élé- 
mentaires par l’invariance des règles de quanti- 
fication des champs vis-à-vis des inversions. Nous 
remarquons en particulier que, malgré le fait que 
l'opmr = 1, l'opération CPMT, tout comme CPT, 
n’est pas une invariance automatique en théorie 
quantique des champs. Au contraire, l’invariance 
CPMT, comme l’invariance CPT, peut déterminer 
le rapport spin-statistique des particules. En 
d’autres termes, l’invariance CPMT, comme 
linvariance CPT, n’est assurée en théorie quan- 
tique des champs que dans le cas où les particules 
de spin entier obéissent à la statistique de Bose et 
celles de spin demi-entier-impair à la statistique 
de Fermi. 


Manuscrit reçu le 11 octobre 1960. 
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ÉTUDE A HAUTE RÉSOLUTION D’UNE RAIE DE VIBRATION-ROTATION 
DANS L'HYDROGÈNE SOLIDE 


Par H. P. GUSH (*) 
Laboratoire Aimé-Cotton, CG. N. R. $., Bellevue. 


Résumé. — Les forces intermoléculaires donnent lieu à un spectre infrarouge d’absorption dans 
Yhydrogène solide. La raie de vibration-rotation $,(0) a été étudiée à haute résolution dans un 
cristal de parahydrogène pur et dans un cristal contenant des petites concentrations d’ortho- 
hydrogène. Un couplage quadripolaire entre une molécule para excitée et une voisine ortho est 


responsable d’un dédoublement de la raie. 


Abstract. — Intermolecular forces give rise to an infrared absorption spectrum in solid hy dro- 


gen. 


The rotation-vibration line S,(0) has been studied at high resolution in crystals of pure para- 
hydrogen and in crystals containing small impurities of orthohydrogen. 


Quadrupole coupling 


between the excited para molecule and a neighbouring ortho molecule is responsible for a splitting 


of the line. 


Introduction. — L’étude du spectre de l’hydro- 
gène solide dans l’infrarouge a apporté des con- 
clusions intéressantes en ce qui concerne les inté- 
ractions moléculaires dans les cristaux. Les travaux 
expérimentaux de Gush, Hare, Allin et Welsh [1] 
et les travaux théoriques de Van Kranendonk [2] 
ont montré que le spectre de vibration-rotation s2 
compose de deux parties. L’une a pour origine une 
transition dans une molécule simple appartenant 
au réseau. L’autre se produit lorsqu'une transition 
vibrationnelle dans une molécule est accompagnée 
soit par une transition rotationnelle dans une molé- 
cule voisine, soit par l’excitation d’une vibration du 
réseau ; ces transitions s’appellent des transitions 
doubles. 

La présente étude porte uniquement sur la raie 
fine S,(0) qui résulte de la transition 


0er 1; 0 7-2 


dans une molécule simple de parahydrogène appar- 
tenant à un réseau de molécules para. On avait 
déjà montré [1] que cette raie, de nombre d’ondes 
4 486 cm1, est très fine et qu’en outre elle est 
accompagnée de satellites faibles pour de petites 
concentrations d’orthohydrogène dans le cristal. On 
pouvait espérer éclaircir le problème de l’origine de 
ces satellites par des mesures plus précises de leur 
largeur et de leur fréquence. 


Détails expérimentaux. — Les transitions 
simples dans un cristal d'hydrogène donnent lieu à 
des raies fines, contrairement aux transitions 
doubles. Leur étude doit donc se faire à une réso- 
lution élevée ; pour étudier la raie 5,(0) dans le cas 


(*) Boursier du National Research Council of Canada 
pour les années 1957-1959. Adresse actuelle : Department 
of Physics, University of Toronto, Toronto 5, Canada. 


du parahydrogène pur il est nécessaire de travailler 
à une résolution de l’ordre de 100 000. 

Pour atteindre cette résolution élevée on a 
employé un spectromètre Fabry-Perot intégral déjà 
décrit par Chabbal et Noorman [3] et légèrement 
modifié pour ce problème. Il comprend un prémo- 
nochromateur à réseau et un seul étalon de 70 mm 
de diamètre dont les lames sont revêtues de coatings 
diélectriques. Le détecteur est une cellule à sulfure 
de plomb Kodak refroidie, de dimensions 
0,15 X 1,5 mm, employée avec un condenseur cône 
channel [4]. Le balayage se fait par variation de la 
pression d’air. La réalisation de mesures d’absorp- 
tion avec un spectromètre Fabry-Perot pose des 
problèmes de filtrage sérieux ; en effet les ordres 
élevés du prémonochromateur peuvent tomb?r 
dans une région ou la finesse du Fabry-Perot est 
très réduite. On a supprimé la lumière parasite en 
plaçant un second réseau devant le prémono- 
chromateur, en utilisant pour faire le modulateur 
(toujours employé dans les mesures infrarouge) du 
plexiglas, et en filtrant la lumière qui tombe sur la 
cellule à travers un filtre Jena BG 18. 

Lorsque le balayage se fait par variation de pres- 
sion, une variation de une atmosphère ne permet 
d'explorer qu'environ 1 em? à une fréquence de 
5 000 em. Le domaine spectral à étudier ayant 
une largeur de 5 em—t on a fait un balyage en dent 
de scie [5]. L’épaisseur de l’étalon étant connue, 
l'échelle de fréquences est calculée à partir de la 
pression dans l’étalon, cette pression étant con- 
trôlée avec un manomètre de précision. 

Les cristaux d'hydrogène solide sont préparés et 
étudiés dans un eryostat optique refroidi à l’hélium 
liquide (fig. 1). Le cristal est préparé de la façon 
suivante. On remplit à moitié le cryostat d’hélium 
liquide, la cellule d'absorption étant maintenue au- 
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dessus du niveau du liquide ; ensuite on condense 
de: l'hydrogène gazeux dans la cuve Jusqu'à ce 


qu’elle soit pleine. Le solide est obtenu en abaissant 
la cellule vers la surface de l’hélium liquide ; fina- 


Helium 


Cellule 


F1G. 1. — Cryostat d'optique. 


lement le cristal est mis en place devant les fenêtres 
du cryostat et la température de l’hélium est 
abaissée au-dessous du point À pour éliminer les 
bulles d’hélium gazeux dans le chemin optique. 

La cellule est représentée sur la figure 2. Elle 
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comprend deux cellules distinctes. La cellule supé- 
rieure À de 0,3 mm d’épaisseur est utilisée pour 
faire des cristaux de parahydrogène pur. La deu- 
xième cellule est utilisée pour les mélanges de para- 


icm. 


Fic. 2. — Cellules d’absorption. 


orthydrogène. Elle comprend elle-même deux par- 
ties B, et B, d’épaisseurs 0,3 mm et 10 mm. On 
mesure la raie S.,(0) observée dans un cristal de 
parahydrogène pur dans la cellule A ; la raie S,(0) 
observée dans un cristal de para contenant une 


- 


4,0 de 


100 para 


JAY 


| 


98% 


Log,, L/I 


03 


00 


-03 4486 +03 


FiG. 3. — Raie #,(0) pour plusieurs concentrations d’ortho-hydrogène ; température 


0,3 mm. 


\ 
mie 


-03 4486 +03 


-03 4486 +03 Fréquence (C0 


: 2 °0K ; épaisseur du cristal : 


Nas 


impureté d’ortho dans la cellule B; ; ; les satellites 
faibles observés dans un cristal épais, cellule B,, 
ayant la même composition ortho-para. Cette mé- 
thode à permis une mesure précise de l’écart en 
fréquence des satellites de la raie centrale, sans 
avoir établie la fréquence absolue transmise par 
l’étalon. Une connaissance très précise de la fré- 
quence de la raie #,(0) n’était pas nécessaire pour 
la présente étude, et on’a adopté la valeur donnée 
par [1] de 4 486, 0 em, 


Résultats et discussion. — I. LA RAïIE, S,(0). — 
Le profil de la raie telle qu’elle est observée pour 
une épaisseur de 0,3 mm est indiqué dans la 


0,4 
AT 
02 
L | 
(Q) 2 4 6 8 
% d'ortho 
Fic. 4. — Largeur à mi-hauteur, Ao, de la raie #,(0). 


figure 3 pour plusieurs concentrations d’ortho- 
hydrogène. La raie est très fine dans le cas du 
parahydrogène pur et sa largeur à mi-hauteur ne 
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vaut que 0,12cemt - (La largeur de la fonction d’ap- 
pareil est inférieure à 0 05 em! et son effet sur le 
spectre peut être négligé). La largeur croît rapide- 
ment en fonction de la concentration d’ortho et déjà 
pour une concentration de 8 % elle est quatre fois 
plus élevée que pour le cas du para pur (fig. 4). 
On a pu évaluer le coefficient d’absorption inté- 
gré défini par la somme : 
5 fl 1 T(v) 
= à f<0e 1(5) dy 

ou on pose : N — nombre de molécules par centi- 
mètre cube ; { — l’épaisseur de l’échantillon ; v — 
fréquence de la transition ; Z,(v) = l’intensité de 
lumière incidente sur l’échantillon à la fréquence v ; 
I(v) — l'intensité de lumière transmise par l’échan- 
tillon. Il est évident que si la forme de Z(v) est 
changée par l’emploi d’un spectromètre à faible réso- 
lution cette somme ne garde pas sa propre valeur ; 
pour avoir la valeur exacte de & il faudrait corriger 
le spectre observé pour l’effet de la fonction d’appa- 
reil. Dans la présente étude la largeur à mi-hauteur 
de la fonction d’appareil est inférieure à la moitié 
de la largeur de #,(0) ; les corrections sont par 
conséquent très petites et elles ont été négligées. 
La moyenne de trois déterminations donne 
à = 3,27X 10715 cm st pourdla raies, (0)fdans 
le parahydrogène pur à 2 ©K. Ceci est à comparer 
avec la transition Q,(0) + S,(0) pour laquelle & 
Vaut.œe 4,90 xX1081%260m°2 HO) ESrrapporr 
est 1/13 ce qui est en accord raisonnable avec les 
calculs théoriques de Van SÉEREUER qui à prévu 
un rapport 1/15 [2]. 


(1) Cette valeur était déterminée à l’Université de : 
Toronto avec un spectromètre à réseau. 
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Fic. 5. — Raïies satellites observées avec un cristal d'épaisseur 10 mm 
; épaisseur 0,3 mm. Température : 


S,(0) dans le para pur 


| 14860 


14880 Fréquence (cm ) 


: concentration d’ortho 1 %. En pointillé la raie 
2 0K. 
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IT. LES RAIES SATELLITES. — La figure 5 montre 
le spectre observé dans le mélange : 1 ©} ortho, 
99 % para, l’épaisseur du cristal étant 10 mm. On a 
dessiné sur la même graphique la raie S,(0) comme 
elle est observée à travers une épaisseur de 0,3 mm 
de parahydrogène pur. Le spectre comprend une 
forte raie d'absorption centrale, S,(0), plus des 
raies satellites dont l'intensité croît en fonction de 
la concentration d’orthohydrogène ; elles ne sont 
pas présentes dans le parahydrogène pur. Les 
écarts de fréquence entre la raie centrale et les 
satellites sont indiqués dans le tableau I. 


TABLEAU I 


ÉCcARTS DES SATELLITES DE LA RAIE #,(0) 
ÉRas 


SATELLITE (*)} 


ÉcarT 


— 2,492 + 0,01 cm1 
+ 0,868 + 0,01 cmt 
+ 1,804 + 0,01 cm—! 
+ 1,945 + 0,01 cm—t 
+ 2,290 +.0,01 cm1 


® À © 5 » 


(*) Suivant la figure 3. 


Étant donné l'importance des transitions 
doubles dans le spectre de l’hydrogène solide, à 
première vue la présence des raies satellites dans 
les mélanges ortho-para s'explique par des tran- 
sitions pour lesquelles la transition 5,(0) dans une 
molécule de parahydrogène est accompagnée par 
un changement d'orientation d’une molécule voi- 
sine d’orthohydrogène. On a déjà fait appel à e2 
procédé pour expliquer la structure de la branche Q 
pour des faibles concentrations d’ortho [1]. Bien 
que ce point de vue soit partiellement valable dans 
le présent cas, il est certainement simplifié et il ne 
permet pas un calcul de la position exacte des 
satellites. Néanmoins ceux-ci sont une conséquence 
du fait que les molécules ortho et para excitées 
peuvent prendre plusieurs orientations dans le 
réseau. 

L’état final d’une molécule de parahydrogène qui 
subit la transition $S,(0) est v — 1, J — 2. Cet 
état possède cinq niveaux de dégénérescence pro- 
venant du fait que la composante du moment 
angulaire dans une direction fixe, m, peut prendre 
cinq valeurs. Dans un cristal de parahydrogène pur, 
les forces intermoléculaires anisotropiques sont r2s- 
ponsables d’un dédoublement en trois niveaux selon 
le cas: m—0, +1, +2. [2] Seul le niveau 
m = + 2 est actif en absorption dans le para- 
hydrogène pur, et la raie S,(0) est simple. 

Lorsque des molécules d’orthohydrogène sont 
ajoutées au réseau il peut y avoir une perturbation 
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supplémentaire du niveau v = 1 ; J — 2 qui pro- 
vient de l'interaction entre les moments quadri- 
polaires de la molécule para excitée et l’impureté 


d’ortho. Autrement dit, le système de deux molé- » 


cules voisines para-ortho possède une dégénéres- 
cence de quinze pour la transition S,(0) parce que 
la composante du moment angulaire de la molécule 
ortho peut prendre elle même chacune de trois 
valeurs : 0, + 1. Cette dégénérescence est soulevée 
par des forces intermoléculaires et une structure 
compliquée en résulte. à 

III. Le spectre théorique. — Le potentiel d’inter- 
action entre deux molécules d'hydrogène peut être 
mis dans la forme d’une somme : 


ee © 
Pirelli ad 
n 


où À est la distance entre les deux molécules, 
—> — 3 


et &,, w, sont les orientations des molécules rela- 
tives à un axe intermoléculaire. Il est montré par 
Van Kranendonk [2] que les seuls termes aniso- 
tropiques qui sont effectifs dans un cristal d’hydro- 
gène sont V,, et Vos. Le premier représente le 


2 
couplage quadripolaire. Nous avons e(R) — De 
où Q est le moment quadripolaire et az2 —= 1, 
dy = — 4, ay —= 6. Le deuxième terme repré- 


sente un champ effectif cristallin et il peut être mis 
sous la forme : 


_— 
Ve = (716) (4x)/2e YA oi) 


— 
ou & est une constante pour le cristal, et o; est 
l’orientation de la molécule relative à l’axe hexa- 
gonal. 

On a calculé l'effet de ces deux potentiels pertur- 
bateurs. sur lé niveau = 1: J = 2908 
J — 1 d’un système de deux molécules voisines 
para-ortho dans le cristal. La position théorique des 
raies que l’on peut en déduire est indiquée dans la 
figure 6. Il y à un dédoublement en vingt huit 
niveaux distincts dont deux sont deux fois dégé- 
nérés. Il y a deux valeurs possibles de l’angle entre 
l’axe intermoléculaire para-ortho, et l’axe hexa- 
gonal du cristal ; l’énergie du système est différente 
dans chaque cas à cause du champ cristallin. C’est 
la raison qui fait que l’on trouve trente niveaux 
d'énergie au lieu de quinze comme on a prévu par 
un raisonnement simple. Les détails du calcul 
seront présentés dans une autre communication [6]. 

L'accord entre la fréquence calculée des raies et 
la fréquence observée est très satisfaisant et on 
peut en conclure que les satellites sont effecti- 
vement dus à la perturbation d’une molécule para 
par une molécule voisine d’ortho. Certaines raies 
prévues par la théorie ne sont pas observées ; une 
comparaison complète demanderait un calcul des 
intensités et ce calcul est actuellement entrepris. 
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Les constantes e — Q?/5R5 et , sont des para- 
mètres dans les équations pour l’énergie des niveaux 
et on peut faire varier ces paramètres pour avoir le 
meilleur accord avec les observations. Néanmoins 
les valeurs de # et # ne peuvent pas être choisies 
indépendamment. Une observation de l'effet 
Raman de l’hydrogène solide par Bhatnagar [7] a 
montré que la raie de rotation pure, S,(0), se 
compose de trois composantes. Selon les calculs de 
Van Kranendonk [2] l’écart en fréquence entre les 
composantes de ce triplet symétrique est donné 
par À = 5(0,903 & + 1/6 &). La valeur expéri- 
mentale de À est 2,00 cm1. € et € ont été choisis 
telle qu’elles satisfassent cette relation et on trouve 
Docs —=,——0,45 cmt. 

La valeur de € qu’on trouve ici est en excellent 
agrément avec la valeur qu’on peut prévoir à partir 
du moment quadripolaire théorique et une con- 
naissance de la distance intermoléculaire dans le 
cristal. La meilleure valeur du moment quadri- 
polaire peut être évaluée à partir des calculs récents 
de Kolos et Roothan [8] : on trouve Q — 0,49 ea. 
La distance intermoléculaire est connue à partir 


de la densité de l’hydrogène solide : R — 7,09 a,. 
On a donc € — 0,59 cmt, 

Le spectre théorique comprend un groupe de 
raies centré sur la fréquence 4 486 em! dans une 
bande d’environ 1 cmt de large (fig. 6). Il est 
vraisemblable que l’élargissement de S,(0) en fonc- 
tion de la concentration d’orthohydrogène comme 
il est montré dans les figures 3 et 4, soit dû à ces 
raies. 

Jusqu'ici il n’y a pas de théorie permettant de 
calculer la largeur de la raie S.,(0) dans le para- 
hydrogène pur. Il est remarquable que cette raie 
soit aussi fine. Un traitement théorique de ce pro- 
blème serait très intéressant. 

En terminant, je tiens à remercier M. P. Jac- 
quinot, Directeur du Laboratoire Aimé-Cotton, qui 
a bien voulu m’y accueillir pour effectuer ce travail : 
ainsi que d’autres recherches, M. R. Chabbal, qui 
a mis à ma disposition le spectromètre Fabry- 
Perot qu'il avait mis au point, ainsi que M. J. Van 
Kranendonk, avec qui j’ai eu des discussions sur 
les points théoriques. 


Manuscrit reçu le 15 novembre 1960. 
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CALCUL DE LA FORME DU PROFIL 
DE LA RAIF DE RÉSONANCE DU MERCURE ÉLARGIE PAR L’ARGON 
BASÉ SUR LA THÉORIE DE JABLONSKI 


Par S. LEGOWSKI, 


Laboratoire de Luminescence de l’Académie Polonaise des Sciences, 
Institut de Physique de l’Université de Torun, Pologne. 


Résumé. — Se basant sur la théorie de Jablonski, on a calculé la distribution des intensités 
de la raie de résonance du mercure élargie par la pression de l’argon, supposant que le potentiel 
d'interactions des atomes du mercure et de l’argon est du type Lennard-Jones [12-6]. Les para- 
mètres du potentiel ont été calculés sur les données thermodynamiques et spectroscopiques trouvées 
dans la littérature. Le profil déterminé de la raie a été comparé avec les résultats expérimentaux 
obtenus par Michels et ses collaborateurs. On a constaté le bon accord des calculs théoriques avec 
l'expérience. 


Abstract. — The shape of the mercury resonance line broadened by argon was calculated on 
the basis of the theory given by A. Jablonski. Assumption was made that the interaction poten- 
tial is of the Lennard-Jones [12-6] type. Potential parameters were determined from thermo- 
dynamic and spectroscopie data given in the literature. The theoretical calculations were com- 


pared with Michel’s experimental data. 
observed one. 


Introduction. — L'idée principale sur laquelle 
s’appuie la théorie de Jablonski [7], [12] de l’élar- 
gissement des raies spectrales par la pression, est la 
similitude du mécanisme pour la formation des 
spectres de bandes lors des passages électroniques 
dans les molécules et du mécanisme pour l’élargis- 
sement des raies spectrales. Ces deux phénomènes 
sont dûs aux mouvements des noyaux des atomes. 
Dans le domaine des spectres de bandes, tout au 
moins dans un des états électroniques entre lesquels 
il y a des passages, il s’agit d’un spectre « discret » 
des mouvements des noyaux (des états de vibration 
de la molécule). Tandis que pour le problème d’élar- 
gissement, les mouvements des noyaux dans les 
deux états électroniques, font partie du spectre 


continu. C’est là la seule différence de ces deux 


mécanismes. Leur analogie fait penser que les deux 
phénomènes devraient être traités suivant des mé- 
thodes théoriques les plus semblables possible. C’est 
d’ailleurs de cette façon qu’ils ont été traités par 
Lenz [18], [20] — aussi bien sa théorie d’élargis- 
sement des raies spectrales (qui est le dévelop- 
pement de la théorie de Lorentz) que la théorie de 
formation du spectre de bandes, se basent sur l’ana- 
lyse de Fouriér des vibrations de l’électron per- 
turbés par les mouvements des noyaux. La théorie 
classique de la formation des spectres de bandes a 
bientôt cédé à la théorie basée sur le principe de 
Franck-Condon. 

En 1931, Jablonski proposa l’application du prin- 
cipe de Franck-Condon, dans sa forme primitive, 
et des courbes du potentiel pour expliquer le phéno- 
mène d’élargissement des raies spectrales par la 
pression. Il attira également l’attention sur le rôle 
des forces de Van der Waals dans ce phénomène (1). 


(1) En plus des théories d’élargissement des raies spec- 
trales par collisions, il existait alors des théories d’érargis- 


The shape determined is in close agreement with the 


Kubhn a calculé la distribution des intensités dans 
la raie élargie, se basant sur la dépendance du 
potentiel de Van der Waals de la distance donnée 


par London, et a mesuré la distribution des inten- 


sités de la raie Hg 2 537 À en absorption, élargie 
par l’argon. Il a constaté le bon accord entre la 
distribution calculée et les observations, pour une 
certaine région de la raie élargie. Les expériences 
de Minkowski [29] sur les raies D du sodium ont 
donné des résultats concordant avec la distribution 
de Kuhn. Dans ses calculs, Kuhn a seulement pris 
en considération les collisions simples. Margenau 
[21], [24] a élargi cette théorie aux collisions mul- 
tiples. Les théories ci-dessus statistiques ou poten- 
tielles, s’appuient sur la forme primitive du prin- 
cipe de Franck-Condon. Toutefois, on ne peut 
attendre des résultats satisfaisants de l’application 
de ce principe qu’en le considérant sous sa forme 
mécanique-quantique. Weisskopf [33] a été le 
premier à appliquer cette forme, dans le but de 
justifier les théories classiques des collisions 
(A. Jablonski considère que ses raisons ne sont pas 
convainquantes [12] A Jablonski [7], [12] a 
donné la théorie générale basée sur le principe de 
Franck-Condon. Le spectre continu de l’énergie du 
mouvement des noyaux se trouve changé en spectre 
€discret » par l’inclusion de l’atome absorbant (ou 
émetteur) dans une sphère, avec l’atome pertur- 
bateur (?). Ensuite, la méthode est la même que 


sement par les champs électriques des ions, des dipôles et 
des quadrupôles de Debye et Holtsmark. Voir les articles 
de Weisskopf, Margenau et Watson [24] et de Breene [1] 
et les références citées par eux. 

(?) Il faudrait en principe traiter les atomes émetteurs 
et tous les perturbateurs comme une seule grande molécule. 
Oùn obtient quand même une bonne approximation si l’on 
considère d’abord un couple d’atomes — un atome émet- 
teur et un atome perturbateur — et si l’on calcule ensuite 
l'effet des atomes perturbateurs de certaine densité. 


PRE de 
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pour la théorie des spectres de bandes. La distri- 
bution des intensités dans le spectre discontinu 
obtenue de cette façon se trouve transformée en 
distribution d’intensités dans le spectre continu, en 
multipliant par la densité des niveaux dans l’état 
final. Différents moments cinétiques du couple en 
collision, doivent être pris en considération (avec 
leur poids statistique propre). L'expression géné- 
rale pour la distribution des intensités dans la raie 
d'absorption, est : 


ea 32 N lmam |. 2 
(AE) — ROME TES 2 G + (21414) Afro dE, 
(1) 
où 
ae fé Var) entr) dr, Zi —1 (2 
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avec NV, nombre des perturbateurs par 1 cm8 ; Y, 
masse réduite de l’atome émetteur et d’un pertur- 
bateur ; À, rayon de la sphère ; Æ,, et Es, énergies 
du mouvement des noyaux des états électroniques 
supérieur etinférieur ; Pia et Win, fonctions propres 
correspondantes ; !, nombre quantique du moment 
cinétique. Dans cette expression, la somme par 
rapport à « ! » a été remplacée par une intégrale. 

L'application de la fonction W. K. B. et de 
certaines autres approximations dans les calculs, 
conduit à une distribution pratiquement identique 
à celle de Kubn : 


&kr K3/n 
L. Lena 


2C2 re) 


re) SEAT (3) 
où les potentiels d'interactions dans l’état excité et 
dans l’état fondamental (branche d’attraction) 


sont : 


K = C' —C"; k, constante de Boltzmann ; 7, 
température absolue. Cette distribution n’est une 
approximation suffisante que pour une certaine 
région de la raie élargie. Elle est divergente au 
centre de la raie. 

Afin d'obtenir une distribution des intensités, 
valable pratiquement pour toute la raie élargie, 1l 
faut appliquer d’autres fonctions propres au lieu 
de la fonction W. K. B. et tenir compte, dans les 
courbes potentielles, non seulement de la branche 
d’attraction, comme on l’a fait jusqu’à présent, 
mais aussi de la branche de répulsion. On aurait 
également dû considérer la distribution des vitesses. 
Mais, pour le moment, les calculs n’ont été effectués 
que pour une vitesse moyenne à une température 
donnée. 

Le but de ce travail est le calcul, le plus précis 
possible, de la distribution des intensités de la raie 
Hg 2 537 À, élargie par l’argon pour la région la 
plus large possible et la comparaison de la distri- 

bution obtenue avec l’expérience. 
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Potentiels d’interactions. — Si l’on désire déter- 
miner la forme du profil de la raie spectrale en se 
basant sur la théorie de Jablonski, il faut avant 
tout définir les potentiels d'interactions de l’atome 
émetteur et de l’atome perturbateur. 

L’atome émetteur se trouve dans le champ d’in- 
teractions des atomes (ou des molécules) envi- 
ronnants. L'énergie potentielle mutuelle des atomes 
du mercure et de l’argon est décrite avec la meilleur 
approximation par le potentiel de Lennard-Jones 


r12-6] 
Pre [E) — ()] (4) 


avec V(r), énergie potentielle de deux atomes dans 
la distance r ; e, valeur minimum d’énergie poten- 
tielle ; 6, valeur de r pour laquelle V = 0. On peut 
définir les potentiels d’interactions entre les atomes 
du mercure et les atomes de l’argon pour l’état 
électronique fondamental et deux états excités, 
d’après les données expérimentales se trouvant 
dans la littérature [2], [6], [13], [26]-[28], [301 1311. 

Le potentiel d'interactions entre des atomes dif- 
férents (quand les atomes sont dans les états fonda- 
mentaux) peut être déterminé expérimentalement 
en tenant compte des propriétés du gaz telles que : 
la viscosité, le deuxième coefficient virial, etc... 
Mais, jusqu’à présent, il n’y à pas de données assez 
précises sur les mélanges de gaz ; la difficulté repo- 
sant en partie sur le fait que les propriétés des mé- 
langes de gaz dépendent, non seulement des forces 
agissant entre les atomes de natures différentes, 
mais aussi des interactions entre les atomes de 
même nature. On peut toutefois déterminer les 
paramètres du potentiel d’interactions des atomes 
différents avec les interactions des atomes iden- 
tiques. Les règles empiriques des combinaisons [5], 
[26] reliant les paramètres des potentiels d’inter- 
actions des atomes différents et les paramètres 
correspondants connus pour les atomes semblables, 
ont la forme : 


il 
OxB —= 2 (5% Lin og) (5) 
et 
Exp — V Ex € (6) 


AVEC : Gus, €ags Constantes de force décrivant des 
interactions de natures différentes, tandis que o,, 
&, Ce, €g sont les constantes de force décrivant les 
interactions de même nature. Profitant des rela- 
tions données ci-dessus et des données pour 
l’argon A, [6], [13], [28], [30], [31] et pour le 
mercure Hg: [2], on a déterminé les valeurs sui- 
vantes des paramètres pour l’état fondamental : 


Heat == 020,20 
CHgA — il 6 À (7) 


(rHgA)m = 3, 
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L’atome du mercure dans l’état excité 6 %P, et 
l’atome de l’argon dans l’état fondamental, for- 
ment, pendant leur collision, un couple d’atomes de 
deux nombres quantiques possibles : Q — 1 et 
Q — 0. Il existe donc deux courbes de potentiel 
pour l’état excité, avec les poids statistiques, res- 
pectivement 2 et 1. 

Jusqu’à présent, on n’a pas effectué de mesures 
thermodynamiques pour les atomes excités. C’est 
pourquoi il faut prendre en considération les 
données spectroscopiques pour déterminer les cons- 
tantes de forces. On peut déterminer ces constantes 
d’après la distribution des intensités dans l’aile du 
côté des courtes longueurs d’ondes de la raie du 
mercure 2 537 À élargie par l’argon, sur la base de 
la théorie statistique de Kuhn [15] et d’après la 
position des bandes des deux côtés de la raie ato- 
mique du mercure. Les constantes calculées pour 
les états excités Q = 1 et Q — O0 sont les sui- 
vantes : 


chHealk = 421,7 0K 
oMra « = 3,18 À 
(rHeA)m — 3,57 À 


cHea[k — 272,7 0K 
cHea — 3,25 À (8) 
(rHga)m = 3,64 À. 


Le bon choix des paramètres des potentiels d’in- 
teractions sera prouvé par la concordance entre le 
déplacement de la raie spectrale calculé et l’expé- 
rience. Ce déplacement a été calculé sur la base de 
la théorie statistique de Margenau [21] comme 
valeur moyenne statistique de l’énergie de la pertur- 
bation. Les calculs numériques du déplacement de 
la raie de résonance du mercure 2 537 À perturbée 
par largon (52 0C — 81,2 amagat) donnent la 
valeur suivante : De. — 18,31 cm1. Cette valeur 
s’accorde bien avec le résultat expérimental donné 
par Michels et ses collaborateurs [27] : 
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Caleul de la forme du profil de la raie. — Connais- 
sant les potentiels d'interactions on peut calculer la 
forme du profil de la raie de résonance du mercure 
en trouvant la solution de l’équation de Schro- 
dinger pour ces potentiels et en calculant les inté- 
grales de recouvrement (de Condon) correspon- 
dantes (2). L’équation de Schrodinger pour les 
potentiels de Lennard-Jones est une équation diffé- 
rentielle qui possède une singularité essentielle pour 
r — OÜet ne possède donc pas des solutions analy- 
tiques. En appliquant la solution de Langer [17], 
[3] on obtient une distribution des intensités 
valable aussi au-delà des intervalles dans lesquelles 
l’approximation W. K. B. donne des résultats 
suffisants. La solution approximative de Langer 
peut s'inscrire, à l’aide de la fonction d’Airy, de la 
façon suivante : 


a V2 


: 


Er = 


ORUE 1/6 
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r (5) 
(ue) + 


où 


Wi= flQrar pour 0?< 0. 


Cette solution est une fonction régulière pour le 
point de retour classique r. Ayant les fonctions. 
propres de l’équation de Schrodinger pour le pro- 
blème étudié, on peut passer au calcul des inté- 
grales de recouvrement correspondantes. Il est pra- 
tiquement impossible de calculer analytiquement 
ces intégrales pour la fonction de Langer, sans de 
grandes hypothèses simplificatrices. Il faut done 
procéder à des calculs numériques, pour obtenir les 
meilleurs résultats possibles. Pour déterminer nu- 
mériquement les valeurs des intégrales de recou- 
vrement, il faut tout d’abord établir des tableaux 
adéquats des fonctions propres, calculer les pro- 
duits correspondants de ces fonctions et ensuite, 
intégrer numériquement. On a calculé les valeurs 
des fonctions propres dans l'intervalle r — 0, 
r — 11 À. Cet intervalle correspond pratiquement 
au rayon d'interactions interatomique du type de 
Van der Waals pour le mercure et pour l’argon. On 
a calculé environ 450 points pour chaque fonction 
propre. Les valeurs de l’argument u de la fonction 
d’Airy : 


ee ON DE. Ti 2 UL +1 
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ont été calculés numériquement à l’aide de la règle 
des trapèzes. Les fonctions propres tabularisées ont 
été normalisées. On a calculé les produits des fonc- 
tions propres et on a intégré numériquement. Dans 
quelques cas, on a vérifié également l’orthogonalité 
des fonctions. Les intégrales des carrés des fonctions 
et des produits de deux fonctions différentes ont été 
déterminées comme suit : dans l’intervalle où les 
fonctions de Langer sont définies, on a intégré 
numériquement, utilisant les formules de 
Weddle [25], tandis que dans l'intervalle r — 11 À, 
R = 10% À, on a appliqué la solution asymptotique 
de l’équation de Schrôüdinger : 


(11) 


V4 = A cos (#2) (12) 
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. 
et les intégrales correspondantes ont été calculées 
analytiquement. La valeur du rayon de la sphère R 
a été choisie de façon à obtenir une densité sufli- 
sante de niveaux énergétiques. Pour le cas étudié, 
la densité des niveaux est 


dv'/dE, = 0,6.101$ erg-1, 


De cette façon, les fonctions propres déterminées 
donnent les valeurs des intégrales des produits des 
fonctions, lesquelles devraient être orthogonales, 
de l’ordre 10%. Ces valeurs sont 100 fois 
plus petites que les valeurs des intégrales de recou- 
vrement. 

Puisque la distribution des intensités calculée 
doit être comparée avec la forme du profil de la 
raie d'absorption, on a calculé les intégrales de 
recouvrement A», pour les passages du niveau 
de translation £" — 3/2 XT de l’état électronique 
fondamental aux différents niveaux des états élec- 
troniques excités Q = 0 et Q — 1. 

On a calculé les valeurs A,» pour quatre pas- 
sages, comme fonction du nombre quantique !, 
définissant le moment cinétique des atomes en colli- 
sion. Les valeurs suivantes des accroissements de 
_ fréquence correspondent à ces quatre passages : 
Av = — 5,06 em-1;—10,12 em-1;— 15,18 cmt ; 
10,12 em-1. Conformément à la formule (1), pour 
une valeur AE donnée, il faut calculer la proba- 
bilité des changements de l’énergie du mouvement 
des noyaux qui accompagnent les passages électro- 
niques à l’état excité Q — 0 et à l’état Q — 1. Cela 
veut dire que, pour un passage électronique donné, 
il faut calculer la probabilité des variations de 
l’énergie du mouvement des noyaux, aussi bien 
pour les valeurs positives que pour les valeurs 
_ négatives des AE. Mais, pratiquement, les passages 
électroniques à l’état excité Q — 0 sont accom- 
pagnés seulement de variations de l’énergie de 
translation AE > 0, puisque la probabilité de la 
variation de l’énergie AE < 0, se rapportant à ces 
passages, est tellement petite qu’elle se trouve dans 
les limites de l’erreur des calculs. Les passages élec- 
troniques à l’état Q — 1 sont accompagnés d’une 
façon analogue, uniquement des variations de 
l’énergie de translation AE < 0. 

Les valeurs obtenues des carrés des intégrales de 
recouvrement, en fonction du nombre quantique 
du mouvement des noyaux, sont représentées sur 
les figures 1, 2,3 et 4. 

Ayant les relations entre la probabilité des varia- 
tions de l’énergie de translation et le nombre quan- 
tique /, on peut calculer, d’une façon numérique 
simple, les valeurs de l’expression (1) et, par là 
même, déterminer les valeurs relatives des inten- 
sités aux points étudiés de la raie de résonance du 
mercure élargie. Dans les intervalles de fréquence 
de — 20 cm1 à — 200 em-t et de 8 cm! à 
20 em 1, on a calculé la distribution des intensités 
en appliquant la formule approximative (3). La 
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2 
0,5: À 
04 
0,3 
ot 
010 20 4 40 90 60 70 80 


Fi. 14. — Carré de l’intégrale de recouvrement (de Condon) 
en fonction du nombre quantique / pour la raie du mer- 


cure 2 536,5 À ; Av = — 5,06 em—1. 
A? 
0,3 
0,2 
0,1 
Z 
0 10 20 3000040 0ME DEN 60 70180: "30 
Fi. 2. — Même fonction que sur la figure 1, 
mais pour Av — — 10,12 em—t. 
A? 
0,2 
0,1 
Z 
0 10 20 30 40 50 60 70 80 90 
FiG. 3. — Même fonction que sur la figure 1, 
mais pour AV = — 15,18 cm1. 
> 
9,1 À 
14 


0 20 40 60 80 100 120 140 


F1. 4. — Même fonction sur la figure 1, 
mais pour AY = 10,12 cmt, 


distribution relative des intensités pour la raie de 
résonance du mercure 2537 À, perturbées par 
l’argon, a été déterminée de la manière décrite e1- 
dessus. On a calculé les intensités relatives par 
rapport au centre de la raie déplacée ne tenant 
compte que des collisions de l’atome émetteur avec 
un seul atome perturbateur (le déplacement de la 
raie a été déterminé sur la base de la théorie statis- 
tique). La distribution calculée est présentée sur la 
figure 5 et est comparée avec la distribution expéri- 
mentale donnée par Michels et ses collabora- 
teurs [27]. On a également porté sur la figure 5, la 
courbe de la dispersion avec la même hauteur du 
maximum, le même déplacement et la même demi- 
largeur que pour la distribution expérimentale. 
Sur la figure 5, on voit directement que les points 
calculés numériquement se superposent très bien 
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sur la courbe expérimentale. Par contre, la for- 
mule (2) donne la représentation correcte de la dis- 
tribution des intensités uniquement dans les inter- 
valles de fréquence de — 22 cmt à — 150 em! et 
de 10em-1 à 20cm t du centre de la raie. Les 
limites de ces intervalles se trouvant plus près du 
centre de la raie, sont définies par la limitation 
d'application de l’approximation W. K. B. Les 
limites extérieures des intervalles sont condition- 
nées par les hypothèses simplificatrices faites pour 
le calcul des intégrales de recouvrement. 


2 | 1 [ 
20 410 O -10 -20 -30 -40 -50 -60 -70 -80 -9D 


Le 


-100 -200 -300 
F1G. 5. — Raïe du mercure élargie par l’argon 

ES I = EC OC QE 12 Ame 
a  ! ourhe expérimentale d’après Michels. 
distribution de la dispersion. 
o points calculés par la méthode W. K. B. 
& points calculés numériquement. 


De l’estimation de l'influence de la largeur instru- 
mentale, de l’élargissement de Doppler et des colli- 
sions multiples sur la forme de la raie observée, il 
résulte que, dans la première approximation, on 
peut négliger ces influences. 
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Conclusions. — La méthode de calcul présentée, . 
basée sur l’application de la solution approximative » 
de l’équation de Schrodinger, applicable pour tout. 
l'intervalle de la variable indépendante, et sur les. 
calculs numériques des intégrales de recouvrement, 
a permis de déterminer la distribution des inten- 
sités d’après la formule générale de la théorie de 
Jablonski, sans avoir fait de simplifications limi-” 
tantes. On a exprimé l’opinion que[4]: “ It is clear 
that in any line broadening theory, a treatment : 
which considers only a single perturbing atom and 
then averges over all transitions of this system 
will always yield a line form which diverges at the 
center.” En effet, l'expression (3) donnant la dis- 
tribution asymptotique des intensités dans l’appro- 
ximation W. K. B. devient infinie pour le centre 
de la raie. Cela ne provient pas du fait que l’on a 
pris en considération des collisions simples, mais 
est relié aux approximations appliquées lors de 
l’établissement de cette formule. Par contre, l’appli-. 
cation de la fonction de Langer et les calculs numé- 
riques précis des intégrales de recouvrement, con- 
duisent à des valeurs finies, s’accordant avec l’expé- 
rience, aux fréquences pour lesquelles la formule (3) 
prend des valeurs très grandes tendant vers l’infini. 
Ces mêmes hypothèses simplificatrices font que la 
formule (2) ne donne des résultats corrects que pour 
un certain intervalle de la raie élargie. En appli- 
quant la formule générale (1), on obtient l’accord 
avec l’expérience, au delà des limites de ces inter- 
valles. 

L’auteur tient à exprimer sa profonde gratitude 
au P' A. Jablonski pour sa sollicitude, ses conseils 
éclairés et les nombreuses discussions qu’il a bien 
voulu lui accorder pendant ses recherches. 


Manuscrit reçu le 27 octobre 1960. 
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APPLICATION DE LA TRANSFORMATION DE FOLDY-WOUTHUYSEN 
A QUELQUES CAS D'INTERACTIONS ENTRE ÉLECTRON ET CHAMP ÉLECTROMAGNÉTIQUE 


Par M. BAKTAVATSALOU, 


Institut Henri-Poincaré, Paris. 


Résumé. — En utilisant la même méthode que Foldy-Wouthuysen, nous établissons d’abord la 


formule générale de l’hamiltonien de Dirac débarrassé de ses matrices impaires à l’approximation 
(1 fm), puis nous l’appliquons aux trois cas d'interaction : a) champ magnétique statique ; 
b) champ électrique statique ; c) champ électromagnétique variable avec le temps et obéissant 


aux équations de Maxvwell-Lorentz. 


Nous cherchons enfin à interpréter physiquement les différents termes du résultat obtenu. 


Abstract. — À general expression to the 4th approximation (1 /m)* for the Dirac Hamiltonian 
without odd matrices is given, using the Foldy-Wouthuysen method. 

The interaction is then discussed in the following three cases : a) Static magnetic field ; b) Static 
electric field ; e) Electromagnetic field varying with time and satisfying Maxwell-Lorentz equations. 

The physical interpretation of the different terms obtained from the final result is also discussed. 


1. Introduction. -— Les particules de spin 1 12 


(particules de Dirac avec à — c — 1) sont décrites 
par des fonctions d’onde Ÿ à 4 composantes. Pour 
des raisons notamment d’interprétation physique 
des opérateurs de Dirac, Foldy et Wouthuysen [1] 
ont élaboré une transformation permettant de dé- 
composer la fonction d’onde de Dirac et l’équation 
correspondante à l’énergie positive et celles corres- 
pondant à l’énergie négative. Cette transformation 
est essentiellement basée sur la séparation des 
matrices en matrices « paires » et «impaires ». Dans 
le cas où les particules de Dirac sont libres (sans 
être soumises à l’action d’aucun champ extérieur) 
la transformation est simple et la méthode est 
élégante [1], [2]. 

Quand les particules de Dirac sont soumises à 
des champs extérieurs (exprimés classiquement par 
leurs intensités et leurs potentiels), il ne semble pas 
exister une transformation canonique simple et on 
a alors recours à des approximations successives. 
C’est ainsi que Foldy et Wouthuysen [1] ont étudié 
le problème de l'interaction de la particule de 
Dirac avec le champ électromagnétique statique à 
l’approximation (1 /m)°. 

Utilisant la méthode de Foldy et Wouthuysen, 
nous étudions trois cas d'interaction (deux premiers 
cas simples et le troisième très général) à l’approxi- 
mation (1/m<). Nous essayons ainsi, en poussant 
les calculs à des approximations supérieures, de 
trouver des caractéristiques et des propriétés ma- 
gnétiques, électriques et électromagnétiques des 
particules de spin 1/2. 

Rappelons dans cette introduction que le terme 
d'interaction obéit à certaines conditions. Ainsi : 
le terme d’interaction est linéaire par rapport aux 
potentiels (cas des champs faibles) ; le terme d’in- 
teraction ne disparait pas en même temps que le 


moment p ; les distributions de charge et de cou- 
rant sont suffisamment localisées pour que leurs 
interactions avec les termes du champ puissent être 
exprimées, par des séries, en fonction des potentiels 
du champ et des dérivées successives de ces poten- 
tiels par rapport à la position de la particule. 

Ces conditions nous permettent en particulier de 
calculer les différents termes d’interaction de la 
série ; les coefficients de cette série caractérisent 
les propriétés magnétiques, électriques et électro- 
magnétiques de la particule. 


2. Formule générale de l’Hamiltonien de Dirac 
débarrassé des matrices impaires à l’approximation 
de (1/m). — La fonction d’onde de la théorie de 
Dirac est une matrice colonne à 4 composantes 


elle satisfait à l’équation (À — c — 


1 
Fo Ÿ (pi A; y (, 1) 


) 
dt : 


où l'opérateur hamiltonien Æ est une matrice her- 
mitienne à 4 lignes et à 4 colonnes : 


AD | ; 
H = fm + à p + terme d'interaction. 


+ 
Nous prenons pour « et B, la représentation de 
Dirae : 


y 


160 
(HAN EC RARE Er 10000 0 
DROIT AIO El 018780 0 

Ua —= 3 Es 

è LU T0 0 Net (DR EAN eu 
LR ENT ER | NE ARR CT 


Un opérateur matrice est dit « pair » [3] s’il n’a 
pas d’éléments reliant les composantes (4, d.) et 
(Ÿs Ya) de la fonction Ÿ de Lirac tandis qu’un 
opérateur «impair » est celui qui a de tels éléments 
de matrice. 8 commutant avec la matrice paire et 
anticommutant avec la matrice impaire, tout opé- 
rateur À de Dirac s’écrit : 


A=SfA+6AR)+54—648} 


où les crochets du 2° membre représentent la partie 
paire et la partie impaire de l'opérateur À. 

Soit alors S un opérateur hermitien dépendant 
explicitement du temps ; faisons subir à 4 la trans- 
formation canonique 


d — eis 


la fonction Ÿ' satisfaisant à l’équation, 


Pen) = Hy(S) 


> 
dt 
où 


H' = eis He -1S — es à ets, 
L 
Sachant que : 


2 
| 


es He—is = H + À [SH] + a LS, ISA + 


eis 3 e—is POS. | 
dt 


L + 
RCI ET 


Pen 
(2) devient alors : 


+ils is #55] 


OI 


Quand la particule de Dirac est placée dans un 
champ extérieur, l’hamiltonien peut toujours se 
mettre sous la forme : 


H=f$fm+e+oO 


(3) 


(4) 


où € (even), O (odd) sont respectivement les opéra- 
teurs pairs et impairs pouvant dépendre explici- 
tement du temps. Nous supposons en outre que € 
et O ne peuvent être d’ordre inférieur à (1/m)°. 
Considérons alors la transformation canonique 
engendrée par 
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cet opérateur S$ étant considéré comme petit, nous 
pouvons utiliser le développement en série (3). 
Nous calculons ainsi l’hamultonien transformé A 
de (3) dont nous retiendrons seulement les termes 


dans l’approximation (1/m)*. Pour ceci nous cal- . 
culons d’abord les différents commutateurs sui-. 


vants : 


2]: 
» dt , 


am? m? 
EU | 
(SLS, LSAN) = —, 08 — À, BIO[O[Oe]I] — 5 80! ; 
DO II 1000 1IE 


: 1 
gone lOIOIOIOeN] + 273 05; 


PE ee; 


(SES LSTSLS AT & 2 05. 


[SES TS (SAN) = 5 80% + 


La transformée de (3) s’écrira alors : 


2 Ke 1 + 800 

“ A nd 
t d0 41 1 

|| el a (O0 30" 


_ e F [0 [o Sl Te PIOLO(O EN) 


CRE i d0 
— Em PO F SSüm [° [o [o nl 
l 
$ 384m* 
Cette expression peut aussi se mettre sous la 
forme suivante : 


L'on NAS 1 2 1 3 
HONTE 2m ou Fra 
1 4 1 5 

pi rte Fine 


1 il 
| (1! 21)m GE Ci (21 22)mè 


Le As B[O[O[Oe]]] + ne [O[O[O[O<]]]] + … 


1 d0 1 


Fame een 


D? 


N°3 
— rase à [0 [0471] 


+ room [0 Lo [ol] + “ 4 


Cette méthode permet ainsi d'obtenir l’hamil- 
tonien de Dirac débarrassé des matrices impaires à 
une approximation arbitraire de 1 /m. Il suffit pour 
cela qu’on fasse une série de Sue ue dont 


chacune soit engendrée par S — — —  BX (terme 


impair de l’hamiltonien). Nous RE cette 
méthode à diverses types d'interactions. 


3. Particule de Dirac placée dans un champ 
magnétique pou JC k (2) dérivant d’un poten- 


tiel vecteur 4 . 
+ — 
& —= rot À. 


L'hamiltonien d'interaction s’écrit : 


de + 
ra (Ut) 

Suivant la méthode indiquée plus haut, nous 

effectuons d’abord la transformation engendrée par, 

>> 


Tr pe te 


$ À) 
LOTS es 


nous obtenons ainsi, (en remarquant que € — 0, 
> [> — 
LE É — € ie 


m8f(p—e4) ect 


môMES 
_. 2 (Ne 20) EEE p4 Peer) 
D hs 
lp ei) per L'EAU Fe 


En continuant ainsi jusqu’à l’élimination com- 
plète des matrices impaires à l’approximation de 
(1/m)#, nous obtenons finalement, 


Hu = êm+ 6(p—e A) ps 


. 21 


La série ainsi obtenue est convergente et il est 
facile de montrer que 


—>\2 > —> |2 
— A) —ecK 


Hu = Bi me + (CAEN) HR &ei} 
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En prenant 8 — 1 correspondant aux énergies 
positives, l’hamiltonien de Dirac débarrassé des 
matrices impaires à tout ordre de 1 /m, est : 


Hu = m? + NRA) es ee 


Ce résultat simple peut aussi être obtenu par une 
méthode élégante due à Case [4]. En effet cet 
auteur montre que si un opérateur O se met sous 
la forme : 


O = p10; + p30 


où les opérateurs O, et ©, commutent entre eux et 
commutent avec les matrices p, et p4 introduites 


par Dirac, 
> — 


&X — © O1; B = p3 


la transformée de © aura la forme : 
(g\ —_ p(02}1/2, 


L'application de cette méthode à notre premier 
cas d'interaction lui donne le même résultat que 
nous : 

1/2 


Hy=pim+|(p-e4) ef à 


Malgré l’élégance et la simplicité de la méthode 
de Case, il ne semble pas qu’elle soit susceptible de 
généralisation. Ainsi nous continuerons à appliquer 
la transformation de Foldy et Wouthuysen dans 
les deux cas d’interaction suivants. 


4, Particule de Dirac placée dans un champ 
> (> 
électrique statique E (pal dérivant d’un potentiel 
ES 
scalaire @ \æx). 
— > 
= — V 9. 


L’hamiltonien d'interaction s’écrit : 


ie 
H = fm + «x p — ce. 


Avec, 
; es 
S=— ee 
>> 
nous bte let 400 x p}. 
HAS 1-72 
re 
>(—> | 
er UE) 
7 8m? 
1 ire 
7 3m Por en L 
>> >> ds > 
te paie Up HN; F) pero tell 
4 (l die rs 
MT MI TIT OISE 
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pe VA Mare | PAPE Dans ce système d'unités nous avons posé c — 1, 
ses À AG div £)—( avé) \ | et les on magnétique et électrique sont 


cie VE p)p+2ievlaiv£)p+ a V(ÉPNp A0) 
+teV\E pJp+ieVVVE /p+eN\EpAp AS? y e—1. pet J sont respectivement les densités 


de charge et de courant. 


FA CMUN ERA S L’hamiltonien d'interaction est : 
>> > 
PAË—E Ap)}} 

Cr Er clerc 
+ 30mi p\x p}, 


— k —1di FA 
ip" de Mere ee pre ep te on 


La première transformée de Æ s’obtient avec 


Avec ; ee Ve 
Sn SA UD eh) 
LD —+ < 2 BA 2m | 
te) 

Nes ri) E = — eo. 
et ainsi de suite. nous continuons jusqu’à l’élimi- 
nation complète des matrices impaires à l’approxi- 
mation de (1/m)*. Nous obtenons finalement, 


Nous obtenons, 


; >> 
Hi; = Êm — eg + 2e B(a.v e) 
’ 4 2 m 
Hg = Êm—eo Dire Bp 


nn 
NÉ etpAË Er) | me 
Les RDS ea mi 
; 7 8m?| de? 5 dt dt r) 
ren Eee Sr er 
8m 8m? 12m — fev Vo—ee(pAVe—VeA pr) 
VF»); >> >> (2 fe Fa 
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+ 2%0(E en Es Cn-6) —e6XR 
F1 4064 si ë >> > 3m 
+ p.Vis PANE Le eV VdivE 2 


p—e) eo [0] 2 8 trop 
8 Ar ec k| 
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FPE eAŸ Re Fa a A 


> > (>> 
PE) or v lai Er Lusev Ép) HNe) 


+ 2e Ÿ (iv Ê) És ne) 


— Ph div FAERIA S A Fe A . ds + For 
_ ; où, 
5. Particule de Dirac placée dans un champ 
électromagnétique variable avec le temps et obéis- | Fo 1 | + +» ve ms 


: ; ù 
sant aux équations de Maxwell-Lorentz : TT PRE NS SR TA RSR Le 
— 
Das = > dE > me AO do > k 4 >> 
J— rot À, TORRES kr J — tie a (2 P = JE r) âe® (2 ee): xp 
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div XÆ — 0, >> >> > > 
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OI EE 
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SERA NE te mon Core) 
_ > > > 
oovaiv(Vo Ar) LépiVV 0 
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D Cain pare a)drare ce(s 12) 
Aller lea) de lv enr 2) 
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—divA24Ve+ io (À AV 8)! 


> x re 
+ sn A10(E A À). (10) 
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En continuant notre calcul avec 
e [> À 

=mtant GT) (re) 

> > > + 
Rp pee 
3m? 


et ainsi de suite jusqu’à l'élimination complète 
des matrices impaires à l’approximation de (1/m)* 
nous obtenons finalement : 


, re 
Him = Êm—ev + 8 (+ 24h — € a 


9 d° 
23 


PER 


SRE 
+ = div £— PE (PNA) 


a RAS AVe—Ve Ap)i 


+ et Lef(p—e a) = R 
A\ 
mel rca) sn p—eA)(É +24) 


nette +e4) A(É +2 se 


A 
HÉNEEENE] 
é Rs Le Lo [o || À in [O[O[O[Oej]}]. 


Les commutateurs figurant dans le second 
membre de Hzm Sont donnés par (9) et (10). 


(11) 


6. Interprétation physique des différents termes 
du résultat. — Notons d’abord que les expressions 
deviennent de plus en plus compliquées quand 
l’approximation augmente. Notons aussi que dans 
le troisième cas général nous retrouvons les termes 
des expressions de #3 et Hy. Nous avons done à 
interpréter les différents termes de Hzym- 

L'interprétation physique des termes de Hz1m 
peut être obtenue en remarquant que la transfor- 
mation de Foldy et Wouthuysen est une trans- 
formation non locale. On transforme, en effet, une 
fonction 4 en une fonction d’ qui s'obtient en 
effectuant une moyenne sur la valeur de 4 dans un 


ed 
volume entourant x dont les dimensions sont de 
l’ordre de la longueur d’onde de Compton : {/m. 
Ainsi donc dans la nouvelle représentation l’opé- 


rateur position n’est plus æ, mais un opérareur de 


= 
« position moyenne » X. Alors l'interaction d’une 
particule de Dirac et d’un champ extérieur va 
s’exprimer en fonction des grandeurs du champ à 
sa position moyenne. C’est ainsi que nous pouvons 
expliquer les termes de correction qui interviennent 
dans les expressions de Hy, Hy, Hum VU que la 
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particule occupe, non un point, mais une ie 
_ —\ 2 
p—eÀA 
2m + 
est l’opérateur hamiltonien classique de l’équation 
de Schrodinger ; le terme div Æ est dû à la distri- 
bution de la charge dans la région occupée par 
l’électron ; le coefficient e/8m°? est la correction 
bien connue de Darwin et qui est responsable du 
déplacement du niveau S dans l’atome d’hydro- 
SES 
gène ; le terme e/2m he x) représente l’énergie du 
couplage 252 RSR magnétique ; le terme 


het A Le Ë À 4) peut être interprété 


région de l’espace : l’expression 


N°3 
comme le terme d'énergie provenant du moment 


électrique associé à l’opérateur e/4m? CG Â mi 
p*I8m* est le terme généralement ajouté à l’équa- 
ba de Schrodinger pour rendre compte de la varia- 
tion de la masse avec la vitesse. 
Nous pensons que des recherches ultérieures 
pourront préciser l'interprétation physique des 
termes 


e2 dt Le 
Ë AA), 
Emi © 


>> > 
( À 4) 
qui nous paraissent intéressants. 


Manuscrit reçu le 23 septembre 1960. 
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INTERACTIONS ENTRE GRAVITON ET ÉLECTRON EN THÉORIE FONCTIONNELLE | 


Par PHAM XUAN YEM, 


Physique mathématique, Institut Henri-Poincaré, Paris. 


Résumé. — Les interactions entre un électron ét un graviton se manifestent par les termes de 
couplage qui figurent dans leurs équations. L’intervention de l’électron dans les équations de la 
gravitation est due à son tenseur d’impulsion-énergie. En revanche, les potentiels de gravitation 
associés au graviton agissent sur le mouvement de l’électron. 


Abstract. — The interactions between an electron and a graviton are expressed by the coupling 


terms in the equations of those two types of particles. 


The action of the electron in the 


gravitational equations is represented in terms of its energy-momentum tensor. Conversely, 
the potentials of the gravitational field associated with the graviton will react upon the movement 


of the electron. 


Dans un article paru récemment [1] nous avons 
proposé un modèle d'interaction entre le rayon- 
nement et la matière. Nous nous proposons dans 
cette seconde partie d’étudier un autre modèle 
d'interaction : celui de l’électron et du graviton. 
Nous montrerons que leurs équations décrivent la 
loi de gravitation microscopique. En passant à la 
moyenne d’un grand nombre de gravitons, on 
trouve que l’équation de l’onde moyenne donne 
exactement la loi macroscopique de gravitation 
faible comme en relativité générale. 

En électrodynamique quantique, on suppose 
qu’un électron est toujours entouré d’un nombre 
variable de photons qu’il émet et qu’il absorbe vir- 
tuellement. Il est donc légitime de supposer que 
l’électron, dont la masse a une origine non électro- 
magnétique, est toujours entouré d’un nombre va- 
riable de gravitons qui sont responsables des phé- 
nomènes de gravitation microscopiques. Avant 
d’aborder le problème d’interaction entre l’électron 
et un grand nombre de gravitons, examinons tout 
d’abord comment se présente l’interaction entre un 
électron et un graviton en théorie fonctionnelle. 
Soit z, l'onde physique représentant l’électron et w 
l’onde du graviton ; w et w, obéissent respecti- 
vement aux équations : 


| Le Ue — Qe (I) 
Leg — Qp. 

he : Neon ar r(Ug, Ue) ee Qc IT 

Nous écrivons : die 0. (I) 

où Leo Ue Le co (III) 


—— 


Lgo Ug — Qgo 


désignent respectivement les équations de l’électron 
libre et du graviton libre. 

Les interactions entre les deux corpuscules se 
manifestent par les termes analytiques r(we, we) et 
Ru, ue) qui figurent dans les équations (II). On 
peut garder Iss parties linéaires en w. ou en w, des 
termes d’interaction aux premiers membres des 
équations (11) et on transporte ainsi les parties non 


linéaires aux séconds membres, ce qui revient à 
remplacer les termes non linéaires Q, par Q avec 


Qo F Q. 


Lo Fe — 0 est l’équation de Dirac de l’électron 


mo 


het 


Fa VU + imo. 

Lo Fz — 0 est l’équation du graviton libre, cor- 
puscule de spin 2 obtenu par la fusion de quatre 
corpuscules de spin 1/2. 

Leo Te — 0 se divise en deux groupes d’équa- 
tions : les équations d’évolution et les équations de 
condition. 

Lo peut être mis sous la forme de l’opérateur de 
Dirac : (2) 


ITONRE RE 
Leo ES Vo mile 


m et u sont des inverses de la longueur d’onde de 
Compton : 


Fi —\ Me CA, Due of? 


Remarquons que les solutions des équations : 


{ Leo Le = 0 
À So Le = 0 


doivent vérifier l’équation de Klein-Gordon (sans 
que la réciproque soit nécessairement vraie) (l’équa- 
tion de Dirac se déduit de l’équation de Klein- 
Gordon par l’introduction des matrices «). 

Nous avons : 
PE tree 
EFe = 64 Vo 


En théorie non linéaire, l’onde physique x obéit 
à une équation ©.u — Q, la solution de cette équa- 
tion ne vérifie plus l’équation de Klein-Gordon 
mais plutôt une sorte d’équation analogue à celle 
de Klein-Gordon et non linéaire. Nous devons nous 
attendre à une équation de la forme suivante : 


(1) Dj ue = mu + F,'Q.) — F;tr) 
(2) [ue = bug + Fa(Qs) — F3(R) 


(IV) 
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où F,, F, sont des fonctions avec F(0) — 0. (Condi- 
tion de raccordement avec l’équation de Klein- 
Gordon.) | 

Comme w vérifie l’équation (IV, 2), les coeffi- 
cients du développement de # sur la base des 
seize matrices T\ obéissent aussi à des équations 
de même forme. 

Les équations (IV, 2) peuvent être considérées 
comme équations de propagation de l’onde w (à 
cause du dalembertien [ ]). Notre raisonnement 
sera confirmé par les calculs. Surtout dans une 
théorie linéaire du graviton comme celle de 
Mme Tonnelat [9] il sera inutile d’effectuer les cal- 
culs pour aboutir à l’équation de propagation 
[lu — u? uw. Mais dans une théorie non linéaire, 
la forme de F(Q) n’est obtenue que par des calculs 
détaillés. En remplaçant les ue, R(ue, Ua), Q{(Ue, Us, V) 
par leurs développements sur la base des ma- 
trices T, citées plus haut, on constate que les coef- 
ficients de ce développement se répartissent en six 
groupes d'équations [2]. Seul le groupe S, retient 
notre attention car il nous conduit aux lois de gra- 
vitation microscopique. 

Soit DA, Q1, RA les coefficients du dévelop- 
pement de w, Q, R; sur la matrice de base T, ; 
posons ya = Da + Q41 — R,, le groupe S, s’ex- 
plicite en [2]: 


du Dive) — dy Duo) eat ILVI0 
d Droulv = — Xuv (V) 
du Drpolv — d Dioolu — X{([uv] (pol)- 


Ce sont les équations du graviton et nous devons 
nous attendre à ce que ces équations donnent les 
lois de gravitation microphysique. Pour rattacher 
le corpuscule de spin 2 à la gravitation il nous faut 
définir un tenseur symétrique de rang 2 comme 
tenseur des potentiels de gravitation, qui lui-même 
représente la métrique de l’espace microphysique. 
(Comme en relativité générale, le tenseur des poten- 
tiels de gravitation n’est autre chose que la mé- 
trique de l'Univers.) 

D’après les raisonnements de M. Destouches [2], 
nous posons : 


guy = A(UDueltven + Quueitven + Qruv). 


Nous avons d’une part le graviton qui se propage 
dans un univers pseudo-euclidien de Minkowski, 
d’autre part à son onde physique w se trouve 
associé un tenseur g,,. Nous supposons que le 
champ créé par le graviton est un champ faible — 
ce qui est toujours acceptable au point de vue 
physique. — Le tenseur fondamental g,, ne doit 
pas s’écarter trop de la métrique minkowskienne, 
elle est de la forme : 

uv = uv + Ruv 
avec uv = (—1, —1, —1, +1) 


et },, est petit. 
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Comme il est dit plus haut, tout champ associé à : 
l’onde w, vérifie une équation de propagation de . 
la forme : 


CT Ug = ue? ug + FalQs) — F;(R) 


le champ g,, obéit par conséquent à : 


C] guv = LU? guv + F:(Qg) — F,(R) (VI) 
avec : 
F:(Qg) = — Ay? Qu) a AUD Qreuv + dy Qrele) 
nn Adu(20d6 Q(ew — dy Qpo))- (VT’) 


L’équation (VI) est une équation rigoureuse qui 
dérive du groupe ($,). C’est une équation non li- 
néaire en g,,. Elle constitue la loi de gravitation 
fournie par la théorie non linéaire du graviton. 
Nous constatons que la façon la plus naturelle 
pour rattacher la loi de gravitation due au graviton 
à la relativité générale est de poser : 


2 = — pu? 
ES fee) (vr') 
Re, = F0) 
L’équation (VI) du graviton devient : 
ET gN RE EE Pa ee Tv (VII) 


C’est la non linéarité de cette équation (à cause 
du terme ®,,) qui introduit des quantités ayant 
même forme qu’en relativité générale. Nous pou- 
vons écrire la loi de gravitation sous la forme : 


87rk 


; L 
Guv 3 uv G= — Xgpy + SV Tuv. (VIII) 


Les équations (VII) et (VIII) sont équivalentes. 


REMARQUES. — 1) Au lieu du tenseur g,, on peut 
considérer le tenseur F,, — Rav £us k, et l’'équa- 


tion (VII) devient : 


(IX) 


2) Nous avons posé plus haut (VI”), d’une façon 
arbitraire 

16 zxk 

FAQ) = 


c4 


Luv et FR) = — 


pour la seule raison de rattachement de l’équation 
du graviton à la relativité générale. Mais au point 
de vue mathématique, il faut montrer que F,(@) 
et F,(R) se comportent comme des tenseurs dans 
l’espace de Minkowski. Or ceci est vrai puisque 
FA@:) et FA(R) se transforment exactement 
comme Quws do Qrouy --. Comme on le constate 
aisément dans la relation (VI'). 

Puisque F,(R) est un tenseur, nous le désignons 
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par le tenseur d’impulsion-énergie de l’électron afin 
de raccorder l’équation du graviton à la relativité 
générale. 

De même pour F,(Q:;) — (16 zkfet) &,, nous 
pouvons attribuer à &,, la notion de tenseur d’im- 
pulsion-énergie propre du graviton. On y voit une 
analogie frappante entre le photon et le graviton : 
l’électron en mouvement crée un champ électro- 
magnétique, les photons associés à ce champ réa- 

_gissent sur le mouvement de l’électron, ce qui 
donne lieu à des forces « damping », des self- 
énergies ainsi que la polarisation du vide, etc. Les 
termes non linéaires en sont la cause. De même, un 
électron est entouré d’un grand nombre de gra- 
vitons, ces derniers créent un champ gravifique qui 
réagit sur le mouvement de l’électron et les termes 
non linéaires Q, en sont responsables. C’est pour 
cela que nous considérons F,(Q@:) à une constante 
près comme le tenseur d’impulsion-énergie propre 
du graviton. 

Revenons maintenant au tenseur d’impulsion- 
énergie 0,, de l’électron [3]. On connaît bien son 
expression en fonction de l’onde w. 


ACER L 
Ouv = 9; (ue Yu dv Ue — dy Ue Yu Ue) 


0,,n’est pas symétrique mais on peut toujours le 
symétriser. En cas général, la symétrisation con- 
siste à ajouter à 0,, une certaine fonction, mais 
dans le cas de l’électron, Pauli [3] a montré qu’on 
obtient un tenseur T,, symétrique et équivalent 
à 0,, en considérant : 


Tuv = (1/2) (Buv + vu). 


A partir du tenseur 7, on définit le tenseur con- 
travariant 74=— 718 4% 7°. et puis le quadrivecteur 
impulsion-énergie P4, avec : 


Pv = [Tue d° æ. 


La 4€ composante de P# retient notre attention 
car elle intervient dans le calcul si nous nous bor- 
nons à considérer seulement l’équation : 


167k 46Tk 
MY 2 ges + ci T ya 2e Taa 
En effet : 
"AG HU: Dite 
Tas Le Cdt dr 
Mais 
= (av — Tai) ue + Que) — rl 4) 
dTa \ 


Tia = uë(— inc œ + me Ca) ue + D (Q(Ue) — rl Me)] 
Avec l’approximation suivante Q = 0,r = 0 et 
la vitesse de l’électron est petite (cu; & ue € 1), 


nous avons : 
T'aa © Ue Me C? O4 Ue 
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et l’équation [1 d,, = — RE 


Ta donne la loi de 


gravitation de Newton [4]. 
Équation de l’onde moyenne. — Partant de 
l’équation non linéaire du corpuscule de spin 2 : 
Leo Ug + Rue, Ug) = Qg 


nous arrivons à établir la loi de gravitation non 
linéaire liée au graviton. Cette loi est : 


äl 8Tk 
Guv 3 uv G = — Xguy + a Tu 
ou bien : 
__167k 4167rk 


ANNEE € 


Es] uv = 2)guv n Tv 


c4 

Nous allons montrer qu’en passant à la moyenne, 
l’équation (IX) devient une loi macroscopique de 
gravitation. En effet, considérons un système de 
plusieurs gravitons (V gravitons) dont l’onde 
moyenne est définie par [5]: ‘ 


1 
UM mr? Uj. 
J 


Le potentiel moyen gym associé à 
moyenne 4 est défini comme toujours par : 


l’onde 


Euvm = A(UDuotvenM + QueltvenM + Quww). 


Si les quatre conditions de M. Destouches sont 
remplies, um obéit à une équation linéaire (Qy 
tend vers zéro). De l’équation linéaire : 


Soum + R(um, ue) = 0 


nous déduisons l’équation en g,,m qui est aussi elle- 
même linéaire : 

167k É 
C1] Suvm = — 2Xguvm + DA Tv. (X) 


Le terme F(Qm) = ue 


avec Q; (puisque F,(0) = 0). 

C’est l’équation du champ moyen déduite d’une 
façon formelle de l'équation & uy = (. 

Nous allons montrer qu’elle s’identifie avec la 
loi macroscopique de gravitation. Pour cela nous 
devons supposer que le champ moyen du système 
est un champ gravifique faible. Or dans un champ 
faible en relativité générale, l’espace-temps est 
muni d’une métrique presque euclidienne : 


Tuvm tend vers zéro 


Suv — Yu ci Ruv avec Ru petit. 


Alors, si la condition suivante est remplie [4] 


d Pix dal 
K=1 ÔT% 


0 


. $ l 
le tenseur géométrique d’Einstein [4] G,, — 3 8, G 


est égal à = EE 
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La loi de gravitation en relativité générale : 


1 8Trk 
Guv LE Luv G = — Aguv + Ty Tu 
devient [4] : 
- 167rk 
CT Suv = — 2Xguv + ns Tv (XI) 


L’équation de gravitation faible est linéaire. 

Les équations (X) et (XI) sont identiques. Le 
passage du microscopique au macroscopique est 
donc résolu si l’on suppose que le champ est faible 
et si l’on remplace les grandeurs à valeurs com- 
plexes par les grandeurs à valeurs réelles. 


Équation de l’électron. — Quelle est l’équation 
d’un électron en interaction avec un graviton ? 
On a fait de nombreuses tentatives pour relier la 
microphysique à la relativité générale en consi- 
dérant que les corpuscules se meuvent dans un 
espace riemannien. On remplace les dérivés ordi- 
naires par les dérivées covariantes, les matrices y, 
vérifient les relations d’anticommutation 


Vu Yo Yo Ya = 2Buv 
Malheureusement ce processus présente beaucoup 
de difficultés formelles [6]. On peut se trouver tenté 
de supposer avec Moshinsky [7] que l’électron se 
meut toujours dans l’espace euclidien. L’inter- 
action du champ gravifique avec l’électron se mani- 
feste par le terme 71,, g**, terme analogue à l’inter- 
action 7, A* du champ électromagnétique avec 
l’électron. Dans ce dernier cas, on doit remplacer 
les opérateurs p* par p# +(efc) A". De même, si 
l’on prend pour l'interaction l'expression T,, g*’on 


doit remplacer [7] p“ par: pi guv PA Fbrié 


(notons que p*— p, en coordonnées orthonor- 
males). 

Mais la façon la plus simple pour obtenir l’équa- 
tion de l’électron plongé dans un champ gravifique 
est d’établir l’hamiltonien $ de ce corpuscule. Une 
fois connu $, l’équation de l’électron s’écrit : 


(s+ in) u = 


(postulat de la théorie fonctionnelle). 
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En théorie classique des champs, l’hamiltonien 
relativiste d’une particule libre est 


H, = Vp°ce + mc. 
Pour une particule plongée dans un champ gravi- 
fique de potentiel newtonien U, l’hamiltonien 


est : [8] 
H = /ue + me e*) (1 +2). 


En relativité générale, l’hamiltonien du cor- 
puscule plongé dans un champ gravifique dépend 
non seulement du potentiel scalaire g,, mais aussi 
du potentiel vecteur ga (t — 1, 2, 3). Mais on peut 
toujours choisir un certain repère particulier par 
rapport auquel les g,, sont nuls de sorte que l’ha- 
miltonien À peut s’écrire sous la forme (XTT) [8]. 

On constate aisément que dans le cas où la 
vitesse de l’électron est beaucoup plus petite que c, 
H diffère peu de l’hamiltonien non relativiste H1. 


(XII) 


p° 
FR == Ver LD 4 mU. 
En effet, 


Or nous savons que : 


84aa — — (1 de 2): donc H = V{(p? ce? + m? ci) (— gaa). 


En passant à la théorie quantique, l’hamilto- 
nien $% de l’électron devient : 


$ = (— iic a. + me c? «4) v= £a 


Cela nous permet d’écrire l’équation de l’électron 
dans un champ gravifique. 


[ À + (— irca.V + me C2?) V— sul Ur O}, 
L’interaction du graviton sur l’électron se mani- 


feste par le terme multiplicatif /— g,,, avec : 


NOUS AZ (BD(Ho1taen + Quuetsen + Qaa)- 
Manuscrit reçu le 11 juillet 1960, 
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ÉTUDE DU TYPE D'INVARIANT DE L'INTERACTION GAMOW-TELLER 
EN DÉSINTÉGRATION $- DE 6He 


Par TosniKko YUASA, 


Laboratoire de Physique Nucléaire d'Orsay, Faculté des Sciences de Paris. 


Résumé. — Le spectre énergétique des atomes de recul lors des désintégrations f— de He a été 
obtenu à l’aide d’une chambre à brouillard de Wilson à très basse pression. 
D’après l’analyse de ce spectre, le type d’invariant de l’interaction G-T pure ne nous semble 


pas être T', mais À. 


Abstract. — The energy spectrum of recoil atoms due to f—decay of $He has been obtained 


using a very low pressure cloud chamber. 


The analysis of this spectrum indicates that the type of invariant of the G-T interaction 


should not be T7, but À. 


I. Introduetion. — On sait qu’on peut déter- 
miner les invariants dans l’interaction $ par l’étude 
des corrélations angulaires $-neutrinos, sans utiliser 
des noyaux orientés. Théoriquement ce serait la 
meilleure méthode pour déterminer les rapports des 
constantes |Cx| (où X — $, V, T, À ou P) de 
l'interaction, 


J = > Cx(Ÿ5 Ox Ya) (Le Ox ds) + h.c 


étant le produit des champs quantifiés corres- 
pondant à l’émission ou à l’absorption des difié- 
rentes particules impliquées dans les processus, 
= f J(æ) d°x étant l’hamiltonien d'interaction. 


+ Depuis la découverte de la non conservation de la 
parité, on sait qu’il faut ajouter à J, le pseudo- 
invariant : 


J' = > Cx(Ÿ5 Ox Ÿn) (V8 Ox Vs Yv). 


- Mais l’expérience montrant que la polarisation 
longitudinale des $ a la valeur maximum v/c pré- 
vue par la théorie, il en résulte nécessairement [1] : 


Cy= Cr; Cs=—Cs 
Cie Cr Cr 
(Cp = — CP) 


On est donc toujours ramené à la détermination 
des |Cxl. 

Pour des noyaux non orientés, la loi de corré- 
lation B-neutrino reste la même qu'avant la non- 
conservation de la parité : 


1 + 2vfe cos 0 ... (2) 
. avec 


x= nf fol (ICrl?— 1C4/? 
+|f: 


= if: Fticer? IcrPTE TORCE |Cal®). 


2([Cr|?— (Csl] : 


En particulier pour une interaction G-T pure: 
1 
Xe AGREE (Cal FCPI TE ICAIE NS. (3) 


Pour chaque type d’invariant séparément, À 
prend les valeurs suivantes : 


T-., 44/8 
eriA 4/3 
CPE A : 
Frs dt 


Par conséquent si l’on observe la direction d’un 
rayon B émis et celle du noyau de recul associé, on 
doit déterminer la valeur de X et en conséquence le 
type d’invariant. Expérimentalement, 3 types de 
mesures ont été Jusqu’ici réalisées. 

19 Dans le triangle (fig. 1) traduisant la relation 


8 


HTC 


- 


Pe + 4 + R—0, on mesure À, c’est-à-dire l’énergie 
de l’atome de recul R?/2M, ps, c’est-à-dire l’énergie 
du 8, æ = (p$ + 1)42? q en résulte : g =W, —- W 
(en négligeant l’énergie de l’atome de recul), et on 
peut construire le triangle connaissant ses 
3 côtés [2]. 

En pratique on a des bandes d’énergie pour 
l’atome de recul et pour le 8, la bande pouvant 
s'étendre parfois à tout le spectre $. 

20 On mesure pg et 6,. Il y correspond en général 
deux triangles [3]. 
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30 Mesure du spectre de À ou R?/2M [4]. 

Lauterjung, Schimmer et Maier-Leibnitz [5] ont 
effectué un autre type des mesures : mesures de 
corrélation d'énergie entre les 8 de 

8Li — $8Be + B— + v 
et les deux particules « de $Be — 2(*He) qui suivent 
cette désintégration $ considérée comme G-T pure. 

Parmi ces résultats, il y a certaines divergences : 
les expériences sur Ne (Ridley [2]), sur $He 
(Rustad et al. [3], Szalay et al. [3] et le deuxième 
travail de Allen et al. [2]) indiquent, soit la prédo- 
minance de 7’, soit le mélange A-7, tandis que les 
expériences sur He (le premier travail de Allen 
et al. [2] et Pleasonton et al. [4]) et sur SLi (Lau- 
terjung et al. [5] et Barnes et al. [5]) indiquent À 
avec le mélange à 10 % au plus de 7. Les résul- 
tats de Rustad et Ruby avec $He ont été mis en 
doute par Wu et Schwarzschild en ce qui concerne 
leurs conditions expérimentales ; les corrections 
d’ailleurs semblent indiquer une prédominance de 
A et ces auteurs ont déclaré qu’il existe certaines 
incertitudes expérimentales (*). 

D'autre part il y a l’expérience de Goldhaber et 
al. [6] sur l’hélicité du neutrino dans une désin- 
tégration 8 (capture X) du type G-T qui indique 
que l'interaction G-T est À ; néanmoins 1l y a une 
incertitude de 15 % au plus. 

En résumé la plupart des résultats expéri- 
mentaux sont en faveur de À mais il n’y a aucun 
résultat démontrant que l’invariant À figure seul, 
sans mélange de T. 

La première méthode, mesure des p4 et À, a en 
général un défaut : faible précision de la mesure 
de À, tandis que la deuxième méthode, mesure 
de pget 6,,a d’autres défauts : précision limitée des 
statistiques à cause des coïncidences fortuites et à 
cause de l’incertitude des mesures des angles 0,. 

Pour éviter ces défauts de la deuxième méthode 
Allen et al. [4] ont récemment étudié uniquement 
le spectre énergétique des atomes de recul pour un 
grand nombre d’éléments dont le type d'interaction 
est G-T ou F + (G-T). Ils en ont déduit le coef- 
ficient de corrélations angulaires À par l’intermé- 
diaire de l'expression P(w, w,) dw dw, (dw, : éner- 
gie cinétique de l’atome de recul) déduite de la 
formule pour la probabilité P(w, 0) dw dQ. Ils ont 
ainsi obtenu les résultats suivants : 


COEFFICIENT ForME 


Lo ; 
ISOTOPE D SCT TION DE D’INTER- 
- CORRÉ LATION ACTION 
22A  (6+) grande partie : F + 0,97 + 0,14 , PT où VA 
22Né (B+) Fet GT 0,00 + 0,08 ST ou FA 
22Ne (E—) G-T — 0,37 + 0,04 A 
2He (B—) G-T — 0,37 = 0,05 A 


(*) ‘* Post deadline paper ”’ the 1958 New York meetings 
of American Physical Society. 


JOURNAL DE PHYSIQUE 


Ages Cole 


4 


No 3! 
À 
Ces auteurs ont conclu que la forme d'interaction. 
pour G-T est À (T < 10 %) en comparant le 
spectre énergétique des atomes de recul avec ceux, 
théoriques pour À et T dans le domaine d’énergie, 
> 700 eV. Au-dessous de cette énergie ils ont 
observé un écart entre les points expérimentaux et. 
les courbes théoriques de À et T pour He et Ne. 
Ils l’expliquent par l'existence des atomes de recul 
à charge 2. | 
En admettant leurs hypothèses, nous avons 
effectué une correction approchée ; nous avons. 
supposé que, comme cas extrême, l’écart entre les 
points expérimentaux et la courbe théorique T est 
dû aux atomes de recul ayant une charge 2 et par 
conséquent ayant une énergie deux fois plus grande, 
et nous avons obtenu la courbe C représentée dans. 
la figure 2. 


N 


\s Er 
0 700 1400 eV 
FiG. 2. — Courbes théoriques et points expérimentaux 


pour la répartition énergétique des atomes de recul de $He 
(Allen et al. Phys. Rev., 1959, 116, 134, fig. 4) et courbe 
expérimentale corrigée par nous pour effet des atomes. 
de recul ayant une charge 2. 
À : invariant À. 
T : invariant 7. 
C : courbe corrigée par nous. 


La courbe C représente une estimation extrême 
pour être favorable à T. Cette courbe se trouve 
entre les courbes T' et À, ce qui signifie, d’après. 
cette correction, que l’interaction ne serait peut- 
être pas À pure, mais qu’elle n’est certainement pas 
T pure (*). 


(*) L'auteur remercie M. le Pr Allen qui lui a commu- 
niqué son opinion sur cette correction. Il juge que cette 
correction serait justifiée à condition que la probabilité 
d'existence des atomes de recul ayant une charge 2 tombe 
brusquement au-dessous de 700 eV ; autrement dit cette 
correction a une tendance à favoriser l'interaction T. Le 
Pr Allen l’a informé dans la même occasion que le Dr Ridley 
de Harwelle avait effectué les mesures du spectre des 
atomes de recul de He à l’aide du temps de vol et qu'il 
avait obtenu À = — 0,35 + 0,02. Mais le Pr Allen se 
demande si son spectre ne serait pas aussi déformé par des 
atomes de recul ayant une charge 2. 
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Pour éviter cette incertitude causée par les 


atomes de recul ayant une charge 2, Snell et al. [4] 


ont étudié le spectre des quantités de mouvement 


des atomes de recul de 6He avec un spectromètre à 


double focalisation (magnétique et électrique) en 
introduisant la source de 6He dans un tube conique 
dont le bout constitut le point de la source. Ils ont 
obtenu, pour ces deux espèces d’atomes de recul, un 
spectre concordant à la courbe théorique À 
(+ 5 % près). Ils ont trouvé que le rapport des 
atomes de recul ayant une charge 2 au total est 
9,5 + 0,5 % à 1 000 eV. 

L'emploi d’une chambre à brouillard de Wilson 
pour l’étude des corrélations angulaires B-atomes 
de recul de He nous a paru intéressant bien que 


cette méthode ait certains défauts, en particulier : 


diffusion des atomes de recul par le gaz, effet 
d'échange de charge, déformation des trajectoires 
dans le temps, faibles statistiques dues à la diffi- 
culté d'introduire une quantité suffisante de $He 
ayant une période de 0,83 sec dans la chambre, et 
à la diminution de l'efficacité de la chambre avec la 
diminution de pression de gaz, etc car cette 
méthode a des avantages que les autres n’ont pas : 
le point de désintégration est bien défini, par consé- 
quént la mesure des angles 6, est précise ; les 
atomes de recul ne subissent aucune perturbation 
éventuelle à cause de leur accélération par un 
champ extérieur ; et enfin on peut exclure certains 
bruits de fond par l’observation directe. 

C’est pourquoi nous avons envisagé en 1955, 
d'étudier les corrélations angulaires $B-atomes de 
recul de $He à l’aide d’une chambre à brouillard de 
Wilson construite dans ce but pouvant fonctionner 
soit en autocommande soit en synchronisation et 
avec de très basses pressions (jusqu’à 1 cm Hg) [7] 
à [10]. 

Au cours de nos préparatifs, Szalay et al. [3] ont 
publié les résultats de leur étude des corrélations 
angulaires B-atome de recul de $He, à l’aide d’une 
chambre à brouillard remplie de 

(H2 + H,0 + C,H;0H) 
à une pression totale de 20 cm Hg. Leurs résultats, 
jusqu’à présent, sont favorables à l'interaction T (1). 

Puisqu’il existe encore des désaccords concernant 
la forme d'interaction G-T, surtout pour ‘He, 
nous avons voulu reprendre ce problème à l’aide 
de la chambre à brouillard précitée à très basse 
pression ; en faisant déclencher la chambre par un 
rayon 8 émis par $He au sein de la chambre, et 
compté par un compteur à scintillations, on peut 
photographier la trajectoire de l’atome de recul 
associé à ce $. 

Malheureusement l’état actuel de l'intensité des 
neutrons produits par le faisceau focalisé de neu- 


(1) D’après une récente communication privée, une publi- 
cation de ces auteurs qui met en doute ces conclusions 
va paraître dans Acta Hungaria. 
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trons de 6,5 MeV émis par le cyclotron de notre 
laboratoire (équivalent à environ 300C de la source 
de neutrons Po + Be) ne nous permet pas de faire 
une telle expérience autocommandée par f, et nous 
sommes obligés pour le moment, de nous contenter 
d'obtenir un spectre énergétique des atomes de 
recul de $He, à l’aide de détentes synchronisées avec 
le fonctionnement du cyclotron. 

Nous n’avons pas encore obtenu le nombre suffi- 
sant d'événements pour faire une statistique satis- 
faisante. (Cette expérience sera reprise après une 
modification des conditions de focalisation permet- 
tant d’obtenir un faisceau de deutons plus intense.) 
Mais puisqu'il s’agit d’une mesure par chambre de 
Wilson des parcours des atomes de recul B avec la 
plus basse pression réalisée jusqu’à présent et que 
ces résultats préliminaires donnent déjà certaines 
indications, ilnous semble intéressant de les publier. 


IT. Expériences. — a) Source DE $He. — Une 
des difficultés expérimentales est d'obtenir et d’in- 
troduire une source intense de $He en temps suffi- 
samment court dans la chambre, étant donné que 
le rendement de l’extraction de He produit est 
faible et que la période de $He est courte. 


Paraffine 


_N 
de deutons 


vers chambre 
NÙ à brouillard 


Fenétres 


Te me: 


Nous avons utilisé les neutrons rapides produits 
par la réaction (d, n) : le faisceau de deutons foca- 
lisé à l’aide de deux lentilles électrostatiques tombe 
sur une cible de Be métallique. Les neutrons ainsi 
produits irradient des cibles de Be(OH), déposées 
sur des fragments de lames couvre-objet. L’inten- 
sité des neutrons est équivalente à une source d’en- 
viron de 300 C (Po + Be) dans le cas où le faisceau 
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de deutons irradiant la cible est d'environ 0,2 uA. 
Le dispositif schématique pour la production de 
la source est reproduit dans la figure 3. 


b) INTRODUCTION DE He DANS LA CHAMBRE A 
BROUILLARD. — Pour introduire 6He dans la 
chambre, nous avons procédé de la façon suivante : 

19 On évacue la chambre à brouillard, le réci- 
pient de Be(OH), et les canalisations qui les relient 
en ouvrant les deux électrovannes e. v. ete. v. IT 
(fig. 3) et un robinet à aiguille situé sur le couvercle 
de la chambre. 

20 On ferme ces électrovannes et on règle l’ouver- 
ture du robinet à aiguille. . | 

30 On bloque électroniquement le système de 
détente de la chambre et simultanément on établit 
le courant du déflecteur du cyclotron. 

40 On irradie Be(OH), un temps chronométré 
(— 3 sec) à l’aide d’une minuterie. Au bout de ce 
temps, le courant du déflecteur étant coupé, les 
deux électrovannes s’ouvrent introduisant He ordi- 
naire sous pression dans le récipient de Be(OH), et 
entraînant $He vers la chambre à brouillard à 
travers le robinet à aiguille. 

5° Au bout de 1 à 2 secondes contrôlées par la 
deuxième minuterie, les électrovannes se ferment 
et le système de détente de la chambre à brouillard 
est débloqué. 

50 La troisième minuterie règle l’attente et le 
déclenchement de la détente. 

70 On mesure, après chaque détente, la pression 
du gaz dans la chambre qui est réglée à environ 
13 cm Hg (11,25 cm Hg de He + 1,75 cm Hg de 
H,0). 


c) CHAMBRE A BROUILLARD. — La description de 
la chambre et les caractéristiques de fonction- 
nement à de très basses pressions ont été déjà 
publiées [7], [10]. 

Dans le présent travail, la pression initiale était 
environ 13 em Hg. Mais étant donné que le remplis- 
sage de He est instantané la pression initiale réelle 
de chaqué détente n’est pas très précise ; cette 
incertitude serait inférieure à 5 %. D’autre part la 
pression finale n’est pas calculable étant donné que 
la détente est intermédiaire entre l’adiabatique et 
l’isotherme. En admettant que la pression initiale 
soit 143 cm Hg (He + H,0), on obtient, comme 
pression finale, soit 8,2 cm Hg (7,0 cm Hg de 
He + 1,2 cm Hg de H,0 ; adiabatique) soit 
9,6 cm Hg (8,5 cm Hg de He + 1,1 cm Hg de H,0 ; 
isothermique) pour le degré de détente de 1,353. 
Avec cette pression et avec He + H,0, la détente 
est plus proche de l’isothermique [10]. 

Toutes les mesures de parcours sont normalisées 
à la pression initiele de 13 em Hg. 

On n’a pas baissé plus la pression dans la 
chambre à cause du faible rendement des phéno- 
mèênes. 
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Comme contrôle de la détente et comme réfé-. 
rence pour la position en profondeur des trajec-. 
toires des atomes de recul lors du dépouillement. 
sommaire, on a placé une source de Po à l’intérieur. 
de la chambre. On peut dépister la profondeur de la. 
trajectoire intéressée grâce aux trajectoires des. 
rayons « qui se trouvent en éventail dans la zone 
éclairée de la chambre. 

Les clichés stéréoscopiques ont été pris en utili- 
sant uu tube de décharge à Xe (Phillips PF900X) 
dont la durée du flash est prolongée à l’aide d’un 
self de 1 mH relié en série jusqu’à environ 107$ sec, 
pour augmenter la sensibilité photographique. On a. 
utilisé les films HPS d’llford de 65 mm de largeur 
pour avoir la plus grande précision possible. 


d) MESURE DES TRAJECTOIRES DES ATOMES DE, 
RECUL. — Les trajectoires des atomes de recul dont 
l'énergie maximum est 1,405 keV, étant très courte, 
il nous a fallu les mesurer avec une méthode aussi. 
précise que possible. Nous avons donc utilisé le” 
stéréocomparateur de Pülfrich qui donne des pré 
cisions mécaniques de 1/100 mm pour les coot-" 
données x et y ; 1/1 000 mm pour la parallaxe. 


III. Résultats expérimentaux. — Nous avons 
retenu environ 300 paires de clichés dans lesquels” 
nous avons observé au dernier tri fait à l’aide d’un 
simple stéréoscope, environ 160 cas susceptibles 
d’être des atomes de recul de $He. Maïs après un 
examen plus précis fait avec le stéréo-comparateur 
de Puülfrich, nous n’en avons retenu que 50. 

L'aspect de ces trajectoires des atomes de recul 
est assez caractéristique et distinct de celles des 
rayons Ô des « du Po qui se trouvent également 
dans la chambre. Les premières sont en général 
très noires, très courtes et très souvent droites, 
tandis que les dernières sont très sinueuses et leurs 
longueurs plus variées. 

Nous montrons une de ces trajectoires dans la 
photographie. Les clichés À et A’ montrent un 
atome de recul avec les agrandissements différents, 
tandis que le cliché B représente un rayon « du Po 
avec quelques rayons Ô. Même l’agrandissement 
de 112 fois ne permet pas de voir les ions séparés les 
uns des autres dans la trajectoire des atomes de 
recul. 

La répartition des parcours normalisés à la 
pression initiale de 13 em Hg (He + H,0) est 
montré dans la figure 4. 

Le parcours le plus court mesuré était 0,26 mm, 
et celui maximum était 5,8 mm dans la pression 
précitée. Pour déduire les parcours dans l’air saturé 
de H,0 à NPT, il faut multiplier par un facteur 
d'environ 0,03 (entre 0,027 et 0,032) c’est-à-dire 
que le parcours maximum observé correspond à un 
parcours de 0,17 mm dans l’air saturé de H,0 
à NPT, 

Pour une énergie aussi basse, la relation parcours- 
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énergie n’est pas très bien étudiée n1 théoriquement 
ni expérimentalement. 

Nous avons étudié [11] la relation parcours- 
énergie des protons projetés lors des désintégra- 
tions « de An et AcA dont l’énergie minimum 
mesurée était de 0,5 keV à l’aide d’une chanbre à 
brouillard à basse pression (1 cm Hg H,0). 

En se basant sur ces mesures nous pouvons esti- 
mer très approximativement le parcours d’un 
atome de recul de $He de 1,4 keV dans l’air saturé 
de H,0 à NPT qui était de 0,02 mm. 

Plus tard, Mills [12] a étudié la relation parcours- 


PA 


atome de recul 


Cuioné B. Chiché À. énergie des rayons & jusqu'à 100 keV en utilisant 

Un « du Po et quelques 5. Un atome de recul. les rayons « projetés par les neutrons à l’aide d’une 

(agrandissement =. (agrandissement — chambre à brouillard remplie de He + H,0 à une 
grandeur naturelle X 9,7). grandeur naturelle X 9,7). 


CLICHÉ A’. — Le même atome de recul 
(agrandissement — grandeur naturelle X 112). 


30+ W 


0 P 700 1400 eV 


F16. 5. — Courbe I : Spectre énergétique des atomes de 
recul de He obtenu en admettant que Ec parcours : 
20 a; : Composante supérieure, 
a; : composante inférieure, 
a, : Spectre obtenu en supprimant les points se trou- 
vant en dehors du point P. 
Courbe II : Spectre énergétique des atomes de recul 
de $He obtenu en admettant que Ec (parcours). 


10 F pression de 4,5 cm Hg. En extrapolant sa courbe on 

obtient très approximativement 0,03 mm comme 

parcours d’un atome de recul de $He de 1,4 keV 

dans l’air saturé de H,0 à NPT. | 

A partir de notre parcours maximum observé, on 

obtient le parcours de 0,17 mm dass l’air saturé 

TRE 2 3 4 Ë TA de H,0 à NPT comme nous l’avons dit plus haut. 

(11,25 cm Hg de He +1,75 cm Hg deH,0) Si l’on admet que le parcours maximum correspond 

Fi à l’énergie maximum, il est 7 fois plus grand que 

pic ire Répartition des parcours des atomes de recul lui estimé d’après [11] et 5 fois plus que celui 
“Bo réduire les parcours dans l’air saturé de H,0 estimé d’après [12]. 

à NPT, il faut multiplier par 0,03). Si l’on considère que les atomes de recul, ayant 
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un parcours plus long que 2,7 mm (fig. 4) qui 
correspond au point P dans le spectre énergétique 
(fig. 5) ont subi l’échange de charge au cours de 
leur trajectoire, et par conséquent, leur parcours 
est plus ou moins prolongé, le rapport du nouveau 
parcours maximum (0,08 mm dans l’air saturé de 
H,0 à NPT) et celui estimé d’après la courbe de 
Mills devient de 2,65. C’est une hypothèse possible 
et nous allons en tenir compte plus loin lors de 
l’analyse des résultats. 

Enfin l’examen d’une des paires de clichés pris 
par Szalay et al. (*) montre que, s’il s’agit bien 
d’une trajectoire d’atome de recul, son parcours est 
au moins 2,5 fois plus long que notre parcours 
maximum observé en le déduisant dans l'air 
à NPT. 

Puisque nous ne pouvons pas déduire l’énergie 
d’un atome de recul en partant de son parcours, 
nous les avons calculés en admettant les deux rela- 
tions approximatives parcours-énergie : 10 Ecc par- 
cours ; 20 Ecx (parcours)#3. Avec la première hypo- 
thèse et en admettant que le parcours maximum 
observé correspond à l’énergie maximum, 1,4 keV 
et en normalisant la surface de cette courbe à celle 
de la courbe T dans la figure 2, on obtient la 
courbe I dans la figure 5. 

Avec la deuxième hypothèse, on obtient la 
courbe IT dans la même figure. 


IV. Analyse des résultats expérimentaux. — A 
cause de la faible statistique, on ne peut pas encore 
conclure d’une façon précise. Mais nous avons 
essayé, pour savoir l’ordre de grandeur de précision 
nécessaire dans l’avenir, d'analyser nos résultats 
expérimentaux en admettant les différentes hypo- 
thèses possibles. 

Nos spectres I et II dans la figure 5 montrent 
qu'il y a deux parties assez distinctes : partie 
d’énergie supérieure et partie d’énergie inférieure. 
Nous supposons d’abord qu’il existe deux sortes 
d’atomes de recul : ceux ayant une charge 1 et 
ceux ayant une charge 2. 

Nous avons décomposé le spectre I de la figure 5 
en deux composantes : spectre a, et spectre a, et 
nous supposons que le spectre a, correspond aux 
atomes de recul ayant une charge 1 et spectre a, 
à ceux ayant une charge 2. On admet que l’énergie 
des atomes de recul ayant une charge 2 est deux 
fois plus grande que celle déduite de son parcours. 

La courbe an, dans la figure 6 représente le 
spectre a, dans la figure 5 dont la surface est nor- 
malisée à celle de T dans la même figure qui est 
rapportée de !a figure 2. 

Si l’on suit le même procédé pour le spectre a, 


(*) M. le Pr Szalay nous a prété quelques clichés, et nous 
tenons à l’en remercier. Il nous a également indiqué, après 
le dépôt de cet article, que M. Csikai, son collaborateur, 
avait réexaminé ce parcours maximum et qu’il le trouvait 
1,3 fois plus long que le nôtre. 
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mais d’abord en doublant l’énergie de chaque point. 
et ensuite en normalisant à la surface 7, on obtient. 
la courbe à, dont on pourrait constater la même. 
allure que celle de la courbe a; Mais plus éloignée 
de la courbe T.. | 

Ensuite, nous avons admis l’hypothèse du pro- 
longement de parcours de certains atomes de recul. 
à cause de l’échange de charge, et on élimine les 
atomes de recul ayant un parcours plus long que: 
2,7 mm correspondant au point P dans la courbe I. 
dans la figure 5 et on obtient la courbe a, dans la 
figure 6 en faisant correspondre le parcours de 
2,7 mm à l’énergie maximum de 1,4 keV. 
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F1G. 6. — Analyse du spectre I représenté dans la figure 5, 
et les courbes théoriques À et T de la figure 2. 
Courbe a; n : courbe «a; ; surface normalisée à 7. 
Courbe &;n: Courbe a, ; corrigée pour l'énergie et 
surface normalisée à T', 
Courbe &;,n : Courbe a, ; corrigée pour l’énergie et 
surface normalisée à T. 


En examinant ces trois courbes dans la figure 6, 
on constaterait que l’effet de la charge des atomes 
de recul déformerait considérablement la forme du 
spectre. Cependant aucune courbe n’a d’allure de 
la courbe théorique 7. | 

Quant à la courbe IT dans la figure 5, elle est 
prédominante dans la partie de faible énergie et en 
conséquence 1l est certain que ses comppsantes 
seront plus écartées de la courbe T que pour la 
courbe I. 


V. Conclusions. — D’après nos résultats préli- 
minaires concernant le spectre énergétique des 
atomes de recul de fHe, on peut constater que la 
forme du spectre est très sensible à leur charge et 
qu’on ne peut pas avoir de spectre suffisamment 
précis tant qu’on ne peut pas diminuer encore la . 


TT 
; 
| 
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pression dans la chambre à brouillard et augmenter 
la statistique. 

De toute façon l'examen de la forme du spectre 
énergétique des atomes de recul seule n’est pas 
suffisamment sensible pour déterminer le type 
d'interaction ; il faudrait effectuer, en plus, les 
corrélations angulaires directes. 

Néanmoins, on peut déjà dire que la forme du 
spectre expérimental ne concorde pas avec l’inter- 
action 7° ; elle se rapproche de l'interaction À. 
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DÉTERMINATION DE L'INTENSITÉ D'’OSCILLATEUR 
DES RAIES DE LA SÉRIE VERTE DE Cu,0O AUX BASSES TEMPÉRATURES 


Par J. B. GRUN, M. SIESKIND et S. NIKITINE, 


Laboratoire de Spectroscopie et d’Optique du Corps Solide, 
Institut de Physique, Université de Strasbourg. 


Résumé. — Étude spectrophotométrique de l’absorption de la cuprite au voisinage de la « série 
verte » aux très basses températures. On a pu évaluer le facteur d’intensité d’oscillateur de la 
raie n — 2 de la «série verte », le coefficient maximum d’absorption et la largeur de cette raie et 
comparer avec la raie correspondante de la «série jaune». 


Abstract. — The green series of absorption lines in Cu,0 crystals has been studied by spectro- 
photometrie methods at temperatures down to liquid helium. The oscillator strength, the 
maximum of the coefficient of absorption and the width of the first lines (n — 2) of the green 
series have been determinated and compared with the corresponding results for the yellow series. 


Le spectre d'absorption de la cuprite, aux très 
basses températures, se compose dans l’ordredes 
longueurs d’onde décroissantes, de « l’absorption 
continue rouge », puis de la « série jaune » que suit 
« l'absorption continue jaune » et enfin la « série 


TT fond conliny Touge celcile 
Derritr. Sona Continu Jeune calculé 


verte ». Dans deux précédents mémoires [1], [2], 
nous avons étudié les absorptions continues rouge 
et jaune et la série Jaune. Nous nous proposons 
d'exposer dans ce mémoire quelques résultats obte- 
nus sur la série verte. 
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F1G. 1. — Courbe d’absorption d’un échantillon de Cu,O de 10 x d’épaisseur à 4,2 °K. 
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F1G. 2. — Courbe d’absorption d’un échantillon de CuO, 
de 10 y: d'épaisseur à 4,2 °K (série verte). 


Evaluation des intensités d’oscillateur. — Nous 
avons montré précédemment que l’absorption con- 
tinue rouge, formée de deux fonds continus, corres- 
pond aux transitions excitoniques indirectes pré- 
vues par la théorie d’Elliott [3] lorsque Ak %Z 0. 
Le carré des coefficients d’absorption respectifs des 
deux fonds continus varie linéairement avec le 
nombre d’onde v. Nous avons done pu prolonger 
ces fonds continus dans tout le spectre et les sous- 
traire de l’absorption totale. De même, la série 
jaune et l'absorption continue jaune correspondent 
aux transitions excitoniques interdites lorsque 
Ak — 0 et k Æ 0, mais petit. La puissance 2/3 du 
coefficient d'absorption du fond continu jaune varie 
linéairement avec le nombre d’onde v. On peut 
donc prolonger ce fond continu jusque vers la série 
verte (fig. 1). La formule admise n’est d’ailleurs 
qu’approximative. 


N°3 


Par ailleurs, ce spectre continu doit se raccorder 
à celui que l’on observe au delà de la série verte, 
du côté des grands nombres d’onde. La théorie ne 
donnant aucune indication quant au raccord de ces 
deux spectres, nous avons raccordé à l’estime selon 
deux tracés extrêmes (1) et (II) Nous avons 
ensuite soustrait cette absorption continue (I) 
et (II) de absorption totale et obtenu ainsi deux 
tracés de l’absorption dans les raies (fig. 2). 

La méthode de Krawetz nous a permis d'évaluer 
les facteurs d'intensité d’oscillateur par cellule élé- 
mentaire des raies de la série verte à partir de 


l'intégrale S 44 K dy de la courbe d’absorption 
de chaque raie étudiée : 


cellule élémentaire 
exp ÉL 


2,45.10—108, 


Nous avons pu évaluer les facteurs f pour les 
tracés (1) et (IT) (tableau 1). On constate que seule 
TABLEAU I 


FACTEURS D'INTENSITÉ D’OSCILLATEUR PAR CELLULE ÉLÉ- 
MENTAIRE DE LA SÉRIE VERTE A 4 OK POUR LES TRACÉS (I) 


ET (11) DU FOND CONTINU D’ABSORPTION. ÉCHANTILLONS 
DE DIFFÉRENTES ÉPAISSEURS. 
= een BP NE 
PT TOP LOS 1 OS fr, 1 65 PT Or 0 
2 3 3,1 3,3 3,4 3,1 3,9 
3 0,4 0,8 0,43 0,81 0,44 0,82 


la raie n — 2 donne, pour les deux tracés des résul- 
tats concordants ; pour les raies de nombre quan- 
tique n > 2, on obtient uniquement un ordre de 
grandeur. À 4 °K, on a ainsi : 


2 — 32.40 


Ân=3 = 0,6.10—5. 


On remarque que le facteur f de la raie n —2 


ne varie pas avec la température, dans la limite 
de la précision donnée par les deux tracés. 


Comparaison des intensités d’oscillateur des raies 
n — 2 des séries jaune et verte. — On peut obser- 
ver que la raie n» —2 de la série verte est plus 


INTENSITÉ D’'OSCILLATEUR DES RAIES 


1 
intense que la raie correspondante de la série jaune. 
En effet : 

HAS 2 0 0e 
fn=2 


m=olfh=s 11. 


= 3,2.10—5 


Ceci peut s'expliquer par la différence entre les 
constantes de Rydberg des séries jaune et verte. 
En effet, la formule théorique du facteur f s’écrit : 


f = Cuime)$ [(n?— 1) fn] 


où est la masse excitonique effective, e la cons- 
tante diélectrique, m la masse de l’électron, n le 
nombre quantique des raies et C'une constante. 

D'autre part, les deux séries peuvent être repré- 
sentées par une formule hydrogénoïde : 


V = Voo — (Rexc/n?) 
et 

Rex = Fal(u/me?) 
où À est la constante de Rydberg de l’hydrogène, 
donc : 

cRexe[Rx — ufme. 


Le facteur f peut alors s'exprimer en fonction de 
la constante de Rydberg : 


1 = C[eRexe [Ra] (n° —1)fn$. 


D’après les travaux de Haken [4], la constante 
diélectrique qui intervient est la constante statique. 
Elle est donc la même pour les deux séries. Le 
rapport des facteurs f est donc égal : 


fn=olhze © 8xef1 bre] < [1 214 /786]5 = 9. 


On obtient ainsi l’ordre de grandeur du rapport 
des f expérimentaux. 


Évaluation du coefficient maximum d’absorp- 
tion et de la largeur des raies. — Nous indiquons, 
dans le tableau I, le coefficient maximum d’absorp- 
tion et la largeur de la raie nr — 2, aux tempéra- 
tures de l’He, l'H, et VPN, liquides, dans le cas du 


TABLEAU II 


FACTEUR D'INTENSITÉ D’OSCILLATEUR PAR CELLULE ÉLÉMENTAIRE, 
COEFFICIENT MAXIMUM D’ABSORPTION ET LARGEUR DE LA RAIE N = 2 DE LA SÉRIE VERTE, 
AUX TEMPÉRATURES DE L'HE, L'H, ET L’N> LIQUIDES, DANS LE CAS DU TRACÉ (II) DU FOND CONTINU 


= _ —— == + = — S 

TROUS T0 Kmax 10cm > ANS CM ANA Citer AVES CT Ô 
[A SL 8,9 74,5 52 91 0,27 
20 3,4 8,5 78 58,5 97,5 0,25 
97 3,4 79 90 72 108 0,20 
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tracé (IT) du fond continu (fig. 3). On remarqueque 


la raie n — 2 est très dissymétrique comme les raies 
de la série jaune. Pour caractériser cette dissy- 
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F16. 3. — Profil de la raie n — 2 de la série verte 
AL 2 0RE (p.108): 


métrie, nous avons mesuré, à partir de l’abecisse du 
maximum de la raie, une demi-largeur vers les 
grandes valeurs de v: Av, et une demi-largeur 
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versles petites valeurs de v : Av. Nous avons pris 
pour paramètre exprimant la dissymétrie, la quan- 
tité : 
ss Avip — A vijo 
Avi + Avi 


Conelusions. — Cette étude nous a donc permis 
de déterminer les valeurs expérimentales de l’inten- 
sité d’oscillateur des raies de la série verte de Cu,0. 
Cette détermination a été faite avec une bonne 
précision pour la raie n — 2 ; pour la raie n —3, 
nous n’avons pu donner qu’un ordre de grandeur. 

Nous avons pu comparer les facteurs f des raies 
n — 2 de la série jaune et de la série verte et expli- 
quer leur différence par la différence des constantes 
de Rydberg des deux séries. 

Nous avons également indiqué le coefficient 
maximum et la demi-largeur, à différentes tempé- 
ratures, de la raie 7% — 2 et montré la dissymétrie 
de cette raie. 


e 
Manuscrit reçu le 17 décembre 1960. 
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EXPOSÉ ET MISE AU POINT BIBLIOGRAPHIQUE 


CAUSALITÉ ET RELATIONS DE KRAMERS-KRONIG (1) 


Par N. G. VAN KAMPEN, 
Institut de Physique théorique de l’Université d’Utrecht. 


traduit par François LURÇAT 
Institut de Physique de la Faculté des Sciences de Lille. 


Résumé. — À l’aide d'exemples empruntés à diverses branches de la physique (élasticité, théorie 
des oscillations, optique, etc.) l’auteur montre que les relations de Kramers-Kronig (relations entre 
les parties réelle et imaginaire de l’indice de réfraction, ou d’une grandeur analogue) sont une 
conséquence du principe de causalité, selon lequel l’effet ne peut se produire avant la cause. Il 
montre ensuite comment ce principe peut être appliqué à la microphysique et donne des relations 
entre les parties réelle et imaginaire des éléments de la matrice S. Il discute la diffusion d’un champ 
électromagnétique classique et celle de particules de Schrôdinger, puis il indique comment ces 
résultats peuvent être étendus au cas de la diffusion de champs quantifiés. 


Abstract — Taking examples in various fields of physics (elasticity, theory of oscillations, 
optics,.…)the author showsthat the Kramers-Kronig relations between real and imaginary parts of 
the index of refraction (or analogous quantity) are a consequence of the principle of causality, i. 
e. the principle that the effect cannot take place before the cause. He then shows how this prin- 
ciple can be applied to microscopic physics, giving relations between the real and imaginary 
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parts of the élements of the S matrix. 


The scattering of (classical) electromagnetic field and of 


Schrôdinger particles are discussed, and the principle of generalization to the scattering of quan- 


tum fields is given. 


1. Introduction. — Quand la lumière se ‘propage 
dans un milieu, il y a dispersion : l'indice de réfrac- 
tion n varie avec la fréquence «© de la lumière. Une 
for mule qui exprime r en fonction de « s’appelle une 
formule de dispersion. En 1872 Sellmeier [1] écrivit la 
formule de dispersion suivante, à partir de l’hypo- 
thèse que l'influence du milieu peut être attribuée à 
différentes sortes de particules élastiquement liées : 


n(ow) —1 = À 
| 5 


Les «, sont les fréquences propres des particules qui 
entrent en résonance. Au voisinage des fréquences de 
résonance, cette formule de dispersion n’est plus une 
bonne approximation. Chaque constante 4, est égale 
au coefficient d'absorption totale dans la raie de réso- 
nance correspondante. Pour les rayons X, les réso- 
nances sont trop serrées pour qu’on puisse les consi- 
dérer séparément ; c’est pourquoi vers 1927 Kramers 
et Kronig [2] proposèrent presque simultanément de 
remplacer la formule de Sellmeier par : 


n(ow) — 1 = Ait lès dr. u) 


2 — w°? 


L'indice de réfraction est-donc déterminé dès que la 
fonction a(w,), c’est-à-dire l’absorption en fonction de 
la fréquence, est connue. Pour que l’intégrale dans (1) 


) Article paru dans Nederlands Tidschrift voor Natuur- 
kunde, 1958, 24, 1-14 et 29-42. 


ait un sens, on convient de prendre pour la singula- 
rité « la valeur principale (?). 

Plus tard 1l apparut que la relation de Kramers et 
Kronig (1) est vraie dans des conditions beaucoup plus 
générales que ce qui précède ne le laisse supposer. Elle 
est entièrement indépendante de tout modèle parti- 
culier pour l'interaction de la lumière avec le milieu, et 
peut s’obtenir d’un point de vue général de la façon 
suivante. Supposons qu'un court train d’ondes lumi- 
neuses tombe sur un milieu. Pour étudier la propa- 
gation à l’intérieur de ce milieu, on doit se représenter 
le train d’ondes comme formé d’ondes monochro- 
matiques (composantes de Fourier) ; chaque onde 
monochromatique individuelle, de fréquence w, se pro- - 
page avec sa propre vitesse c/n(&), et subit un certain 
affaiblissement, caractérisé par un coefficient d’absorp- 
tion a(w). Ceci modifiera notablement la forme du 
train d’ondes ; en général il s’étalera quelque peu. Les 
détails de cette déformation dépendent de la forme 
précise des fonctions n(w) et a(w) (donc de la structure 
du milieu en question), mais il y a une condition à 
laquelle ils doivent nécessairement satisfaire : le front 
ne peut se déplacer avec une vitesse supérieure à c. On 
peut alors montrer de façon purement mathématique 
qu'il est nécessaire et suffisant que la relation (1) soit 


(2) C'est-à-dire que l’on enlève d’abord du domaine 
d'intégration un intervalle (© — €, w + €) symélrique par 
rapport à «, et qu'après intégration on passe à la limite 
EU. 
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vérifiée, pour que cette « condition de causalité » soit 
satisfaite. 

Une condition de causalité analogue existe pour 
toutes sortes de systèmes physiques. Une formation 
générale et peu précise est : l'effet ne peut se produire 
avant la cause. Pour chaque application, cette formu- 
lation doit naturellement être précisée. Ces dernières 
années, ces conditions, et les relations de Kramers- 
Kronig qui s’en déduisent, se sont avérées utiles dans 
des domaines très différents mais surtout pour l’analyse 
de l'interaction entre les nucléons et les mésons x (?). 
Comme dans ce dernier cas on peut à peine recon- 
naître les idées parmi les calculs compliqués, je mon- 
trerai ici sur des exemples simples comment on parvient 
aux relations de Kramers et Kronig. Au paragraphe 2, 
à l’aide d’un exemple tiré de la théorie de l’élasticité, on 
montrera le lien mathématique entre la causalité et 
relations de Kramers et Kronig. Au paragraphe 3 on, 
énumèrera un Certain nombre d'applications et 
d'extensions à d’autres domaines. Le paragraphe 4 est 
consacré à la propagation de la lumière dans un milieu, 
berceau des relations de Kramers et Kronig. Les résul- 
tats seront utilisés au paragraphe 5 pour obtenir des 
renseignements sur la diffusion de la lumière par les 
centres diffuseurs individuels du milieu diffusant. Au 
paragraphe 6 on étudiera cette diffusion individuelle à 
l’aide d’une autre méthode, qui au paragraphe 7 sera 
appliquée à la diffusion des particules par les noyaux. 
Au paragraphe 8 on étudiera, de façon quelque peu 
sommaire, la formulation de la condition de causalité 
en théorie des champs, et son application à la physique 
des mésons. 


2. Exemple tiré de la théorie de l’élasticité. — 
Supposons qu’un barreau, sous l’action d’une force 
appliquée X(f), subisse un allongement x(t). Supposons 
en outre qu’il se produise une hystérésis élastique, 
c’est-à-dire que x(t) dépende non seulement de X({) 
mais encore de «l’histoire antérieure » donc des valeurs 
X(t’) de la force aux instants antérieurs (f < +). Dans 
le cas de petits allongements, cette relation est encore 
linéaire, en sorte qu’elle est de la forme [3] 


= 


On a naturellement 


Gr) = 0 


Gi — 1) X(e) FPE Le Ce are de 


(2) 


TE Ol0 


pour (3) 
C’est dans ce cas particulier la condition de causalité: 
La forme triviale que prend ici la condition de causa- 

lité est due, non pas tant à la simplicité de l’exemple 

choisi, mais au fait que la cause X(?) et l'effet æ(1) sont 
exprimés en fonction du temps. En général cependant 

il sera préférable de décrire les phénomènes à laide 

des fréquences, car la relation entre x et X sera plus 

simple. 
C'est-à-dire qu’au lieu de G{t), X(t), x(t) on intro- 


(%) Dans ce domaine on parle habituellement de « rela- 
tion de dispersion », ce qui crée souvent une confusion 
avec la notion définie ci-dessus de « formule de disper- 
sion». Ici j’utilise toujours les mots « relations de Kramers- 
et Kronig ». 
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duira leurs transformées de Fourier G(w), 54 (w), x (w) 
définies par exemple par : 


(4) 


On obtient alors au lieu de l’ « équation de mouve- 
ment » (2) 


Co 1 Ge) eiot de. 


(5) 


La fonction complexe G(w) décrit les propriétés 
élastiques de la tige, de la manière suivante : si on 
applique à la tige une force périodique de fréquence «, 


il s'établit une oscillation, dont l’amplitude est donnée : 


par |G(w)l et la différence de phase par rapport à la 
force, par l'argument de G(«). Le travail produit par 


Fic. 1. — Tige viscoélastique. 


la force dans l’unité de temps est égal au carré de 
l'amplitude de la force, multiplié par un « coefficient 
d'absorption », 


oIm G(o)]. (6) 
Im désigne la partie imaginaire. La partie réelle 


Rel G (w)] n’a en général pas de signification directe ; 
mais lorsque l’énergie absorbée, et par conséquent la 
différence de phase, sont faibles, on a 


Re G(w)] — |G(w)l 
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donc la partie réelle mesure l’élasticité (à la fré- | 


quence &). 

La transformation de Fourier a sans doute permis 
de mettre l’équation de mouvement sous la forme 
simple (5), mais la condition de causalité ne peut plus 


s’écrire aussi simplement que (3). Comment peut-on … 


exprimer, à l’aide de G(w), que G(t) n’est pas une 
fonction quelconque dans l'intervalle (— co + co 
mais qu’elle doit vérifier la condition (3) ? Il y a là 
une question purement mathématique, à laquelle 
répond le théorème 95 du livre de Titchmarsh sur les 
intégrales de Fourier [4]. Avant de formuler exac- 
tement ce théorème, donnons d’abord une démons- 
tration sommaire : ; 

La première étape de la démonstration consiste à 
remplacer dans (4) la variable réelle «© par un nombre 
complexe w — u + iv. Cela est sûrement légitime 
pour D > 0, car alors Pintégrant est multiplié par le 
facteur exp (—- vi), ce qui ne peut qu’améliorer la con- 
vergence de l'intégrale pour { — + co. Évidemment 
pour { — — co ce facteur serait en général désastreux, 
mais grâce à (3) il ne fera pas de tort. On peut donc 
en posant 


Êle) = [ Gt) cit de () 


Ts 
étendre la définition de G à la moitié supérieure du 


À Da Vrt-orp hpfat 


N°3 


plan complexe, Comme (7) converge uniformément, 
elle définit une fonction analytique sans singularités 
(€ holomorphe »). Donc : par suite de la condition de 


causalité (3), G(«) possède un prolongement holomorphe 
dans la moitié supérieur du plan complexe. En outre, 
on peut voir facilement que (7) tend vers zéro lorsque 
D — co. 

Comme deuxième étape nous choisissons un chemin 
d'intégration fermé, comprenant l'intervalle (— A, 
+ R) sur l’axe réel et un demi-cercle de rayon À? dans 
le demi-plan complexe supérieur. La théorie de Cauchy 
donne pour tout point w situé à l’intérieur du contour : 


“= 4 +R G 
G{w) — e , SL © + 1 40 ; 
2T1J—R © — 2T1J/ (demi-cerle) 


oi On passe à la limite À — co, le second terme 
s’annule, et l’on obtient 


A LT Go) 
GP} Ni —co D. — W (es (8) 


Réciproquement cette formule implique que G(w) 


Fic. 2. — Le plan complexe w. 


est holomorphe dans le demi-plan supérieur et nulle à 
l'infini. 

Comme troisième étape nous mettrons (8) sous une 
forme utilisable pratiquement. Pour cela faisons 
tendre w — u + iv vers une valeur réelle, c’est-à-dire 
v — + 0. L’intégrant acquiert alors une singularité 
pour w — u, et le chemin d'intégration doit passer en 
dessous de ce point. On obtient donc la moitié du résidu 
en ce point, plus l’intégrale sur le reste de l’axe réel 
qui est une valeur principale. Le résultat est : 


Le 


Zs 1 1 co  (G(w) 


.où P signifie que pour la singularité il faut prendre la 
valeur principale. En séparant les parties réelles et 
imaginaires, on obtient 


Rd&u)] = © P [© Iml G(w)] 


do (94) 


—09 OU 


Im Cl] = —1p (© Rd Gla)] 4. 


T 00 0 U 


(9) 


Ces relations ne concernent que les valeurs prises 


par G sur l’axe réel, et relient par conséquent des 
grandeurs physiquement mesurables, D’après notre 
démonstration, elles constituent des conditions néces- 
saires pour qu'il existe un prolongement holomorphe 
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de G(&), qui satisfasse à (8). Réciproquement, on peut 
montrer que chacune des relations (9a), (9b) est aussi 
une condition suflisante. 

Le résultat peut être résumé dans le : 


THÉORÈME PRINCIPAL. — Les quatre énoncés 


. . « . > 
suivants relatifs à une fonction G(«w) sont des formu- 
lations de la même propriété : 


(i) La transformée de Fourier de G (w) vérifie (3). 


(ii) G(w) possède un prolongement holomorphe dans 
la moitié supérieure du plan complexe, qui vérifie (8). 


(ii) Les parties réelle et imaginaire de G(w) sont 
liées par (9a). 
(iv) Les parties réelle et imaginaire sont liées par 


La moitié supérieure du plan complexe de la 
variable w sera par la suite notée 7,. 
La relation (9) est appelée transformation de Hilbert: 


les fonctions Re[ G(w)] et Im] G(w)] sont chacune la 
transformée de Hilbert de l’autre. On appelle aussi (9) 
formules de Plemelj [5]. En physique on les nomme 
relations de Kramers-Kronig. On utilise en particulier 
(9a) sous la forme : 

Rel G(u)] RP. SR do (10) 


(qui résulte de (9a) parce que G(— &) — G(«)). Cette 
relation exprime l’élasticité en fonction du coefficient 
d'absorption (6). Pour les applications pratiques 
c’est (10) qui est utile, mais en théorie, il vaut mieux 
utiliser (9). 


Nous n’avons jusqu'ici — comme c’est l’usage en phy- 
sique — prêté aucune attention aux propriétés mathé- 
matiques que doit posséder la fonction Gt), et en parti- 
culier à son comportement pour { — co. Pour que le théo- 
rème principal soit valable, il faut que G{t) soit de carré 
sommable c’est-à-dire que 


CO k 

[ |G(H)l? de existe et soit finie. (1) 
—C9 

Ceci garantit l’existence des intégrales utilisées, et des 
transformées de Fourier [4]. L’étude des exemples suivants 
(paragraphes 4, 6, 7) montrera que cette condition est 
souvent physiquement acceptable, car elle exprime le fait 
que l’énergie totale doit être finie. Dans l’exemple qui pré- 
cède, cependant, l'intégrale (11) n’a aucune signification 
physique, et la condition n’est donc pas naturelle. Elle 
n’est d’ailleurs pas vérifiée dans le cas le plus simple 
(xt) = cste X ÆX(t)). Sans doute il n’est pas difficile de 
généraliser la démonstration de façon à ce qu’elle com- 
prenne aussi ce cas là ; mais on ne peut pas être sûr de ne 
pas rencontrer de nouvelles exceptions. La seule solution 
satisfaisante est d'adapter la démonstration de telle façon 
qu’on n’impose à G{t) aucune autre condition que celles qui 
résultent de sa signification physique. Mais pour l'exemple 
considéré on n’a pas encore réussi à faire cela. 

Une autre propriété de la fonction G(«) doit encore être 
mentionnée, car elle est indispensable pour un traitement 
mathématique rigoureux. Du fait que w G{w) est holo- 
morphe dans Z+ et que Im[ow G(w)] > 0 sur l’axe réel 
(à cause de sa signification physique de coefficient d’absorp- 
tion), il résulte que Zm[w G(w)] > 0 dans tout +. w G(w) 
appartient donc à la classe des fonctions « quasi-bornées » 
(et même à celle des fonctions « imaginaires positives »), 
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dont on connaît beaucoup de propriétés mathématiques [6]. 
Le fait que l’énergie dissipée soit positive résulte d’autre 
part de la thermodynamique, donc de la mécanique statis- 
tique. Nous rencontrerons d’autres exemples montrant le 
na étroit de la causalité avecles considérations statistiques 
CES 0) . 


3. Autres applications. — Bien que l’exemple traité 
au paragraphe 2 soit très élémentaire il contient les 
traits essentiels que l’on retrouve dans chaque appli- 
cation de la condition de causalité. On peut les for- 
muler ainsi. On donne d’abord une « cause » et un 
« effet », et la relation entre eux écrite sous forme 
générale (cf. (2)). La condition de causalité dit alors 
que cette relation doit être telle que l'effet ne peut 
être influencé par le comportement de la cause aux 
temps ultérieurs. Tant que la cause et l’effet sont écrits 
comme des fonctions de #, la formulation dela condition 
sera la plupart du temps triviale (pas dans tous les 
cas cependant) (cf. (3)). Mais lorsqu'on passe aux trans- 
formées de Fourier, c’est-à-dire lorsqu'on introduit la 
fréquence comme variable indépendante, la condition 
devient moins simple. Elle prend alors la forme d’une 
proposition relative aux propriétés analytiques de la 
«réponse » ou «admittance » pour des fréquences com- 
plexes. En dernière analyse cette proposition peut être 
exprimée comme une relation à laquelle doit satisfaire 
l’admittance pour des fréquences réelles. Cette relation 
exprime la partie réelle de l’admittance comme une 
intégrale portant sur la partie imaginaire ou inver- 
sement : ce sont les relations de Kramers et Kronig (9), 
resp. (10). Pour choisir les grandeurs qui seront consi- 
dérées comme cause et comme effet, et pour écrire la 
relation générale entre ces grandeurs, il suffit de con- 
naître les propriétés physique du système considéré ; 
de même que pour choisir la forme de la condition de 
causalité dans la représentation en fonction du temps. 
Mais la traduction de cette condition dans la repré- 
sentation en fonction des fréquences, et la démons- 
tration des relations de Kramers et Kronig sont de 
pures mathématiques, et sont à leur place dans un 
traité sur les intégrales de Fourier. 

Cette recette est applicable à tout système linéaire 
qui peut être étudié à l’aide de la notion de fréquence. 
Un exemple mécanique simple est l’oscullateur amorti 
auquel on applique la force X(t). L'élongation x vérifie 
l’équation suivante (masse — 1, fréquence propre «, 
amortissement y) : 


©? æ — X(f) (12) 


x + 2yz + 


la méthode habituelle donne la solution : 


PEN Sin art ss 3 
æ(t) SA ee À vtt X(1) de 
— CO Or 


our — VATE — y? = w,. Ceci est la relation (2) 
avec un G(t) qui peut dans ce cas être donné explici- 
tement. On en déduit (1) : 


"+ COSin Gr T , fl 
— OT OC a er 9 
4 j SPAE $ © — (w + 1y)? 4 


(4) On peut aussi trouver ce résultat plus directement 
en prenant dans (12) une force périodique X{(t) = À cos wt 
et en cherchant la solution stationnaire correspondante. 
L'existence de cette méthode plus directe est une des 
raisons pour lesquelles le langage des fréquences est si 


N°3 


Il est manifeste que G(w) est holomorphe dans le 
demi-plan supérieur, et tend vers 0 à l'infini. On peut 
aussi vérifier que les relations (9) sont satisfaites. De 
plus G(w) a deux pôles en dessous de l’axe réel, aux 
points —- «7 — ty. Ces pôles expriment l'existence 
d’une résonance de fréquence w, et de largeur à mi- 
hauteur 2y. La figure 3 représente la variation de 


Re G(w)], Im[ G(w)] et [G(w)| au voisinage de «x. 


F1c. 3. — L’admittance G(o) 
au voisinage de la résonance. 


x x 
x x 


résonances 


uo1}0xD/24 


F1G. 4. — Les pôles de la fonction G{w). 


D'une façon générale, on peut dire que : un pôle de 
l’admittance en dessous de l’axe réel représente une réso- 
nance ; la partie réelle du pôle est (approximativement) la 
fréquence de résonance, la partie imaginaire est la largeur 
de la courbe de résonance. 

Un autre exemple est fourni par la polarisation diélec- 
trique. Ici le champ électrique appliqué joue le rôle de 
la force X ; tandis que la polarisation est l’effet. La 


constante diélectrique e(w) (ou plutôt —) est 


l’admittance G(w). Elle doit donc être holomorphe 
dans le demi-plan complexe supérieur de la variable w, 
et les parties réelle et imaginaire doivent vérifier la 
relation (10) [7]. Un exemple important dans le cas 
présent est donné par la fonction : 


e(@) — 1 = af1 + 1œot)—1. 


Elle représente la relaxation diélectrique avec le 
temps de relaxation +, comme l’a montré Debye [8]. 
D'une façon générale : un pôle sur la partie négative de 


souvent utilisé, en mécanique classique comme en méca- 
nique quantique. Une autre raison est qu’en pratique les 


phénomènes rapidement variables se présentent presque 


toujours sous forme d’oscillations stationnaires. 


No 3 


laxe imaginaire représente un phénomène de relaxation ; 


la distance à l'axe réel est l'inverse du temps de relaxation. 


La fonction citée, appelée fonction de Debye, a à peu 
près la même allure que les courbes de résonance de la 
figure 3, sauf que les parties réelle et imaginaire sont 
inversées. 

L'application de la condition de causalité à la suscep- 
tubulité magnétique () y — y’ + iy” est entièrement 
analogue. La relation de Kramers et Kronig s'écrit 
alors [9] (les relations de Kramers et Kronig contenant 
toujours une valeur principale, nous omettons désor- 
mais le signe P) : 

2 fo w’#”(w’) 
Tr Jo 


LOI do" + x(c0). 


D 02 

On a ajouté la constante y(co) car y n’est pas néces- 
sarement nulle pour les hautes fréquences ©. Il peut 
y avoir des phénomènes de résonance et de rela- 
xation [10]. Pour concrétiser les relations de Kramers 
et Kronig, on peut faire apparaître les courbes repré- 
sentant y’ et y” directement sur un oscillographe [11]. 
D'ailleurs, dans ces applications on n'utilise pas tou- 
jours les relations sous forme explicite, mais on les 
sous-entend dans la mesure ou on n’admet pour 
X' + 14” que des fonctions analytiques de « qui soient 
holomorphes dans Z,. 

On obtient une généralisation des relations de 
Kramers et Kronig en considérant un corps élastique 
auquel sont appliquées n forces X, X,... X,. Soient 
T1 T2 ... Xn les déplacements des points d'application : 


l’admittance sera une matrice n X n, Calc), qui 
vérifie les relations de Kramers et Kronig (2). Les 
relations de Kramers et Kronig pour des fonctions 
matricielles interviennent souvent aussi dans d’autres 
cas ; en outre, on peut tenir compte des effets ther- 
miques en ajoutant des variables X, (température) 
et x, (qualité de chaleur fournie) [12]. 

Un cas important est celui des réseaux électriques [13]. 
La tension appliquée à un « dipôle » est considérée 
comme la cause, le courant comme l'effet. L’admit- 
tance 1/Z(w) possède donc les propriétés qui ont été 


démontrées plus haut pour G(w). On peut démontrer 
des propriétés analogues pour les quadripôles et les 
réseaux plus compliqués (pour lesquels l’admittance 
est une matrice). Ces propriétés analytiques de la ma- 
trice admittance ont d’ailleurs été déjà démontrées 
d’une autre façon, en particulier par Foster [14] (?). 


(!) Pour être en accord avec la mécanique quantique 
nous avons pris pour facteur de temps exp (—1ot) et non 
exp (zot). C’est pourquoi nos définitions des qualités com- 
plexes set y s’écartent des définitions usuelles e — e°— 1e, 
resp x = x — ix’. Il en résulte aussi que les moitiés supé- 
rieure et inférieure du plan complexe sont interverties ; 
mais les relations de Kramers et Kronig ne sont pas modi- 
fiées (puisqu'elles ne contiennent que des grandeurs réelles). 

() Et qui de plus est symétrique par suite des relations 
d’Onsager. 

() La démonstration de celui-ci repose sur le fait que 
tout réseau peut être décomposé en éléments d’ont l’impé- 
dance est de la forme À — iwL — 1 fiwC. La propriété de 
la partie imaginaire mentionnée à la fin du précédent para- 
graphe joue ici un rôle indispensable. On admet généra- 
lement aussi la possibilité d’éléments sans résistance 
ohmique (ou plutôt de résistance négligeable). Dans ce cas, 
certains pôles, qui autrement se trouveraient en dessous de 
l’axe réel, peuvent être situés sur cet axe. 
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La raison pour laquelle on étudie ces propriétés géné- 
rales de l’admittance en théorie des réseaux est qu’elles 
permettent de répondre à la question : quelles sont les 
fonctions admittances, données à priort, qui peuvent 
être réalisées à l’aide d’éléments convenablement 
choisis ? 

Considérons maintenant un exemple élémentaire : la 
diffraction de la lumière par un réseau. Supposons que 
tombe sur le réseau une courte impulsion lumineuse, 
c’est-à-dire une superposition d’ondes planes avec des 
amplitudes et des phases telles qu’elles forment un 
« paquet d’ondes » limité dans l’espace, par exemple 
une fonction à. Le spectroscopiste qui observe sous 
un angle Ü reçoit une seule composante de Fourier, 
donc une onde sinusoïdale, qui est par définition infinie. 
Ceci est manifestement en contradiction flagrante avec 
la causalité. 

Pour résoudre ce paradoxe, il faut remarquer que la 
lumière reçue dans la direction 0 n’est pas strictement 
monochromatique, mais occupe une bande de fré- 
quences de largeur Av — c(B sin 0)-1 où B est la 
largeur du réseau [15]. Le paquet d’ondes émis dans 
cette direction n’a donc pas de raison d’être infiniment 
long ; on peut voir facilement que sa longueur n’est 
pas inférieure à B sin 0. Mais un tel élargissement d’un 
paquet d’ondes incident en fonction à n’a rien de 


Fi. 5. — L’étalement d’un paquet d’ondes 
à la traversée d’un réseau. 


paradoxal ; sa longueur est précisément la différence 
des chemins parcourus par la lumière qui atteint 
l'œil via le trait supérieur du réseau, et par celle qui 
est passée par le trait inférieur. D'une façon générale 
on peut voir ainsi que dans un interféromètre de 
pouvoir de résolution Av il doit y avoir une différence 
de marche d’au moins c/Av. On peut comparer ce 
résultat à la relation d'incertitude de la mécanique 
quantique. 


4. Propagation de la lumière dans un milieu. — 
Considérons une plaque formée d’un milieu diélec- 
trique d’épaisseur D. Un paquet d'ondes électroma- 
gnétiques arrive de la gauche et tombe sur la plaque 
sous incidence normale. Chaque composante de Fourier 
subit en traversant le milieu un affaiblissement dû à 
l'absorption, et un déphasage correspondant au chemin 
optique parcouru dans le milieu. Par suite, l'amplitude 
de He onde monochromatique sera, après son 
passage dans la matière, multipliée par les facteurs 
exp (— xD/2) et exp (iwnD|c), où a(w) est le coef- 
ficient d'absorption et n(w) l'indice de réfraction. Ces 
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grandeurs peuvent être groupées en un nombre d'ondes 
complexe 

x(&) = wn(a) fe +3 ia(o) = 
où N(@) — n(w) + 1(c/20) a«(w) est l'indice de réfrac- 
tion complexe usuel. Si on écrit pour le champ élec- 
trique incident sur la face gauche de la plaque 


@N(o) [6 


Eialt) = f£nto) etot do 
on a pour le champ sur la face droite 


Ésut(o) = étx(o)D En(w). (14) 


a ——— 


D 


Fi. 6. — Déformation d’un paquet d’ondes 
au passage dans un milieu. 


La condition de causalité exige alors que exp (1x(«)D) 
possède un prolongement analytique dans Z,. Ceci 
conduit à la relation de Kramers et Kronig mentionnée 
au paragraphe 1 : 


c fS «lo 
n'!o) 1 L do”. 
\ 2 2 
Tr Jo &?— 0 


(15) 


La démonstration est la suivante [16]. On doit donc 
exiger que, Si Ein't) = 0 pour t < 0, Évutlt) = 0 pour 
t < DJC. Donc: si Æ£inlw) est holomorphe dans +, 


Écut'o) exp (— 16D]C) 


doit l’être également. Donc exp | {i4(o) — io fe)D } doit 
être holomorphe dans Z+. Il n’en résulte pas encore que 
(wo) — wfc doive être holomorphe, car les zéros de 


exp{ ix{w) — iwfc)D } 
devraient conduire à des singularités logarithmique de x{w). 
Mais ceci est exclu, car D est arbitraire. De plus X(&) — w/c 
elle même n’est pas nulle à l'infini ; il faut diviser par © 
pour obtenir uae fonction nulle à l'infini. Pour cette der- 


nière fonction, les formules de Plemelj sont valables, d’où 
résulte l’équation (15). Il faut encore remarquer que 


Einlo)|? do représente l'énergie totale du paquet 


d’ondes incident, que l’on peut supposer finie ; la théorie 
de la transformation de Hilbert peut donc s’appliquer ici 
de façon tout à fait rigoureuse. 


La relation obtenue est satisfaite en particulier par 
les formules de dispersion connues. Lorentz obtenait 
par sa théorie des électrons [17] 


2e? 
3m 


e2]m 
NO 1? le re 
Tr 


OË — ©— 1yYeË ©’ (16) 


(les &, sont les différentes fréquences de résonance ; 
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les NV, sont les nombres d'électrons correspondants). Le 


premier membre contient l'indice de réfraction com- 


plexe N(w), en sorte que (16) exprime aussi l'indice de 
réfraction ordinaire r(«) et le coefficient d’absorp- 
tion œ(w) en fonction de grandeurs relatives au milieu. 
On peut vérifier que ces expressions de 72(&) et «(w) 
vérifient la relation (15), mais il est plus simple de 
remarquer que (16) est holomorphe dans 74 et nulle à 
l'infini, done que l’on doit bien avoir (15). Il en résulte 
que dans un tel milieu aucun paquet d’ondes ne peut 
se propager à une vitesse supérieure à c. Ceci a été 
vérifié explicitement autrefois par Sommerfeld et 
Brillouin [18] pour résoudre le paradoxe suivant. Au 
voisinage de &-, l'indice de réfraction n est une fonc- 
tion décroissante de « (dispersion anomale), et peut 
même être inférieur à 1 ; la vitesse de propagation c/n 
devrait être supérieure à c, en contradiction avec la 
théorie de la relativité. Un traitement approximatif du 
problème montre que la plus grande partie du paquet 
d’ondes se propage en eflet avec la vitesse c/n ; mais 


précisément cette approximation n’est plus valable. 


dans le domaine ou N(w) est rapidement variable, 
En même temps il apparut que le front d’un paquet 
d’ondes, dans un milieu pour lequel est vérifiée la 
relation (16), ne peut jamais se déplacer à une vitesse 
supérieure à c, ce qui est évident pour nous aujour- 
d’hui. 


Une particularité intéressante est que la théorie de 
Lorentz ne donne pas en fait un terme d’amortissement 
— iYor &,,mais — ya. Pour toutes les valeurs raison- 
nables de «w, cela revient au même, mais en toute rigueur 
il en résulte que chaque dénominateur dans (16) possède un 
troisième zéro, voisin de w — if, donc dans 7+. Par 


suite la propagation de la lumière dans un milieu d’élec- : 


trons de Lorentz n’est pas strictement causale, mais il 


apparaît un « précurseur» qui est en avance d’une longueur. 


égale au rayon de l’électron. Heureusement pour Sommer- 
feld et Brillouin, ils utilisaient la formule approchée (16) ! 

Un tel pôle acausal apparaît aussi pour un électron de 
Lorentz libre ; il est lié au fait que l’équation de mouve- 
ment de Lorentz possède une solution anormale pour 
laquelle (indépendamment d’une éventuelle force appliquée) 
la vitesse de l’électron croît exponentiellement comme 
exp (t/y) (solution « auto-accélérée »). Lorentz ne s’inquié- 
tait pas de cela parce qu’il considérait le terme d’amortisse- 
ment simplement comme une approximation, valable seule- 
ment pour |w| << 1/y. Cependant Dirac [20] montra que 
les mêmes difficultés se produisent dans la théorie relati- 
viste rigoureuse des particules chargées, et qu’elles sont 
fondamentalement liées à l’électrodynamique de Maxwell- 
Lorentz. En ce qui concerne la mécanique quantique, on 
crut longtemps qu’ilne peut y apparaître de telles solutions. 
Celà parce que l’effet de l'interaction entre la particule et le 
champ ne peut être calculé qu’en faisant un développement 
en série de puissances de la charge e; il est clair qu’ainsi on 
ne peut trouver une solution où y figure au dénomina- 
teur [21]. Cependant dans un modèle simplifié qui peut être 
traité exactement, il apparut qu'il existe un analogue quan- 
tique de la solution anomale [22]. Pauli trouva alors une 
anomalie analogue dans le modèle simplifié de Lee pour la 
théorie mésonique [23] ; on soupçonne que cette anomalie 
est liée aux difficultés fondamentales de la théorie des 
champs et c’est pourquoi les « solutions fantômes » corres- 
pondantes sont venues au centre de l’intérêt. Il est remar- 
quable que Delsarte avait trouvé des solutions anormales 
de ce genre dans un cas très élémentaire, celui du calcul du 
champ électromagnétique dans deux milieux contigus qui 
possèdent des constantes diélectriques et des perméabilités 
différentes, mais indépendantes de « [24]. 


“ 


La formule de dispersion de Sellmeier peut s’obtenir 
à partir de (16) comme approximation, en négligeant 
le terme d'amortissement au dénominateur, Mais alors 
les pôles qui se trouvaient dans la moitié inférieure du 
plan complexe viennent sur l’axe réel, et par suite la 
formule de Sellmeier ne vérifie pas les relations de 
Kramers et Kronig (15) (bien qu’elle ait amené à la 
découverte de ces relations (1))}. La formule de dis- 
persion quantique de Kramers-Heisenberg ne vérifie 
pas non plus (15), parce qu’elle est également en 
défaut dans les régions de résonance. Mais si on intro- 
duit — avec Weisskopf et Wigner — des termes 
d'amortissement convenables dans les dénominateurs 
de résonance, alors (15) est satisfaite [25, 26]. 


Le mouvement des électrons dans un cristal peut être 
traité d’une façon analogue. Si une onde électronique 
« monochromatique » d'énergie Æ tombe normalement sur 
le cristal, elle se propage à l’intérieur du cristal comme 
une onde de Bloch wx{x) exp ? 1kx —1E{(k) t?, où k estun 
paramètre et wx(x) une fonction périodique de même 
période que le réseau. La fonction ÆE(k), dont la forme 
exacte est déterminée par-les propriétés du cristal, à 
l'allure générale indiquée par la figure 7 [27]. La relation 


Fic. 7. — La formule de dispersion pour les électrons 
dans un réseau. 


entre l’onde incidente et l’onde qui sort par l’autre face du 
milieu sera donnée, exactement comme dans (14), par le 
facteur exp (1k(E) D). La causalité laisse prévoir que (FE) 


_ possède un prolongement holomorphe dans la moitié supé- 


rieure du plan £. À première vue cela semble impossible : la 


_ fonction k(E) est formée de branches séparées par des 
_ bandes interdites telles que (Æ1, E,). Mais une étude appro- 


fondie montre qu’on a bien le résultat voulu ; dans les 
bandes interdites k est imaginaire pur, et ces valeurs imagi- 
naires, avec les valeurs réelles indiquées sur la figure 7 
forment une fonction que l’on peut prolonger. Cette fonction 
est holomorphe dans Z+ et possède des points de bran- 
chement en Æ1, FE,, etc. Geci peut être démontré mathé- 


-matiquement à partir des propriétés de l’équation de Schrô- 


dinger pour les électrons dans un potentiel périodique. Dans 
ce cas la condition de causalité ne nous apprend rien de 
nouveau, Car nous disposons grâce à l'équation de Schrô- 
dinger d’une description beaucoup plus détaillée. 


5. Processus individuels de diffusion. —— On sait que 
l'indice de réfraction n et le coefficient d'absorption & 
d’un milieu s'expliquent par la diffusion multiple de la 
lumière par un grand nombre de centres diffuseurs 


() En même temps que Kramers et Kronig, Kallman 
et Marc (Ann. Phys., 1927, 82, 585) écrivirent aussi une 
formule de dispersion pour le cas ou les fréquences de 
résonance sont réparties de façon continue. Ils partaient 
également de (16) mais ne négligeaient pas le terme d’amor- 
sement au dénominateur. C’est pourquoi ils ne trouvaient 
pas les relations de Kramers et Kronig. 
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individuels contenus dans ce milieu. C’est pourquoi la 
relation (15) peut aussi être interprétée comme une 
propriété de ces processus de diffusion individuelle. 
C’est ce point de vue que nous allons développer. La 
diffusion par un centre diffuseur donné (que nous 
appellerons dans la suite « noyau ») peut être décrite 
comme suit. L’onde incidente (dans la direction z) sera 
représentée, en mettant à part le facteur de temps 
eXp (— iot), par exp (1#z), où 4 — w/c. L’onde diffusée 
est, du moins à grande distance, une onde sphérique 
sortante (1/r) exp (ikr), avec une amplitude f qui 
dépend de la direction. L’onde est donc asympto- 
tiquement 


eirr 
ds + f(o, 0). (17) 

L’ « amplitude de diffusion » f(w, 0) ainsi définie 
est une fonction de k, et de l’angle 0 avec la direction 
de l’onde incidente. Le module de f est une mesure (1) 
de l’intensité de l’onde diffusée, et l’argument de f 
est le déphasage dû à la diffusion [28]. 

La relation entre l’amplitude de diffusion f et les 
propriétés de l’onde diffusée n’est toutefois valable 
que pour 0 +0; dans la direction en avant, 0 — 0, 
l'interprétation physique de la grandeur complexe f est 
radicalement différente, En effet, dans cette direction 
l’interférence entre les deux termes de (17) devient 
importante, en sorte qu’une séparation entre ondes 
incidente et diffusée n’a plus de sens. On constate que 
précisément l’amplitude de diffusion en avant f(k, 0) 
est liée aux grandeurs n(&) et «(w) qui expriment les 
propriétés macroscopiques du milieu ; on a 


DEN ne 0 ne die N TE Refflk, 0)] 


E (18) 


N est le nombre de « noyaux » par unité de volume du 
milieu. 


Pour donner une brève démonstration de (18), consi- 
dérors la situation représentée par la figure 8 : une onde 


F1G. 8. — Démonstration du théorème optique. 


plane monochromatique exp (ikz) tombe normalement sur 
une couche de matière d'épaisseur très faible D. L'état de 
l’onde en un point situé à la distance z après la couche est 
donné par : 


gite ei(N—1KD — cikz { 4 + i(n — 1) XD k. 
D’autre part il résulte de la superposition de l’onde inci- 
dente et de toutes les ondes diffusées individuelles ; comme 


(:) En effet, il est lié à la section efficace différentielle 
par la relation do /dQ = |f}?. 
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il y a ND centres diffuseurs par unité de surface, cela 
donne 


eikz + jé r—1 eikrf(X, 0) ND d0 
= gr + 27 ND fe eikr f[k, arc cos (z/r)] dr. 
2 
Comme kz > 1, le facteur f[ ] varie beaucoup plus 


lentement que exp (1kr) donc on peut poser (méthode dela 
phase stationnaire). 


. CA . 
eikz L 92% ND(K, 0) dé ir 
A 


= ekrf4 + 2x NDf(k,0) (— ik) 1. 
En comparant les deux expressions on trouve : 
N—1—=27rk 2? Nf{(k, 0) 


(19) 


ce qui est identique à (18). 

Bien que les relations (18) soient fréquemment utilisées 
en liaison avec la condition de causalité, il ne faut pas 
perdre de vue qu’elles sont d’un caractère tout différent, 
et n’ont en principe aucun rapport avec cette condition. 
Cela est apparent si on considère qu’elles ne concernent 
qu'une seule valeur de «. La première des deux for- 
mules (18) est appelée théorème optique, ou aussi formule 
de Bohr-Peierls-Placzek [29]. Remarquons en passant que 
par suite de (19) [f(k, 0)|/k est bornée ; la limite supérieure 
est (27 ND)-1 où ND (nombre de centres diffuseurs 
par cm?) doit être choisi de façon à ce que les centres 
diffuseurs soient sans interaction mutuelle. 


La relation (18) nous permet maintenant de for- 
muler la relation (15) comme une propriété du pro- 
cessus de diffusion individuel. On trouve que f(w, Ü) 
possède un prolongement holomorphe dans 7} et que 
sur l’axe réel la partie réelle est liée à la partie imagi- 
naire par 


2 ©? 


T 0 


co ]m{f(w’, 0)] do’ 


2 Ta 


Re[f(k, 0)] = 


D ) (20) 

On peut aussi démontrer la relation de Kramers et 
Kronig (20) sans faire le détour consitant à passer par 
le milieu, mais en considérant dès le début la diffusion 
par un noyau. La condition de causalité implique alors 
qu’en un point situé à une grande distance au delà du 
noyau, il ne peut arriver aucune onde diffusée avant 
l’onde incidente. Pour démontrer des relations de 
Kramers et Kronig analogues à (20), pour les directions 
autres que la direction en avant, 0+ 0, la même condi- 
tion nécessaire s’appliquera aux processus de diffusion 
individuels, bien qu’elle n’ait plus alors de relation 
avec les propriétés macroscopiques du milieu. Mais il 
apparaît en même temps que ce n’était pas un hasard 
si la direction de propagation vers l’avant occupait une 
place à part. 

Pour voir cela, considérons un paquet d’ondes planes 
qui se propagent dans la direction z et sont diffusées 
par un noyau (fig. 9). Si ce paquet a un front raide 
qui se propage suivant z — ct, un observateur placé 
à la distance r à droite du noyau ne recevra sûrement 
aucune onde avant l'instant { — r/c. Cependant, pour 
un.observateur dans une autre direction, 0 + O, 
Pinstant où le front de l’onde diffusée l’atteindra 
dépend de la dimension du noyau. La condition de 
causalité perd ainsi son caractère fondamental, et ge 
dégrade en une propriété de la diffusion par un centre 
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diffuseur de dimension finie. Néanmoins elle est utile 


aussi sous cette forme, par exemple pour la théorie de 


la diffusion par les noyaux atomiques, parce que dans … 


N3 


ce cas la condition d’une dimension finie est remplie à . 


une très grande approximation. 


Supposons qu’on peut considérer le noyau comme . 
une sphère de rayon a ; de sorte que pour r > a les 


équations du champ libre sont vérifiées (fig. 9). Un 


2a sin ARE 


Fic. 9. — Condition de causalité 
et dimensions du centre diffuseur. 


paquet d’ondes planes tombe sur le noyau ; l’onde « 


totale s'écrit donc 
(e,e] 
i, A(o) eio(z/o—0 + 71 ciotr/o—0 fo, 6) } do. 
—O0 


Choisissons pour A(w) une fonction qui possède les 
propriétés cités dans le théorème principal ; le paquet 
incident d’ondes planes est alors nul pour z/c —t > 0, 
il a donc un front raide qui se propage selon z = ct. 
Un observateur placé à une grande distance r, dans 
la direction 0 ne peut recevoir aucune onde diffusée 
tant que : 


t € (r — 2a sin 0/2)Je. (21) 


Donc la quantité 


rer Lee A(o) flo, 8) eietr—0 do 
—co 


doit être nulle si £ vérifie (21), et par suite 
A(o) f(w) exp [2i(w/c) sin 0/2] 


doit posséder les propriétés du théorème principal. Il 
suffit pour cela que f(w, 0) soit holomorphe dans la 
moitié supérieure du plan complexe « et que f(w, 0) 
exp [2iw(a/c) sin 0/2] soit bornée dans cette région. 
L’adjonction du facteur exponentiel fait que cette 
condition est nettement plus faible que celle imposée 
à f(w, 0), parce que le facteur exponentiel déeroît lui- 
même rapidement lorsque la partie imaginaire de « 
croit. La relation de Kramers et Kronig est main- 
tenant : 


Ref, 0) e#elt0)] 


. ES D do 
1 J 00 (0/2 — vw?) 


ou {(0) — (afc) sin (0/2). La même formule (22) est 
valable aussi pour des noyaux de forme quelconque; 
si on choisit convenablement /(0). (22) est moins élé- 
gante et moins utilisable que (20), mais on peut vérifier 


, Imi(e, 0) etiet(0)] | (22) 


f 


sg 4 


| 
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sur l’exemple de la diffraction de la lumière par un 
réseau ($ 5) qu’il n’est pas possible de supprimer le 
facteur exponentiel. 


Dans ce qui précède on a supposé implicitement que la 
particule diffusante est fixée dans l’espace, ce qui est une 
bonne approximation pour un noyau atomique. Si on tient 
compte du recul du noyau on trouve, d’après un calcul 
connu pour l’effet Compton, que la lumière diffusée dans la 
direction Ô possède une fréquence ©; inférieure à ©, et 
donnée par : 

Pr ls () 


O7 o Mc? (23) 


ou M est la masse du noyau. L’onde diffusée est donc 
1 [ A(o) f{, 9) cixtre—b de 


et c’est elle qui doit être nulle si (21) est vérifié. Donc, 
A(o) f(w,0) e2ie1(@ dodo: 


possède en tant que fonction de ‘w les propriétés du théo- 
rème principal. Puisque le demi-plan supérieur © est 
appliqué par (23) de façon biunivoque sur le demi-plan 
supérieur &w, nous pouvons conclure à nouveau que f(w, 0) 
est holomorphe dans Z+. Avant de pouvoir tirer de là 
une relation de Kramers et Kronig il faut encore examiner 
l’allure de la fonction pour © — co et w1 — co. On trouve 


alors, par exemple pour un diffuseur ponctuel (/(0) = 0), 
1 6 
2 Raf(o, 0)] = 2 [25 0 PO qu 


ou 1/v est une abréviation pour le second terme de (23). 


6. Diffuseur de rayon fini à symétrie sphérique. — Si 
le noyau diffuseur possède une structure à symétrie 
sphérique, il est avantageux de décomposer l’onde 
incidente en ondes multipolaires, parce que dans ce 
cas chaque onde multipolaire est diffusée séparément, 
indépendamment de la présence des autres. Pour sim- 
plifier nous ne considérerons que la première, l’onde 
monopolaire ou onde s. On peut se la représenter 
comme une onde sphérique, qui tombe sur le noyau en 
provenance de toutes les directions, le traverse, puis se 
propage comme une onde sphérique jusqu’à l'infini. 
L'influence du noyau se traduit par un déphasage que 
l'onde subit en traversant le centre. S'il n’y a pas 
d'absorption, l’amplitude de l’onde sortante est égale 
à celle de l’onde entrante, 

Pour simplifier les formules, nous admettrons que 
le champ de l’onde est décrit par une fonction com- 
plexe ufr, t) qui satisfait en dehors du noyau (r > à) à 
l'équation d’onde 

Au — ü = 0 (24) 

Ceci est valable en particulier pour le champ électro- 
magnétique, si on fait abstraction de son caractère 
vectoriel et qu’on pose c — 1. En l’absence de diffu- 
seur, l’onde s sera représentée par 


1 ; ; 
er 12 oi — — —t 
Bt — "7 sin or ete — = CR ARO RERO (TE 4 
(25) 


La fréquence « est arbitraire, et doit parcourir 
toutes les valeurs de — oo à + co pour que les solu- 
tions (25) forment un système complet. Contrairement 
au cas de la description par ondes planes, le nombre 
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d’ondes est ici un scalaire qui prend des valeurs néga- 
tives ou positives. Le premier terme de (25) repré- 
sente l’onde sphérique qui arrive avec la vitesse 1 ; le 
second terme est l’onde sortante. Si maintenant on 
introduit un noyau diffuseur, l’onde sortante est mul- 
tipliée par un facteur de phase e2int@) — S{w), en 
sorte que le champ de l’onde est maintenant : 
ufr, t) — Le—ior — S{w) er } e—iot, (26) 
La relation avec la section efficace est donnée par 
une formule bien connue [28] 


oo) = 4rw—°? sin? n = 2rw—°(1 — Re [S]). (27) 


Un état de diffusion non stationnaire quelconque, à 
symétrie sphérique, peut toujours s’écrire comme une 
superposition de solutions stationnaires (26) : 


ru(r, t) = [7 Alto) { e—ier — S(o) ior }e—iot do (28) 


A(o) est une fonction arbitraire, qui décrit la forme 
du paquet d’ondes incident. Choisissons encore pour 
A(o) une fonction qui vérifie les conditions du théo- 
rème principal. Alors le paquet d’ondes incident est 
nul pour r + t < 0, donc pour t < —r. La condition 
de causalité exige alors que l’onde soit nulle pour 
t < r — 2a. Donc 


A(o) S(ow) e*iwa 


doit également satisfaire aux conditions du théorème 
principal. De là on peut déduire rigoureusement [30] 
que S est holomorphe dans la moitié supérieure du 
plan complexe w, et que dans cette région on a : 


LS(«) e’ioal < 1. (29) 


De là résulte la relation de Kramers et Kronig 


20 fo Re[S(w'’)eo'a] , 
Im{S(w) eva] — € f eee do”. 
En particulier on peut en déduire pour un diffuseur 
ponctuel (a — 0) à l’aide de (27) : 
2 9 /2 / 
SU tee on 


2T 


Per à do, 
TT GC) — D 

En d’autres termes : dès que o, donc la partie réelle 
de $, est connue pour toutes les valeurs de «w, on en 
déduit S elle même, donc aussi la partie imaginaire. 
Ilest vrai que dans le cas présent ceci ne donne rien de 
très nouveau, car comme [S|? = 1, la partie imaginaire 
est déterminée (au signe près) par la partie réelle ; 
mais 1l est souvent important que le signe puisse être 
aussi déterminé. 

Par exemple dans le cas d’une diffusion faible, ce 
signe permet de déterminer si la force qui produit la 
diffusion est attractive ou répulsive. Au paragraphe 8 
nous verrons un autre exemple emprunté à la physique 
des mésons,. 


La détermination du signe n’exige d’ailleurs pas toujours 
la relation de Kramers et Kronig complète. En effet, 
de (29) résulte l’inégalité [30], [31] 

dnfdw > — a. 

Il suit de là que n ne peut pas décroître trop vite lorsque 

la fréquence croît. Si on fait une estimation raisonnable de a, 


188 


et qu’on puisse trouver un domaine de fréquences où |n]| 
varie fortement, cette inégalité peut être utilisée pour 
exclure une des possibilités n = + |n|. 

Pour la diffusion des particules, on a l'inégalité 
dnfdk > — a, à condition qu’il n’y ait pas d’états liés. 
Dans le cas ou il y a des états liés, on peut encore écrire une 
inégalité du même genre, mais dont le second membre 
contient des termes relatifs aux états liés. 


7. Diffusion des particules. — Jusqu'ici nous avons 
considéré un champ qui obéissait à l'équation 
d’ondes (24),et dont le nombre d’ondes était par con- 
séquent lié à la fréquence par la formule de dispersion 
simple + © — k. Pour des particules obéissant à l’équa- 
ton de de Schrôdinger non relativiste, on a au contraire 
& — hk?/2m. Le paquet d’ondes (28) doit maintenant 
être remplacé par 


VO 

rqlr, t) ik AUR) Êe—itr — 5 (R) ettr } e—iotmt dk (30) 
ou À et «(t) prennent seulement des valeurs positives. 
Il en résulte qu’on ne peut plus choisir A(k) de façon 
à ce que le paquet d’ondes incident ait un front raide ; 
il y à inévitablement un précurseur infiniment long 
bien que faible. Par suite il y a toujours une petite 
chance de trouver la particule loin en avant. La condi- 
tion de causalité ne peut donc plus être formulée de 
la même façon que jusqu'ici. Il est raisonnable de 
poser maintenant la condition suivante : à chaque 
instant t, la probabilité de trouver une particule sortante 
n'est pas supérieure à la probabilité que la particule 
. incidente ait atteint le noyau. Cela revient d’ailleurs 
simplement à ceci que la probabilité de trouver une 
particule en dehors du noyau ne peut jamais être supé- 
rieure à 1. 

En partant de cette condition de causalité, on peut 
démontrer de façon entièrement rigoureuse [32] que 
D(w) est holomorphe dans 7}, mais peut avoir des 
singularités pour « réelle et négative. — En fait on 
savait déjà, dans le cas où le noyau consiste en un 
puits de potentiel, que les états liés donnent lieu à des 
pôles sur l’axe réel négatif. Cependant l'explication 
que fournit la condition de causalité est d’une part plus 
générale, car elle ne suppose rien sur le mécanisme de 
la diffusion ; mais d’autre part elle est moins détaillée, 
parce que l'existence de singularités autres que les 
pôles ne peut être exclue. 

Grâce au fait que S(w) est holomorphe dans 7}, il 
est encore possible de démontrer des relations entre 
Re[S(w)] et /mS(w)]. Mais malheureusement il inter- 
vient dans ces relations des intégrales sur l’axe réel 
tout entier (cf. (9)), alors que seules correspondent à 
des états diffusion physiquement réalisables les valeurs 
positives de ©. Pour les valeurs négatives de «, la 
fonction S(w) est définie par prolongement analytique. 
Pour utiliser ces relations dans la pratique, on doit 
obtenir des renseignements d’une autre façon sur 
l'allure de ce prolongement analytique sur l’axe négatif. 
Dans tous les cas où S(w) a été calculée explicitement, 
il apparaît que S(w) est réelle pour © < 0. On peut 
donc exprimer Re[S(«)] comme une intégrale portant 
sur Im[S(w)] pour les seules valeurs positives de «, 
plus une somme de demi-résidus provenant des pôles 
sur l’axe « négatif. Cependant ce résultat ne peut 
être déduit de la seule condition de causalité. 

Pour les particules relativistes également il y a une 
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«région non physique » analogue, celle comprise entre 
— mc? et me?. Comme ces valeurs de l’énergie n’inter- 
viennent pas dans les états stationnaires de diffusion, 
elles créent premièrement une difficulté pour formuler 
la condition de causalité, et deuxièment font apparaître 
un élément indéterminé dans l’application des formules 
de Plemelj. Cependant là encore on peut soupçonner 
que S(w) est réelle dans ce domaine, en sorte qu’on a 
seulement à tenir compte des résidus correspondants 
aux états liés. 

Jusqu'ici nous n’avons considéré que la diffusion 
pure, c’est-à-dire les cas où le noyau reste après la 
diffusion dans son état initial. On peut aussi étudier 


6 émet 


les réactions, où le noyau réémet une particule diffé- 


rente de la particule incidente (ou la même, mais 
avec une direction de spin différente). La matrice S 
est alors une véritable matrice S,8, où «& et B désignent 
les «voies », c’est-à-dire les différentes possibilités pour 
l’état initial et final. On n’a pas encore traité ce pro- 
blème de façon générale. Comme résultat partiel, on 
peut cependant mentionner que la formule connue de 
Breit-Wigner, qui donne la section efficace 6,8 au voisi- 
nage de la résonance, résulte directement de la condi- 
tion de causalité. 


À ce point nous devons dire quelques mots sur d’autres 
méthodes pour trouver les propriétés générales de la ma- 
trice S. Nous n’aborderons pas les démonstrations qui uti- 
lisent des modèles détaillés pour le noyau [33]. 

Heisenberg [34] tirait des conclusions sur les propriétés 
analytiques de S(«) en partant de l’hypothèse que les 
états stationnaires de diffusion (30) doivent former pour 


t — 0 un système complet de fonctions de r à l'extérieur | 


du noyau. Mais cette démonstration est inexacte ; on peut 
montrer que pour que le système soit complet il suffit 
que S(w) possède des propriétés beaucoup plus faibles [35]. 

Wigner et Eisenbud [36] supposent que l'interaction à 
l’intérieur du noyau est représentée par un hamiltonien 
hermitique, qui possède un système complet de fonctions 
propres à l’intérieur du noyau. Ils considèrent non seule- 
ment la diffusion, mais aussi les réactions. Au lieu de la 
matrice #, ils préfèrent utiliser la « matrice dérivée » R, 
qui est liée à S$ par la relation 


etika S — (1 + ikR)[(1 — ikR). 


Ils trouvent que R(«) est une fonction matricielle de la 
forme : 


(31) 


Les wr sont les fréquences de résonance, et les yŸ déter- 
minent les probabilités des différentes voies. Si on retraduit 
les propriétés analytiques de (31) en propriétés de la matrice 
S', on trouve les mêmes que celles qui résultent de la condi- 
tion de causalité, et en plus que S est réelle pour © < 0 
[37]. 

L’inconvénient de cette approche est que la relation 
de &r et PL avec les fréquences de résonance et probabilités 
de transition observées n’est pas stricte, mais dépend 
encore du choix arbitraire de a. Si on considère, dans l’es- 
pace vide, une petite sphère mathématique de rayon a 


comme un centre diffuseur, on trouvera une matrice À 


de la forme (31) bien qu’il ne soit pas question ici de réso- 
nance. 

Schutzer et Tiomno [38] ont essayé de déduire ces mêmes 
propriétés de la condition de causalité, mais la difficulté 
mentionnée ci-dessus, pour formuler cette condition pour 
des particules, leur a échappé. 

Enfin Wigner a proposé une nouvelle approche, où il part 
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également de la condition que la probabilité de trouver la 
particule en dehors du noyau ne doit pas être supérieure à 1. 
Cette méthode n’est malheureusement publiée qu’en 
partie [39]. 


8. Les relations de Kramers et Kronig en théorie 
des champs. — Jusqu'ici nous avons toujours formulé 
la condition de causalité à l’aide des ondes incidente et 
émergente. Cela est nécessaire, tant qu'on ne veut faire 
aucune hypothèse spécifique sur l’intérieur du milieu, 
resp. sur le centre diffuseur, Dans le cas de la lumière, 
cette méthode conséquente ne donnait lieu à aucune 
difficulté ; mais dans le cas de particules, le traitement 
présente des difficultés considérables. On peut le rendre 
plus aisé en admettant que l’état à l’intérieur de la 
région où a lieu l'interaction peut lui aussi être décrit 
à l’aide de grandeurs de champ. Com me condition de 
causalité, on peut alors exiger : le champ en R au 
temps { ne dépend pas du champ en R,Ç au temps anté- 
rieur to, Si [R — Rol > C(t — to). 

En théorie des champs, on formule cela ainsi. Suppo- 
sons que les grandeurs de champ (dans la représen- 
tation de Heisenberg) soient décrites par les opéra- 
teurs o,(R, t). Dans le cas du champ électromagnétique, 
les ©, représentent les quatre composantes du poten- 
tiel-vecteur ; dans le cas de particules chargées, on 
prend les composantes du quadrivecteur courant. La 


condition d'indépendance mentionnée plus haut 
s'exprime alors par la relation de commutation 
[ou(Ro, to); vir, t)] = Osi(R— 6)? > ct — t0)*. (32) 


Il est clair que (32) n’est pas modifiée par le passage 
à un autre système de référence. 


Cette condition peut aussi s’interpréter autrement. La 
relation de commutation (32), prise pour t = to, exprime 
que œu(r, to) et œv(r, ts) commutent tant que Ro Æ R. 
Ceci garantit la possibilité, à un instant #,, de prendre des 
valeurs quelconques pour les grandeurs de champ en tous 
les points de l’espace. La relation (32) pour #, exprime que 
la même possibilité doit exister pour chaque observateur, 
indépendamment du système de référence ou il se trouve. 
Selon cette interprétation (32) sera considérée comme une 
condition pour que les équations de champ possèdent une 
solution pour des valeurs initiales arbitrairement choisies. 
De ce point de vue, c’est donc une condition pour la consis- 
tance de la théorie. En outre, nous avons supposé impli- 
citement dans ce qui précède que la solution est possible 
par les méthodes habituelles de la théorie des champs, en 
particulier le développement en ondes planes. Cela implique 
qu'aucun terme exponentiellement croissant, du type des 
solutions anormales citées au paragraphe 4, ne péut figurer 
dans ce développement. Plus précisément : on suppose que 
les solutions qui tendent, lorsqu'on « débranche » l’inter- 
action, vers les solutions ondes planes de l’équation de 
Schrôdinger non perturbée, forment un système complet ; 
par superposition de ces solutions on peut alors obtenir- 
toutes les valeurs des grandeurs de champ au temps to (?). 


On ne peut ici qu’indiquer le traitement ultérieur de 
la condition (32). A l’aide de la théorie des champs on 
calcule l'amplitude de diffusion en avant et on obtient 
pour cette amplitude (ou du moins pour une autre 
grandeur liée avec celle-ci) une expression qui en 


() Cet énoncé diffère de l'hypothèse d’Heisenberg ($ 7) 
que les états stâtionnaires forment un système complet 
dans la région extérieure au noyau. 
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raison de (32) est holomorphe dans Z,. Il en résulte des 
relations de Kramers et Kronig qui lient la partie 
« dispersive » de l'amplitude de diffusion avec la partie 
«absorptive ». Cette méthode a été développée surtout 
par Goldberger [40], et appliquée pour la première fois 
à la diffusion des rayons + par les noyaux. Ces deux 
dernières années, c’est surtout l’application à la diffu- 
sion des mésons x par les protons qui a attiré l’atten- 
tion [41]. Jusque vers 250 MeV, le processus de diffu- 
sion peut être décrit à l’aide d’ondes s et p. Les dépha- 
sages correspondants sont notés 1, @3, Et Gi, Mia, ay, 
%33. &34 (onde p, spin isobarique 3/2, moment cinétique 
total 3/2) présente un pic de résonance vers 200 MeV 
(dans le système du laboratoire) tandis que les autres & 
varient moins rapidement avec l’énergie [42]. Cepen- 
dant il y a une difficulté : les sections efficaces obser- 
vées ne permettent pas de déterminer sans ambiguïté 
les déphasages (par exemple, les sections efficaces ne 
sont pas modifiées si on change le signe de tous les 
déphasages ; cf. la fin du $ 6). Mais les déphasages 
sont liés de façon univoque à l’amplitude de diffusion. 
En cherchant quelle solution pour les déphasages 
conduit à une amplitude de diffusion qui vérifie les 
relations de Kramers et Kronig, on peut faire un choix 
entre les différentes solutions possibles [43] (1). 

Pour obtenir ces relations de Kramers et Kronig, on 
doit intégrer la partie imaginaire de l’amplitude de 
diffusion sur toutes les énergies, y compris l’intervalle 
non physique entre — m, €? et mx €? (cf. $ 7). Le pro- 
longement analytique de l'amplitude de diffusion est 
réel dans cet intervalle, et cependant celui ci fournit 
une contribution. En effet, il y a un pôle lié à la possi- 
bilité d’un état lié virtuel du pion avec le proton. Le 
résidu de ce pôle contribue à l'intégrale de la relation 
de Kramers et Kronig, et il se trouve qu'il contient la 
constante de couplage f du proton avec le champ 
mésique. Ceci permet de déterminer la valeur de cette 
constante à partir des expériences de diffusion ; Haber- 
Schaim a trouvé ainsi f? — 0,082 + 0,015 [45]. 


Les grandeurs qui décrivent le proton dans l’hamiltonien 
concernent le « proton nu», c’est-à-dire le proton tel qu’il 
existerait s’il n’y avait pas de champ mésique. Pour décrire 
l'interaction avec le champ mésique, on ajoute à l’hamil- 
tonien un terme d'interaction, avec une « constante de 
couplage non renormalisée » f,. Par suite le champ mésique 
au voisinage du proton sera polarisé ou autrement dit le 
proton s’entoure d’un nuage de mésons virtuels. Ce qu’on 
observe dans une expérience est l'interaction de ce « proton 
habillé » avec un méson réel. La force de cette interaction 
modifiée sera décrite par une nouvelle constante, la cons- 
tante de couplage renormalisée f. C’est donc cette grandeur 
que l’on doit porter dans la relation de Kramers et Kronig 
lorsqu'on veut la comparer avec l’expérience. D'ailleurs la 
théorie de Chew et Low conduit à des résultats très ana- 
logues [46]. 


Lors de la diffusion du méson +, il peut arriver que 
le spin du proton se renverse ; ceci donne un terme 
« spin-fip » dans amplitude de diffusion. Ce terme 
contient un facteur sin 0, doncil est nul pour la direction 
en avant et n'intervient pas dans la relation de 
Kramers et Kronig: C’est pourquoi on a démontré une 


(1) Choix qui aboutit à la solution dite de Fermi, que 
d’ailleurs Bethe, de Hoffmann e. a. [44] avaient rendue 
plausible pour d’autres raisons. 
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nouvelle relation de Kramers et Kronig pour 
[df(w, 0)/d0]8 = 0 [47] dans laquelle le terme avec ren- 
versement du spin est important. En portant dans 
cette relation les données expérimentales, on a pu 
confirmer à nouveau la valeur de f?, et du même coup 
le choix entre les différentes possibilités pour les dépha- 
sages [48]. Plus généralement, on a réussi à écrire des 
relations de Kramers et Kronig pour f(w, 0) avec 
0 -£ 0. Dans ce cas ce n’est pas l’angle de diffusion 
qui doit être gardé constant dans l'intégrale sur ©, 
mais le transfert d’impulsion [49]. Les formules obte- 
nues sont compliquées et difficiles à appliquer. 

Ces succès des relations de Kramers et Kronig appli- 
quée à la diffusion pion-nucléon ont été récemment 
quelque peu perturbés par Puppi [50]. Cet auteur a 
trouvé que les données expérimentales sur la diffusion 
des x ne satisfont à la relation de Kramers et Kronig 
correspondante que si on prend pour f? la valeur 0,04 
en contradiction avec la valeur trouvée ci-dessus. Doit- 
on en conclure que la condition de causalité est violée ? 
Cela me semble prématuré parce que dans la démons- 
tration de la relation de Kramers et Kronig concernée, 
en plus de la condition de causalité on utilise aussi des 
hypothèses plus spécifiques concernant l'interaction. 
Déjà l’intervention d’une constante de couplage — 
qui provient toujours de l’hamiltonien — montre que 
les relations de Kramers et Kronig utilisées ici sont 
moins purement phénoménologiques que dans les para- 
graphes précédents. 


9. Conclusion. — La condition de causalité est-elle 
une loi inéluctable ? On ne saurait l’affirmer, car le 
contenu précis de la condition dépend de la façon dont 
elle est formulée pour une situation donnée. Pour la 
formuler, on doit définir une cause et un effet, c’est-à- 
dire qu’on doit pouvoir circonscrire un système phy- 
sique et fixer ce qui sera considéré comme une action 
de l'extérieur, et ce qui sera considéré comme la réac- 
tion du système. En macrophysique cette séparation 
peut être menée jusqu’au bout sans ambiguïté, pourvu 
que l’on puisse faire abstraction des fluctuations statis- 
tiques. Il est évident que les considérations statistiques 
interviennent inévitablement à propos de la condition 
de causalité puisqu'il faut distinguer entre le futur et 
le passé. Déjà dans l’exemple mécanique simple de 
l’oscillateur amorti du $ 3, il est essentiel que le coef- 
ficient de frottement y soit positif, ce qui rompt la 
symétrie entre les deux directions du temps. La rela- 
tion entre la mécanique statistique et les propriétés 
causales du champ électromagnétique a été expliquée 
par Wheeler et Feynman [51]. 

En microphysique, il n’apparaît pas de difficultés 
nouvelles, tant qu’on applique la condition de causa- 
lité aux seuls instruments de mesure macroscopiques. 
Par exemple dans une expérience de diffusion, la cause 
est l’émission des particules en direction de la cible ; 
l'effet est le fonctionnement du compteur qui reçoit les 
particules diffusées. Il faut seulement remarquer 
qu’une acausalité pendant un temps très court (10—?3s, 
pour l’électron de Lorentz) n’est pas en contradiction 
avec lexpérience. 

Cependant cette formulation orthodoxe est beaucoup 
trop faible pour donner des conséquences intéressantes. 
C’est pourquoi aux paragraphes 6 et 7 nous avons 
posé une condition plus forte en exigeant que la causa- 
lité soit valable également pour les fonctions d’onde x 
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resp. V elles mêmes. Cela est hautement spéculatif 
puisque les valeurs de w et de 4 (ou du moins |4|?) en 
un point donné ne peuvent pas être effectivement 
mesurées. Cependant il est conforme à l’esprit de la 
mécanique quantique de traiter w et [b[? comme des 
grandeurs observables, en sorte que les résultats obte- 
nus sont dignes de confiance. Au paragraphe 8 tou- 
tefois, il s'agissait de grandeurs de champ quanti- 
fiées +,. Celles-ci ne sont sûrement pas mesurables 
localement, parce que leurs interactions avec les autres 
champs ne peuvent être éliminées. C’est pourquoi nous 
considérerons (32) plutôt comme une condition pour 
la cohérence de la théorie des champs que comme 
l'expression d’un fait d'expérience. 

Il semble que les difficultés de la théorie des champs 
actuelle sont étroitement liées à l’existence d’acau- 
salités. Ceci est illustré par les tentatives (qui n’ont pas 
abouti) d'améliorer la théorie en introduisant une 
longueur élémentaire ; elles se ramènent toujours à 
abandonner la causalité pour les faibles distances. Il 
est possible que la nature soit réellement acausale aux 
faibles distances, c’est-à-dire qu'aucune description 
causale avec des grandeurs définies localement ne soit 
possible. Il est également possible que ce soit seule- 
ment la théorie des champs utilisée aujourd’hui qui 
ne soit pas la bonne, et qu’il existe une description 
causale avec des grandeurs de champ convenablement 
définies (probablement d’une façon étroitement liée 
aux résultats d'observations). 

Je voudrais me permettre encore quelques remarques 
d'ordre général. La démonstration des relations. de 
Kramers et Kronig n’est pas simple. En premier lieu, 
il faut formuler la condition de causalité avec précision, 


‘et se rendre compte clairement de sa signification. Il 


n’est guère satisfaisant de déduire de belles formules 
d’une hypothèse, si la justesse de l'hypothèse est moins 
évidente que celle des formules qui en résultent. En 
second lieu, la déduction mathématique exige plus de 
précision qu’il n’est habituel en physique, car dans le 
domaine des fonctions analytiques l'intuition physique 
est un guide peu sûr. Par exemple il ne sert à rien de 
savoir d’une fonction possède un prolongement holo- 
morphe dans Z,,si on ne sait rien sur le comportement 
à l'infini de ce prolongement. Autre exemple : il 
semble séduisant de déduire les relations de Kramers 
et Kronig à partir de l’hypothèse que les états station- 
naires de diffusion forment un système complet de 
fonctions à l'extérieur du noyau ; comme on l’a remar- 
qué au paragraphe 7, ceci est impossible. Il est égale- 
ment impossible de les obtenir à l’aide d’une théorie 
des perturbations, qui fournit un développement de 
amplitude de diffusion ou de la matrice $ en puis- 
sance de la constante de couplage ; car il est essentiel 
qu’un terme d’amortissement figure dans le dénomi- 
nateur, pour que le pôle du dénominateur de résonance 
s’écarte de l’axe réel vers le bas [52]. Enfin, il faut 
toujours se rappeler que si les relations de Kramers 
et Kronig constituent bien une condition à laquelle 
doivent satisfaire les grandeurs qui décrivent la diffu- 
sion, elles ne permettent nullement de déterminer 
celles-ci sans ambiguïté. Elles ne peuvent done pas 
remplacer une formule de dispersion spécifique, mais 
seulement restreindre les possibilités pour cette 
formule. 


Manuscrit reçu le 26 octobre 1960. 
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EFFETS D'EMPILEMENTS 
ET DE TEMPS-MORT 
DANS LES SÉLECTEURS D’AMPLITUDE 


Par P. Moarri, 


Fosition du problème. — Dans un sélecteur d’ampli- 
tude, le canal, de largeur Av, enregistre An impulsions 
de hauteur v, alors qu’il en arrive en réalité An, pures. 
Chaque impulsion comptée déclenche un temps mort 
de durée 7, sur le canal de niveau v ou sur plusieurs 
canaux. Le temps actif étant +, quelle est la relation 
entre An et An, ? 


Résultats théoriques. — 1) Le point d'application 
du temps-mort dans la chaîne électronique n’a aucune 
importance, lorsqu'il est unique. 

2) Le nombre de coups comptés An doit être rem- 
placé par : 

TRE ee 


Ang = Em. nolesr — 1) A(v) Av (1) 


avec : 


ñy.— An pour un sélecteur à un canal. 
KE 


D 2 Ani pour un sélecteur à Æ canaux à T. M. syn- 
chronisés. 
Ny = À pour un sélecteur couvrant tout le spectre à 


T. M. synchronisés. 
nr est le nombre de coups comptés déclenchant le 
temps mort T, et 
ES AV 
TE ee 


fe 
4 not 


ce qui permet de calculer le nombre total de coups 
vrais 29, La fonction A(v) est la fonction d’empilement, 
déjà étudiée. 


Caleuls pratiques. — Dans l’approximation usuelle, 
fol on a: 


TT 


9 
.An no 
LT = 


An ‘(2) 


T.h{v) Av. 


La correction de T. M. est rigoureuse et immédiate, 
quelle que soit la valeur de T.. 

Rappelons que la correction d’empilement n’est né- 
gligeable que si 
An 1 


h(v) € == 
() 10 r Av 


S 


Les effets de temps mort et d’empilement sont don 
assez nettement séparés. 


Conclusions. — 1) La correction de T. M. est nulle : 


pour un sélecteur total ; la forme du spectre n’est pas 


modifiée, mais seulement le nombre de coups comptés. " 
2) Elle est très faible pour un sélecteur à un canal, 
surtout si sa largeur est petite. Elle sera souvent 


négligeable si T n’est pas très grand. 


3) Elle est maxima pour un sélecteur multicanal ne « 


couvrant qu’une fraction notable du spectre, le 1/5 par 


exemple. 
4) On doit vérifier que : 


An 
il TER Try Fa Anmax = 


est indépendant de la valeur de T, Ansaz étant le 
nombre de coups comptés sans temps-mort. 

5) Le «pouvoir de résolution» du système électro- 
nique est : 0 — + + T, somme du temps actif et du 


temps-mort. C’est le temps au-dessous duquel le sys-. 


tème ne compte qu’une impulsion au maximum. 


Mesure du temps mort. On déduit facilement de la 
relation ci dessus ; 
4 À. 


nr  nrlr=0 


(3) 
Il suffit donc de 2 comptages, l’un avec, l’autre sans 
T. M., pour déterminer T. La méthode exige qu’on 
puisse supprimer le T. M., ce qui est aisé quand il 
est créé sur une voie secondaire. 
Lorsqu'il est impossible de l’annuler, on peut mesu- 


ie - 
Ga 


ED dar (ANS 
ET DE 


rer T par cette méthode en faisant varier le paramètre 


déterminant (capacité de liaison d’un monovibrateur 


par exemple), et en extrapolant le nombre maximum 
de coups comptés. 


rer 


La variation de 1 /#r étant linéaire en fonction de T 


et approximativement linéaire en fonction de C, l’extra- 
polation est assez précise. La différence des ordonnées 
représente exactement le temps mort 7. 


Lettre reçue le 22 octobre 1960. 
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